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An experiment has been performed to measure the correlation between direction of neutron 
and fragment in the fission process. The results are consistent with the isotropic evaporation 


of neutrons from the moving fragments. 


HE simplest picture of the emission process 
of fission neutrons is that these are 
evaporated from the moving fission fragments. 
If this is true, and if the evaporation velocity 
is comparable with that of the fragments, then 
one would expect the neutrons, when measured 
in the laboratory system, to be emitted pre- 
dominantly in the direction of the fragments. 
The present experiment was designed to detect 
and measure the directional correlation of the 
fission fragments and the emitted neutrons. 
Figures 1 and 2 show the experimental ar- 
rangement. The direction of the fission fragments 


was determined by the thin collimator placed 
over the Uss35 foil which was one mg/cm* thick. 
The collimator consisted of a grid of closely 
spaced js in. holes drilled in a $ in. thick block of 
aluminum. The holes were drilled over a square 
area one inch to the side which corresponded to 
the area of the foil. The fission fragments pene- 
trating the collimator were detected by the 
ionization chamber formed from the collimator 
and the ion collection plate } in. away. It was 
possible to rotate this system with respect to a 
neutron counter which consisted of a bundle of 
ten proportional counters 1} inches long and 


TABLE I, Counting rate of the numbers of neutrons emitted at angles 0 and #/2 in the laboratory system. 




















Accidental True coin- 
Recoils Fissions Coincidences Coincidences counts cidence counts 
Time _— -———_ _ —- per 2.56- 10 per 2.56- 10+ per 2.56-10* 

Angle min. 64 1024 4 fissions fissions fissions 

0° 23 472 1941 1.0 2.11 0.29 1,82 

x/2 14 273 1093 6.7 0.61 0.27 0.34 

ba 5 91 403 8.5 2.11 0.26 1.85 
3x/2 6 108 573 3.5 0.61 0.25 0.36 


Av. of number of neutrons emitted per fission at 0° and x; n(0) = (1.8347%)/2.56- 10*. 
Av. of number of neutrons emitted per fission at #/2 and 34/2; n(x/2) = (0.34416%)/2.56- 10° 
n(4x)/n(0) =0.18+0.03. 











*This document is based on work 


performed in 1945 at Los Alamos Scientific Laboratory of the University of 


California under Contract No. W-7405-eng-36 for the Manhattan Project, and the information contained therein will ap- 
pear in Division V of the Manhattan Project Technical Series as part of the contribution of the Los Alamos Laboratory. 
** Now at Cornell University. 
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Fic. 1. Plan view of counter arrangement. 


# in. in diameter. The three-mil wires along the 
axis of all counters were connected in parallel to 
a fast amplifier (Model 500). The inside of each 
cylinder was coated with paraffin. The sensitive 
area of the counter was very nearly a square 
one-inch to the side. The front of the neutron 
counter was 1j in. from the center of the Uses; 
film. The geometry was such that if the neutrons 
were emitted exactly in the direction of the 
fission fragments, then no neutrons would be 
counted by the proportional counter at angles 
greater than 45° with respect to the direction of 
the fragments. Argon at 60 cm Hg plus CO, at 
one cm Hg was used to fill the chamber and was 
thus the counting gas for both the ionization 
chamber and. the proportional counter. A nega- 
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tive potential of 1000 volts was applied to the 
collimator and the cylinders so that electrons 
were collected in each case. A typical number. 
bias curve for the proportional counter is shown 
in Fig. 3. 

The fission pulses from the ionization col. 
lection plate and the neutron pulses from the wire 
were amplified separately and then led to a coin- 
cidence counter. The number of fissions, the 
number of neutrons, and the number of coin- 
cidences were all counted simultaneously when 
the foil was irradiated by the nearly pure 
thermal neutron flux of the water boiler. The 
resolving time of the coincidence system was 
0.254 sec. as determined by placing a strong 
source of fast neutrons at various distances from 
the proportional counter and then measuring the 
accidental counts. 

The first run was made to determine the 
counting rate of the number of neutrons per 
fission in the direction of the fission fragments 
(0°), to the counting rate of neutrons per fission 
emitted perpendicular (in the laboratory system) 
to the fragment direction. The fission counter 
bias was set as low as possible so that all frag- 
ments were registered. Table | summarizes the 
data. The neutron counting rate observed per- 
pendicular to the fragment direction was about 
one-fifth of that observed in the direction of the 
fragments. 

The second run was made to try to determine 
the effect of the energy of the fission fragment 
counted both on the number of neutrons emitted 
in the direction of the fragment (0°) and the 
opposite direction (x). The energy of the frag- 
ment was varied by changing the fission counter 
bias. A typical number-bias curve is shown in 
Fig. 4. The gain of the amplifier was different 


TABLE II. Influence of the fission fragment energy on the numbers of neutrons emitted at angles 0 and -. 


















































Accidental 
Fission Recoils Fissions Coincidences Coincidences counts per True coincidence 
counter Time per 2.56- 10+ 2.56- 10+ gy = 
Angle bias in volts min. 64 1024 4 fissions fissions « 2.56-10 fissions 
0 0 10 250 796 15 1.88 0.35 1.53+0.20 
0 0 10 240 742 12 1.62 0.34 1.28+0.24 
0 40 10 239 356 7.75 2.18 0.34 1.84+0.33 
0 50 10 238 106 2 1.9 0.34 1.6 +0.6 
ca 0 10 218 682 10.5 1.55 0.30 1.55+0.24 
40 10 230 318 4.0 1.26 0.33 0.93 X0.25 
50 10 239 55 0 0 0.35 —+100% 
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FISSION NEUTRONS 


TaBLe III. Angular distribution of the neutrons emitted from fission fragments. 

















Recoils Fissions Coincidences Coincidences Accidental True coincidences 
Time ————— per 2.56- 10+ counts per counts per 
Angle min. 64 1024 4 fissions 2.56- 10* fissions 2.56- 10 fissions 
0 10 331 729 16.25 2.23 0.47 1.76+0.28 
1/16 12 501 883 14 1.58 0.59 0.99+0.21 
1/8 12 417 868 7 0.81 0.49 0.32+0.15 
3x/16 11 450 780 6 0.77 0.57 0.20+0.16 
r/4 13 576 910 6.25 0.68 0.61 0.07+0.13 
3n/8 12 525 $24 5 0.61 0.66 —0,05+0.13 
7/16 13 595 905 9.75 1.07 0.64 0.43+0.17 
© 12 556 831 13 1.56 0.65 0.91+0.22 
g 10 487 775 11.75 1.51 0.68 0.83+40.22 
3/2 10 439 809 5 0.62 0.62 0.00+0,12 
le 12 518 779 12.25 1.58 0.60 0.98+0.23 








at that time. The data are summarized in 
Table II. The statistics are poor, as can be seen 
from the large standard errors given in the last 
column. The energy of the fragment has no 
evident effect on the number of neutrons per 
fission emitted in the forward direction. There 
is a slight indication that fewer neutrons are 
emitted from the low energy fragment because 
there seem to be fewer recoils in the backward 
direction at the higher biases. This effect may be 
instrumental. 

A third run was made to obtain a curve of the 
numbers of neutrons emitted per fission at 
various angles. Table III shows the data, and the 
standard errors listed in the last column indicate 
that the results are too rough to plot. Further- 
more, the sensitivity of the neutron counter 
seemed to change during the course of the run. 





Fic. 2. View looking down on the chamber with enclosing 
shell removed. The rotating fission collimator and counter 
are to the right; the neutron counter is to the left; and the 
preamplifier is below. 





However, anything really unusual might have 
been noticed even with such rough data. 

The first run was made under good conditions, 
so let us see what conclusions can be drawn from 
it about the process of neutron emission. Assume 
the simplified picture that the fission fragments 
all have the same velocity and that the neutrons 
are evaporated isotropically with respect to the 
moving fragment and with a uniform velocity. 

The ratio of the number of neutrons emitted 
in the forward direction to those emitted at 90° 
in the laboratory system can be calculated to be 


N($x)/N(0) =r(r?—1)8/(r+1)°, (1) 


where r =v,/v;, and v, is the speed of the emitted 
neutron with respect to the fragment, and », is 
the speed of the fragments. 

Now the proportional counter is sensitive to 
the energy of the neutrons. The speed of a 
neutron emitted at 0° is (v,+v,), while that of a 
neutron emitted at an angle x/2 is (v,?+v/)!. 
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Fic. 3. A typical number-bias curve for the recoil pro- 
portional counter. : 
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Fig. 4. A typical number-bias curve for the fission counter. 
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Fic. 5. Energy response of the recoil proportional counters. 


The energy sensitivity of a neutron counter of 
the same dimensions and gas filling was taken 
at a different time, using the lithium (p, 2) 
neutrons obtained from a Van de Graaff genera- 
tor. The curve is shown in Fig. 5, and for our 
purpose the sensitivity can be approximated 
roughly by K’(E,)!. If we make this correction 


to Eq. (1) we get 
n(x) /n(0) =r(r —1)/(r +1)? (2) 


where n(}x)/n(0) is the ratio of the numbers of 
neutrons counted at angles 7/2 and 0. In Fig. 6 
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Fic. 6. Theoretical ratio of neutron counts at right angles 
to fission direction to neutron counts in the direction of 
the fission fragment as a function of r, the ratio of the 
speed of the emitted neutron with respect to the fragment 
to the speed of the fragment. 


is plotted ($2) /n(0) as a function of r, and we 
can read a value of r=1.8+0.2 corresponding to 
the experimental value m(}2)/n(0) =0. 18+03 
If the average energy of a fission fragment is 70 
Mev, and the average atomic weight is 100, the 
velocity v; is the same as that of a neutron of 0,7 
Mev; hence the energy of evaporation would 
be 0.7 X (1.8)? = 2 Mev. Had we assumed that the 
detector was insensitive to energy the evapora- 
tion energy would be about one Mev. The 
crudeness of the experiment and its interpreta- 
tion are apparent. 

I had hoped to split the neutron counter into 
two equal parts and count coincidences between 
these neutron counters when they were on 
opposite sides of a piece of uranium exposed to 
thermal neutrons and when the two counters 
were on the same side. This would give infor- 
mation about the mechanism of emission of 
fission neutrons as well as of the quantity 
(v?)—v)m?, where v is the number of neutrons 
per fission. However, the very brief time avail- 
able for the experiment did not permit this 
extension. 
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A Thermal Neutron Velocity Selector and Its Application to the Measurement of the 
Cross Section of Boron 


Argonne National Laboratory,* University of Chicago, Chicago,** Illinois 


(Received April 25, 1947) 





A mechanical velocity selector for the study of moncchromatic neutrons in the range of 


energies below 0.3 ev is described. 


The instrument has been applied to the measurement of the cross section of boron, which 


is found to be 703 X 10-* cm? for neutrons of 2200 meters per second velocity. 








INTRODUCTION 


LOW neutrons emerging from various mod- 

erators with different geometries usually 
have average velocities comparable, but by no 
means equal, to the thermal agitation velocity. 
Large differences, both positive and negative, 
are observed depending on the nature and the 
geometry of the moderating substance. This 
phenomenon has been observed by various 
experimenters.'~* 

In this paper we have collected some typical 
examples of the variations of average velocity of 
slow neutrons using different moderators as 
indicated by changes in the apparent cross 
section of boron. Since boron is often used as a 
standard substance in slow neutron measure- 
ments, its cross section has been determined also 
by use of monochromatic neutrons obtained with 
a velocity selector of new design operated in 
connection with the thermal column of the 
Argonne graphite pile. 

The observed temperatures of the neutrons 
emitted from the various moderators and ar- 
rangements of moderators appear to be in 
accordance with the individual arrangements 
employed. Within the experimental errors of the 
method, the cross section of boron varies as the 
1/v law, and the measured cross section is 


*The information contained in this document will 
wet in Division IV of the MPTS as part of the con- 
tribution of the Argonne National Laboratory. 

** All three authors now at Institute for Nuclear Studies, 
University of Chicago. 

1J. Rainwater and W. W. Havens, Jr., Phys. Rev. 70, 
136 (1946). 

*W. W. Havens, Jr. and J. Rainwater, Phys. Rev. 70, 
154 (1946). 

3 J. H. Manley, L. J. Haworth, and E. A. Luebke, Phys. 
Rev. 69, 405 (1946). 

*R. F. Bacher, C. P. Baker, and B. D. McDaniel, Phys. 
Rev. 69, 443 (1946). 
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703 X 10-* cm? per atom per neutrons of velocity 
2200 meters per second. 


TEMPERATURES OF NEUTRONS FROM 
VARIOUS SOURCES 


With the thermal purification column of the 
graphite pile at the Argonne Laboratory as a 
primary neutron source, a number of measure- 
ments were made of the cross section of boron. 
In all cases the detector was a proportional - 
counter filled with BF; gas. By the use of cad- 
mium diaphragms a neutron beam was obtained 
with small angular dispersion. 

The absorber and detector in these experi- 
ments were both boron and consequently both 
obeyed the 1/v law of neutron absorption. It was 
possible, therefore, to use the correction method 
given by Bethe' to calculate the cross section of 
boron for mono-energic neutrons of energy 
kT where T is the absolute temperature of the 
Maxwellian distribution emitted from the source. 
Since the cross section of 2200 meters per second 
neutrons (kT at 293°K) is known, one can then 
determine the effective temperature of the 
neutron beam. It must be understood that these 
effective temperatures are based on the assump- 
tion that the neutron beam is Maxwellian in 
velocity distribution. This is certainly not strictly 
true for most sources employed. 

The results of these experiments are given in 
Table |. It is quite clear from an inspection of 
the table that the effective temperature of the 
neutron beam depends strongly on the source of 
neutrons. During these experiments the tem- 
perature of the thermal column was in the 
neighborhood of 30°C or 303°K. 


‘ 5H. A. Bethe, Rev. Mod. Phys. 9, 134 (1937). 
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TABLE I. 
Source of neutrons Absorber pened nrnany ta A py 
1. Beam from surface of thermal column Gaseous BF; op=855X10-*4 198 
2. Beam passed through a 3.7-cm slab of paraffin Gaseous BF; 598 x 10-4 408 
3. Beam passed through 7.6c m. of heavy water at Gaseous BF; Corrected to 20.4°C 288 
33.7°C in a container 18-in. diam. op=710X10-* 
4. Beam passed through a 22 cm. column of graphite Pyrex plate calibrated in 2800 x 10-* 18.4 
10 cm square velocity selector 
5. Beam from hole in thermal column 125 cm deep, Gaseous BF; 701 x 10-*4 293 
10 cm square 
6. Beam from a “black hole” in thermal column, a Gaseous BF; 755 x 10-*4 255 


hole 10 cm X10 cm X22 cm high connected to 
surface of thermal column by a 42-cm tube of 
cadmium of internal diam. 2.5 cm 








The source arrangement given under case 1 
produces low temperature neutrons because of 
the filtering action of the graphite in the pile and 
thermal column.* Very slow neutrons whose 
de Broglie wave-lengths are longer than peri- 
odicities encountered in the graphite crystals are 
scattered very little and can penetrate to the 
surface of the column more easily than the 
faster neutrons. In case 2 the slower neutrons are 
removed preferentially because both the ab- 
sorption and scattering cross sections of hydrogen 
are larger and also because scattering in the 
forward direction is preferred at higher energy. 
Heavy water (case 3) acts somewhat in the same 
way because also for deuterium compounds the 
scattering cross section and the coherence of 
successive free paths vary with the energy in 
the same direction as for hydrogen compounds. 
Therefore, the effective temperature of the 
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Fic. 1. Cross section of the shutter of the velocity selector. 


* H. L. Anderson, E. Fermi, and L. Marshall, Phys. Rev. 
70, 815 (1946). 
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neutrons is raised from the initial 198°K to 
288°K. The fact that this last temperature is 
quite close to the actual temperature of the 
heavy water probably is coincidental. In case 4 
the filtering effect of the graphite is shown very 
strongly. Most neutrons that are scattered are 
removed from the beam and the graphite column 
is so long that almost none of the warm neutrons 
can travel the whole distance without being 
scattered. Case 5 gives a rather good approxima- 
tion of the temperature of the source. The neu- 
trons in the beam from the deep hole should be 
a fair sample of the neutrons present at the 
bottom of the hole. Essentially it is a case of 
blackbody radiation from a hole in the wall of a 
furnace. Case 6 was expected to give a good tem- 
perature value, but failed to do so, probably 
because the hole was not deep enough. 


VELOCITY SELECTOR 


The velocity selector makes use of a rotating 
shutter to interrupt the beam of neutrons from 
the thermal column of the pile. The shutter was 
constructed by inserting a multiple sandwich of 
0.004-inch to 0.008-inch cadmium foils and 
#s-inch aluminum sheet tightly into a steel 
cylinder about 1} inch in diameter with walls 
#s inch thick. The shutter was mounted in ball 
bearings on a heavy steel base plate and was belt 
and pulley driven by a Dumore grinder motor. 
Maximum rotational speeds of 15,000 revolu- 
tions per minute were possible. It was con- 
structed in the shops of the Metallurgical 
Laboratory under the direction of Mr. T. J. 
O’Donnell who is responsible for its mechanical 
design. 
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Across section of the shutter is shown in Fig. 1. 
From the thickness of the aluminum spacers 
between the cadmium foils, and from the dimen- 
sions of the shutter, one would estimate that no 
neutrons from a parallel beam would be able to 
get through when the shutter was more than 
1.2° from its full open position. In the experi- 
mental arrangement used, it was impossible to 
use a strictly parallel beam of neutrons. The 
collimators actually used allowed a maximum 
divergence of neutron direction in the beam of 
approximately 3°. Consequently, one would 
expect the shutter to be completely closed during 
each 180° of rotation except for an interval of 
3°+2X1.2°=5.4°. Actually it was found that 
the counters indicated background intensity 
except when the shutter was in a 6° interval. 

Through one end of the shutter was inserted 
asteel rod with its axis perpendicular to the axis 
of the shutter and with a minor surface ground 
and polished perpendicular to its axis at each 
end. Light from a projection lamp and lens 
system was reflected from these surfaces onto 
two photo-cells so placed that each photo-cell 
was illuminated twice during each revolution. 
One of the photo-cells was used with an amplifier 
and scaling circuit as a revolution counter. The 
other, adjustable and calibrated as to angulars 
position, was connected to an electronic switch 
circuit which allowed pulses from the propor- 
tional counter to be recorded only when the 
photo-cell was illuminated. 

BF; filled proportional counters were used as 
the neutron detector. A nest of four was con- 
nected in parallel and mounted at a distance of 
146 cm from the shutter. A thick shield of wood, 
iron, and paraffin was placed between the 
counters and the pile to compensate somewhat 
for the fact that the top shield of the graphite 
pile was not so thick as might be desired. A hole 
in this shield allowed neutrons from the shutter 
to reach the counters. 

The neutron beam between the shutter and 
the counters was collimated to make sure that 
no slow neutrons from sources other than the 
shutter could enter the counters. Slow neutrons 
reflected from the walls and roof of the building 
were eliminated by protecting the sides and back 
of the counters with a }-inch thick layer of boron 
carbide. 





The shutter and an improved velocity selector 
arrangement are to be more fully described in a 
paper by Brill and Lichtenberger. 


DETERMINATION OF BORON CROSS SECTION 
FOR NEUTRONS OF KNOWN VELOCITY 


The cross section of pure BF; at several dif- 
ferent pressures was measured for neutrons from 
the thermal velocity selector for velocities 
ranging from 1700 to 5000 meters per second. 
Within the experimental accuracy of the method 
the cross section of boron varied according to the 
1/v law. After corrections for scattering were 
made, the average cross section of boron for 
neutrons of 2200 meters per second velocity was 
699 x 10-** cm?/atom. 2200 m/sec. is the velocity 
of a neutron of energy kT where T is 293°K. 

In order to verify this value a similar measure- 
ment was made with a different boron compound 
as absorber. Na2B,O; was ignited at about 400° 
and dissolved in heavy water. The solution was 
enclosed in a thin-walled aluminum cell, and a 
second cell of identical wall thickness, which 
contained an amount of heavy water equal to 
that in the solution, was prepared. The trans- 
missions of these two absorbers for neutrons 
from the velocity selector were measured, and 
the value of the boron cross section for 2200 
m/sec. neutrons was found to be 700 X10~ cm? 
corrected for scattering. 

In good agreement with these values was the 
cross section as calculated from measurements at 
the indium resonance energy.’ Transmission 
measurements were made using a collimated 
beam of neutrons from the interior of the 
graphite pile of the Argonne Laboratory. The 
indium foil detectors were protected from 
thermal neutron activation by thick cadmium 
covers. Background measurements were made 














TaBLeE II. 
Measurement ox7(B) at 293°K 
Na2B,O;— D,O Velocity selector 700 < 107*4 cm? 
BF; Velocity selector 699 x 10-** 
BF; In resonance 710x 10-** 
Average 703 X 10-** cm? 








7 J. Marshall, Phys. Rev. 70, 107 (1946). 
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by use of an indium filter. Thus the measure- 
ments were limited in more than one way to 
neutrons absorbed strongly by indium. 

BF; gas in a steel cylinder was interposed in 
the collimated beam. The BF; was highly purified 
(the same gas as used in the thermal neutron 
transmission experiments described above). The 
transmission of the steel container filled with 
BF, at 44 and 68 lb/in.? was compared with the 
transmission of the empty container. The 
density of gas used was determined by weighing 
the cylinder. The pressures used and the length 
of the cylinder (30 cm) were such that the trans- 
missions were in an accurately determinable 
range (approximately a 3 transmission for the 
68-lb sample). 


HALL AND P. G. KOONTZ 


The total cross section of BF; for indiym 
resonance neutrons was measured as 107,] 
x 10-* cm*/atom. Assuming 


Oscattering( F) = 3.7 xX 10-*4 cm?, . 
Oscattering(B) =2X10-* cm’, ‘ 
indium resonance energy = 1.44 ev, 


the boron absorption cross section for neutrons 
at velocity 2200 m/sec. is 710 X 10-** cm?/atom, 

The results of the three measurements are 
given in Table II. 

This report is based on work done at the 
Argonne National Laboratory, the University 
of Chicago, under the auspices of the Manhattan 
District, U.S. Corps of Engineers, War Depart- 
ment. 
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Scattering of Fast Neutrons by Helium* 


T. A. HaLtt anp P. G. Koonrtzt 
University of California, Los Alamos Scientific Laboratory, Santa Fe, New Mexico 


(Received March 28, 1947) 


Measurements have been made of the angular distribution of helium recoils for incident 
neutrons of energies from 0.6 to 1.6 Mev. The distribution curves also permit estimates of the 
total cross section in this energy range. The results confirm the existence of a cross-section 
peak of about 6.8 X 10-* cm* around 1 Mev and indicate, under the assumption of s- and p-wave 
scattering only, that the peak is double. But preliminary attempts to fit the data to Bloch’s 
detailed theory of s- and p-scattering with a split p-level have not been successful, and the 
sign of the postulated splitting is not established. 


INTRODUCTION 


SCATTERING theory for helium must 

satisfy quantitatively two previous sets of 
data as well as that reported here. Barschall and 
Kanner! obtained recoil distribution curves at 
neutron energies of 2.5 and 3.1 Mev, and Staub 
and Tatel* obtained the backward scattering 
cross section as a function of energy around 1 
Mev, finding a peak (possibly double). The 
Barschall-Kanner curves are strongly anisotropic 


. ~ * This pa r will appear in Division V of the Manhattan 
Project ‘Technical Satie a0 part of the contribution of the 
Los Alamos Laboratory of the meng of California. 

t Now at the University of Chicago, . icago, Illinois. 

t Now at Iowa State College, Ames, Iowa. 

1H. H. Barschall and M. H. Kanner, Phys. Rev. 58, 
590 (1940). 

2H. Staub and H. Tatel, Phys. Rev. 58, 820 (1940). 


with a preponderance of forward scattering. 
Wheeler and Barschall’ have shown that the data 
at 2.5 Mev can be fitted by the assumption of 
strong spin-orbit coupling, a p-wave resonance 
around 2.5 Mev, and a weak addition of d-wave 
of a size not unreasonably large for 2.5 Mev. But 
this data at 2.5 and 3.1 Mev cannot be matched 
with a simple theory involving only resonances 
around 1 Mev: an s-wave resonance would be 
isotropic and higher resonances at 1 Mev must 
fall off to insignificance at 2.5 and 3.1 Mev. 
Thus the two sets of data provide essentially two 
different problems, at any rate for analysis in 
terms: of resonances. The Bloch formula for 


3 J. A. Wheeler and H. H. Barschall, Phys. Rev. 58, 682 
(1940). 
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resonant p-wave scattering’ used below can 
satisfy only one set, and we only consider the 
fitting of this formula to the data near 1 Mev 
of Staub and Tatel and of this paper. The main 
contribution from this paper is a group of recoil 
energy distribution curves for neutron energies 
from 0.6 to 1.6 Mev. While the authors have had 
opportunity to make only a preliminary theo- 
retical analysis, it has seemed worth while to 
present these new data. 

Earlier observations of scattering near 1 Mev 
were made by Gaerttner, Pardue, and Streib, 
who found a total cross-section peak around 0.8 
Mev, and by Staub and Stephens,® who found a 
total cross-section peak near 1 Mev. Taken 
together, these studies indicated a double peak, 
but they were not sufficiently detailed to permit 
a quantitative check with theory on this score. 
The question of the double peak has been inves- 
tigated in its theoretical aspects by Dancoff,’? who 
concluded that a Thomas relativistic spin-orbit 
coupling should produce an inverted doublet 
(with size of splitting essentially unknown) while 
the tensor spin-orbit interaction of mesotron 
theory could be expected to produce a normal 
doublet with splitting. of the order of 100 kev. 
None of the experimental data has determined 
whether the splitting is normal or inverted, nor 
is the size of the splitting established beyond the 
fact that it cannot be more than about 600 kev. 


EXPERIMENTAL PROCEDURE 


The measurements were made at the larger 
Van de Graaff generator loaned to the Los 
Alamos Laboratory by the University of Wis- 
consin. The lithium p—mn reaction provided the 
neutrons. A cylindrical proportional counter,® 
with guard rings at the end to define the counting 
volume, was filled with a helium-argon mixture 
and placed in the neutron flux; the helium recoil 
pulses were fed through a preamplifier and 
amplifier of rise-time } microsecond and RC 
decay time 20 microseconds into a discriminator 


‘F. Bloch, Phys. Rev. 58, 829 (1940). 

*E. R. Gaerttner, L. A. Pardue, and J. F. Streib, Phys. 
Rev. 56, 856 (1939). 
ss) Staub and W. E. Stephens, Phys. Rev. 55, 131 

™S. M. Dancoff, Phys. Rev. 58, 326 (1940). 

*P. G. Koontz and T. A. Hall, Rev. Sci. Inst. (to be 
published). 
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Fic. 1. Experimental arrangement for curves of energy 
distribution of recoils. 


which counted simultaneously all pulses of a 
size above a set level (‘‘integral count’’) and all 
pulses of a size within a set channel (“‘differential 
count’’). Both settings were adjustable, and by 
moving the channel along the recoil energy scale 
while keeping its width constant, one obtained 
the energy distribution of the helium recoils. 
(The recoil energy was simply proportional to 
the pulse height.) Two kinds of data were taken: 
first, detailed curves of the energy distribution 
of the helium recoils for a given neutron energy 
(Figs. 2-6), and secondly, brief studies of the 
high energy end of these curves at successive 
neutron energies (Fig. 7). In each case, total 
counts above a certain recoil energy (determined 
by the need to eliminate y-background) were 
recorded, and neutron flux was measured. Also 
proton current in the Van de Graaff generator 
was recorded during the runs. 

In the case of the detailed curves the experi- 
mental arrangement was as follows: The scat- 
tering chamber was placed with its end about 
9 inches from the lithium target, and its axis 
pointing to the target. (This puts the beginning 
of the counting volume about 11 in. from the 
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Fic. 3. Equation of parabolic fit: 
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target.) The chamber axis made an angle of 30° 
with the target tube axis. Also off at an angle 
of 30°, and about 9 cm from the lithium target, 
was the flux monitor which was a U*® fission 
foil placed in an ionization chamber. This ar- 
rangement is sketched in Fig. 1. The Van de 
Graaff generator was then set to give the desired 
neutron energy at 30°. The distance of the 
chambers from the target remained about as 
given throughout the runs, but was measured 
after each run to secure accurate flux measure- 
ment. 

Data were taken with two different gas fillings, 
and under operation both as an ion chamber and 
as a proportional counter with gas multiplication 
of about 14. The two fillings were (3 atmosphere 
of helium+2} atmospheres of argon), and (4 
atmosphere of helium+1 atmosphere of argon), 
the lower pressure being used to permit recording 
of slower recoils without interference from 
y-background. The collecting potential was 1450 
volts for use as an ion chamber and 2500 for use 
as a proportional counter with the higher 
pressure, and 1500 volts for use as a proportional 
counter with the lower pressure. At the higher 
pressure an a-particle of 0.5 Mev has a range of 
0.13 cm, and an a-particle of 1 Mev a range of 
0.22 cm; the corresponding figures for the lower 
pressure are 0.28 and 0.48 cm. 

In the case of the curves covering the high 
energy recoils (Fig. 7), the arrangement of the 
helium chamber was the same, but the flux 
monitor used was a fission detector almost insen- 
sitive to energy, placed about 6 feet from the 


target. 
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EXPERIMENTAL RESULTS 


The results are summarized in Figs. 2-6, Jp 
clarifying their meaning, the following facts 
should be noted: 

The curves obtained are almost the same for 
ion chamber and proportional counter operation, 
This was also the case in previous testing of this 
counter using the nitrogen n—p reaction. In the 
ion chamber case, because of the small rise and 
decay time of the amplifier and the low mobility 
of positive ions, the pulse is due mostly to the 
collection of the electrons, which travel over dif. 
ferent distances from their point of origin to the 
central wire. Nevertheless, since ion chamber and 
proportional counter curves are quite similar, 
we may assume that the pulse is proportional 
to the energy of the recoil in each case® (except 
for wall effect, irregularity of multiplication, 
etc.). Hence the abscissas in Figs. 2-6 may algo 
be considered to represent the recoil energy, 
Furthermore, writing 


E=neutron energy 
E,, =maximum energy of recoil 
E,=energy of recoil 
¢=angle of recoil 
6=scattering angle of neutron, center of gravity 
system 
dw =2zx sinéd@ 
da/dw=cross section per unit solid angle in the center of 
gravity system, 
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Fic. 4. Equation of parabolic fit: 
y = 2225 —92.77x+1.312x*. 


* The effect of “electron collection" is described ~~ 
in reference 8. It is discussed in greater detail in the fo 
coming “Los Alamos Encyclopedia” but the necessary 
— were developed before the Manhattan Project 
gan. 
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and 
En=4MaM,/(Mat+M,)’, 

from which it follows that the curve of recoil 
distribution vs. recoil energy has the same shape 
as the curve of do/dw vs. cosé, with @=0, 
cosé=1 at zero recoil energy and 6=7z, cos@= —1 
at maximum recoil energy. Thus Figs. 2-6 give 
directly the differential cross section in the 
center of gravity system. 

For each of these curves a check was made on 
y-background, runs were made at different bias 
settings with a quartz shield blocking the lithium 
target from the proton beam. This procedure 
showed the effect of y-rays from the bulk of the 
Van de Graaff generator, but eliminated y-rays 
normally produced beyond the quartz shield in 
the target tube and lithium target. However, 
other measurements'® have shown that at a 
distance of about 2 inches from the target, y-rays 
from a target and y-rays from the generator bulk 
are of about equal strength, and those from the 
target fall off as 1/r? with distance from the 


' target; so that the y-effect cannot be significantly 


larger than that observed in the procedure above. 
All data plotted were taken at points where the 
background is negligibly small. 

The theory indicates (see below) that the dif- 
ferential curves should be parabolas. Parabolas 
have been fitted and are drawn in Figs. 2-6. At 
the high energy end of the curves, the experi- 
mental curve is broadened and decreased in 
height by the width of the channel of the dis- 
criminator; the parabolic curve has been ex- 
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” These measurements were made by R. Taschek in con- 
nection with another experiment. 
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Fic. 6. Equation of parabolic fit: 
y = 1206 — 26.26x+-0.1788x*. 


tended to a point where the area under it is 
about the same as that under the experimental 
curve. A further requirement was that the end 
point of the parabolic .curve, corrected for 
channel width, should agree in height with the 
experimental peak. This requirement was met 
in each case to within 15 percent or better of the 
experimental peak height. In this way the ab- 
scissa of the high energy end can be determined 
to about 1 percent of the pulse height. The 
abscissa scale in Figs. 2—6 is arbitrary. The end 
of the parabola should be taken to correspond 
to the recoil of maximum energy, which is 
[4M,M./(M,+M,)?]=16/25 of the energy of 
the incident neutron. 

The ordinate scale in Figs. 2-6 and Fig. 7 is 
also “arbitrary,”’ being the actual number of 
counts observed. Most points represent the 
average of two points, so that twice as many 
counts are involved. When two or three runs 
were made on one point, repetition was on the 
whole well within statistical error. The “curve 
peaks”’ of Fig. 8 are also plotted on an arbitrary 
ordinate scale. These peaks are the same as those 
recorded in Fig. 7. 

Figure 9 is a curve of the total scattering cross 
section vs. energy. The numbers plotted were 
obtained in two steps: (1) calculation of the 
cross section down to the integral bias, which 
depends solely “on measured quantities, and (2) 
multiplication by a factor obtained by extra- 
polating the parabolic fit to zero recoil energy. 
(This factor is simply the ratio of areas under 
the full curve and under the curve down to the 
integral cut-off.) 

The first step involved three quantities: 
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Fic. 7. Backward scattering differential curves. 


neutron flux, number of helium atoms present, 
and integral count. The flux measurements have 
a probable error of about 5 percent, owing to the 
fact that only about 3000 fission monitor counts 
were obtained per run, and to uncertainties in 
relative distances of chamber and monitor from 
the target. The number of helium atoms present 
should be quite well known, the pressure having 
been measured by a mercury manometer, the 
temperature being known to a few degrees, and 
the active volume having been measured to about 
3 percent. The error in the number of integral 
counts is negligible, as one obtains about 15,000 
integral counts with each point on the curve. 
An associated error arises in the determination 
of the energy corresponding to the integral bias. 
This error is caused by the uncertainty in locating 
the exact position of the maximum recoil energy 
on the graph. This uncertainty amounts to about 
1} percent in the total cross section. Thus the 
probable error with perfect extrapolation would 
be 6 percent. 
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Fic. 8. Evidence of increased anisotropy at peak. 


TABLE I. Total cross section vs. neutron energy, 

















E(Mev) a(10-* cm?) 

0.6 2.48 

0.8 4.91, 4.63, 4.28 
1.0 6.42, 6.75, 6.16 
1.1 6.35 

1.25 6.96, 6.70 

1.35 5.66, 5.22 

1.6 4.60 








The extrapolation itself is uncertain. It is 
poorest at lower neutron energies, where the 
area under the extrapolated curve is a larger 
fraction of the total, and where the curve cannot 
be fitted so well because its shape is not so well 
determined. The integral bias is always set as 
low as possible without y-background inter- 
ference; at the higher gas pressures this permits 
settings at about 0.22-Mev recoil energy and at 
the lower gas pressure settings of about 0.13 Mev. 
Thus the integral bias energy ranges from about 
4 to 4 of the maximum recoil enérgy in the 
curves presented. While it is not certain that the 
curves should be simply parabolic, it is very 
unlikely that they will deviate much from a 
parabolic shape. (This would require strong 
d-wave or higher order scattering, which is very 
unlikely for energies as low as 2 Mev where the 
neutron wave-length in the center of gravity 
system is more than 3X10-" cm.) The extra- 
polated area is probably correct to 25 percent or 
better (as a guess), so that the errors arising 
from extrapolation might be expected to range 
from 4 to 8 percent. This leaves the final probable 
error on the cross sections around 10 percent. 
This figure seems reasonable considering the 
spread in Fig. 9 at those energies for which more 
than one run was taken. -Fig. 9 then gives the 
energy dependence of the cross section. The same 
points are tabulated in Table I. 


INTERPRETATION AND COMPARISON WITH 
THEORY 

The Bloch theory is general, introducing few 
assumptions. The scattering near 1 Mev is built 
from s- and resonant p-waves. A_ spin-orbit 
interaction is allowed for; it is assumed that 
there are two resonant compound p-levels, of 
energies E12 and E3,2, $ and 3 being the total 
angular momenta of the two levels. The only 
assumptions introduced are that all other levels 
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of the compound nucleus do not affect the scat- 
tering of the neutrons at the energies used, that 
the spin-orbit forces are small enough so that 
the phases of the scattered waves 6, and the 








r 1 r 


half-width [ are about the same for the two 
levels, and that the phases 4,, /#0, for non- 
resonant components, are negligibly small. 

The Bloch formula reads 





~ =x [sins +e 
dw 


where do/dw is the cross section per unit solid 
angle in the center of gravity system, X=)/2z, 
\ is the neutron wave-length in the center of 
gravity system, 5o is the phase shift of the s-scat- 


‘tered wave, I is the half-width of the He levels, 


E is the energy of the incident neutron, and 
@is the neutron scattering angle in the center of 
gravity system. The quantities T and 49 are 
dependent on the incident neutron energy, the 
dependence of 59 not being well-known apart 
from the fact that it must be slowly varying. 
The quantities T, 50, Es, and Ejy2 are the 
parameters left adjustable in fitting the data, 
for we have used the energy dependence" 


l =30E!/(1+3.1/E), 


with E in Mev and I> a constant. 

According to the Bloch formula (or to any 
superposition of s- and p-waves alone), the dif- 
ferential cross section as a function of cos@ is a 
parabola. Figures 2-6 are actually plots of do/dw 
vs. cosé, with @=0, cos@=1 at the origin and 
§=x, cos?= —1 at the maximum recoil energy. 
The parabolic fit is good from 0.6 to 1.1 Mev, 
and is not too bad from 1.1 to 1.6 Mev. Further- 
more, one can fit any one of the curves by using 
parameters near those suggested by Staub and 
Tatel from their data: [T=0.4 at 1 Mev, 
Ey2=1.45 Mev, Eyj2=1.05 Mev, oo=1.5 barns, 
for example, yield a good fit at 1 Mev 


(oo =4r x sin?do). 


Sets of parameters not very different from Staub 
and Tatel’s may be found to fit any of the curves. 
Also, fits may be obtained both with the assump- 


"We are grateful to Professor V. F. Weisskopf for his 
unpublished derivation of this relation, based on the 
analyses of V. F. Weisskopf and D. H. Ewing, Phys. Rev. 
57, 472 (1940). 
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tion of an inverted doublet (E3;2< E,/2) and with 
the assumption of a normal doublet (£3;2> E£y2). 
But as yet no set of parameters has been found 
which, with the energy dependence of [ and é>» 
mentioned, has provided a reasonably good fit 
for all the curves from 0.8 to 1.6 Mev; those 
parameters which suit the lower energy curves 
fit badly at the higher, and vice versa. As long 
as a good fit is not found, nothing can be said 
from this data as to whether the doublet is 
normal or inverted, or even as to the funda- 
mental nature of the scattering. In particular, 
in addition to resonant p-scattering which should 
have the form of the Bloch formula, there still 
may be a significant contribution from a fairly 
constant p-background of non-resonant nature, 
a possibility which is not studied in this paper. 

The data of Staub and Tatel indicate that the 
p-splitting may be significant, giving a double 
peak. (It is possible to draw a smooth single-peak 
curve through their data going scarcely beyond 
experimental error, but the points themselves 
indicate a double peak.) Figure 9 also presents 
the appearance of a double peak, although the 
accuracy is not good enough to assure that it is 
real. But if one assumes the correctness of the 
Bloch formula, Figs. 2-6 furnish indirect evidence 
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Fic. 9. Helium total scattering cross section. 
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that the peak is double. In the case of zero 
splitting and neutron energy equal to the 
resonant energy, the Bloch formula reduces to 


da /dw = X*(sin?59+9 cos?@+6 cosé@ sin?do). 


Then if we take the ratio of differential cross 
section at #= 7 and at the minimum of the curve, 


we get 
Or 9—5 sin7d, 








The value of sin?é9 is probably about 0.3. (We 
have o,=4k? sin’59; oo has been measured to 
be” 1.25X10-*4 cm? at thermal energies and 
probably remains the same or decreases slowly 
as energy increases; A?~1/E.) But in any event 
the ratio (9—5 sin*5o) /(sin?é9—sin‘d9) cannot be 
less than 25, the minimum value reached at 
sin?d9=0.6. The largest ratio ¢,/¢min observed in 
the curves is only about 3, and the resonance 
cannot be so sharp that all the curves miss it 
badly. Therefore, on the assumption of the Bloch 
theory, it is necessary to have a splitting of the 
p-level in order to bring the theoretical ratio 
down. 

There are two more points of interest in the 
data. The kinks in Figs. 4 and 5 were repeated 
too well to be caused by statistical fluctuation. 
They could be caused by d-waves except for the 
fact that the energy is still too low for strong 
d-scattering. There is no possibility of explaining 
them by simple s- and p-scattering. It is possible 
that such deviations from a smooth curve could 
be caused by the presence of a group of neutrons 
of lower energy, but other experiments have 
given no evidence for the existence of such a 
group. 

A second point of interest is indicated in Fig. 8. 


2 J. Schwinger, Phys. Rev. 58, 1004 (1940). 
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The total cross sections are those of Fig. 9, 
averaged. The peak heights, taken from Fig, 7 
but plotted on an arbitrary scale, indicate the 
variation in the backward scattering crogs 
sections. The Bloch formula predicts that 
forward and backward scattering cross sections 
will rise faster than the total cross section as 
resonances are approached. (This will be so 
unless the p—s interference term is too large, in 
which case either the forward or the backward 
cross section will rise slower than the total, 
while the other will rise still faster. It is unlikely 
that the interference term is this large.) In con- 
firmation of the prediction, the peak heights of 
Fig. 8 rise faster than the integral counts. 


CONCLUSION 


The interpretation given above must be con- 
sidered as preliminary comment on the theory; 
it remains to be seen whether or not a more 
systematic approach will yield a detailed check 
between theory and the experimental distribu- 
tion curves. At present, the data have not deter- 
mined the adequacy of s- and p-scattering with 
spin-orbit interaction in treating the 1-Mey 
helium resonances, but it does indicate, assuming 
this adequacy, that there is a splitting of the 
p-level. 
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Measurements on the distribution of neutrons in the atmosphere have been made to an 
altitude of 36,600 feet. Observations were made with counters filled with boron trifluoride 
containing 82 percent of B'*. Information concerning the spectral distribution of the neutrons 
was obtained by using the counters in various shielding geometries. Counters were used bare 
(no shielding), cadmium-covered, boron-covered, and surrounded with a layer of paraffin. 
The counting rates of all counters varied in the same way with altitude in the pressure range 
from 20 to 50 cm Hg. Within this pressure range the intensity J is given quite accurately by 
the relation J] =Ige~°-°”, where p is expressed in cm of Hg. The counting equipment was 


carried in a B-29 type aircraft. 


1. INTRODUCTION 


EVERAL previous experiments! have been 

performed to study the neutrons associated 
with the cosmic radiation, and some theoretical 
discussions® of their origin have been published. 
The purpose of the experiments described here 
was to observe the dependence upon altitude of 
the total neutron flux and of the general spectral 
distribution of these neutrons, in the range of 
altitudes easily accessible to the B-29 airplane 
in which the equipment was carried. Among 
other quantities, the “cadmium ratio” was 
measured at various altitudes: The cadmium 
ratio is the ratio of the neutron counting rate of 
a bare counter to that of the same counter, 
cadmium-shielded, under otherwise identical 


* Now at the University of Chicago. 
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Phys. Rev. 56, 10 (1939); S. A. Korff, Phys. Rev. 56, 210 
(1939); Rev. Mod. Phys. 11, 211 (1939); Proc. Am. Phil. 
Soc. 84, 589 (1941); S. A. Korff and W. E. Danforth, Phys. 
Rev. 55, 980 (1939); H. v. Halban, M. Magat, and 
L. Kowarski, Comptes rendus 208, 572 (1939); S. A. Korff 
and E. T. Clarke, Phys. Rev. 61, 422 (1942); S. A. Korff 
and M. Kupferberg, Phys. Rev. 65, 253 (A) (1944); M. M. 
Shapiro, Rev. Mod. Phys. 13, 58 (1941); M. Blau and 
H. Wambacher, Akad. Wiss. Wien 146, 469 (1937); A. 
Filippov, A. Jdanov, and I. Gurevich, J. Phys. U.S.S.R. 1, 
51 (1939) ; W. Heitler, C. F. Powell, and H. Heitler, Nature 
146, 65 (1940); S. A. Korff and B. Hamermesh, Phys. Rev. 
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*H. A. Bethe, S. A. Korff, and G. Placzek, Phys. Rev. 
57, 573 (1940); ibid., 69, 158 (1946). E, Bragge, Ann. d. 
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conditions. Thus, this ratio is a measure of the 
fraction of neutrons in the thermal region of 
energy. ‘ 


2. APPARATUS AND EXPERIMENTAL METHOD 


The neutron detectors were proportional 
counters, 1} inches in diameter and 11 inches in 
effective length, constructed with thin brass walls 
and a 0.0025-inch center wire. The counters were 
filled with BF; to a pressure of 20 cm of Hg. 
Some counters were filled with BF; enriched to 
82 percent in the B” isotope, while others were 
filled with gas containing boron of normal iso- 
topic constitution. 

All counters gave counting rates essentially 
independent of voltage when operated in the 
region of 1550 to 1650 volts and, since they con- 
tained equal quantities of boron, counting rates 
in a given flux of slow neutrons were in the ratio 
of the B" contents. Those counters containing 
82 percent of B” had an efficiency of about 2 
percent for thermal neutrons. Within the sta- 
tistical counting errors in the experiments 
described here, all of the counters containing 
enriched boron showed equal efficiencies for the 
detection of slow neutrons in the laboratory. 
The counters containing ordinary BF; were used 
in the flights only to check the assumption that 
all counts above background were produced by 
neutrons and not by other cosmic-ray phenomena 
such as showers. It was assumed that the effects 
of processes other than neutron capture were 
negligible if counters in identical geometries gave 
counting rates proportional to their B' contents 
at all altitudes. 




































Four counters were used in different geometries 
during each flight in order that a rough deter- 
mination of the energy distribution of the neu- 
trons might be made. One counter was used 
without shielding. Its response was expected to 
vary in the same way with neutron energy as the 
B” capture cross section, i.e., the sensitivity was 
expected to be approximately inversely propor- 
tional to the velocity of the neutrons. Another 
counter was covered with 0.030-inch cadmium 
sheet. This shield absorbed neutrons with ener- 
gies less than about 0.4 ev but had little effect 
on the response of the counter to neutrons of 
greater energy. A third counter was surrounded 
by a layer one inch thick of B,C (normal boron) 
packed to a specific gravity of 1.3. Because this 
shield exhibits an absorption inversely propor- 
tional to neutron velocity, the enclosed counter 
was effectively protected against neutrons of 
energies less than several tens of electron volts. 
The fourth counter was enclosed in a cadmium- 
covered cylindrical shield of paraffin with 2-inch 
wall thickness. This counter was primarily a fast 
neutron detector since its sensitivity was roughly 
independent of neutron energy in the range from 
a few electron volts to about a million electron 
volts. 








204 AGNEW, BRIGHT, AND FROMAN 









Precautions were taken to insure that the 
counters were influenced as little as possible by 
their surroundings in the airplane. The bare, 
cadmium-shielded and boron-shielded counters 
were placed in the pressurized tail gunner’s com. 
partment of the B-29. All gunnery equipment 
and hydrogenous material were removed from 
this section, and the armor glass of the observa- 
tion ports was replaced by Duralumin and steel, 
The effects of the counters on each other were 
minimized by spacing them about 2 feet apart, 
A small high voltage battery and preamplifiers 
for the counter signals were the only other pieces 
of apparatus in the tail compartment. Removal 
of the kapok insulating layer from the tail 
resulted in rather low temperatures in this section 
of the plane during flight. In order to insure that 
the counter voltage remained constant, the 
battery was enclosed in a steel box lined with 
hair felt and heated electrically from the plane's 
28-volt line. The counters themselves were 
checked in the laboratory and found to be insen- 
sitive to temperature in the range from —55°C 
to 25°C. 

The paraffin-covered counter could not be 
placed near any of the others during the experi- 
ment because of the effect of the paraffin on 


TABLE I. Neutron counting rates observed during flights. 











Enriched to Cadmium 
Date of Atmospheric Altitude Pressure Counts per minute corrected for background normal ratio ratio 
flight condition (ft.) (cm Hg) A(E) Cd(E) B«&C(E) A(N) BsC(N) P(E) A(E)/A(N) A(E)/Cd(B) 
5/8/46 Clear 36,650 16.6 25.37 = _- — 332 -- 
34,700 18.1 21.57 _ 5.06 _ 301 4.30 
29,700 22.9 15.37 _ 3.86 216 3.98 
*24,.850 28.4 9.67 —_ 2.56 _ 150 3.78 
20,100 34.8 5.45 _ —0.95 _ 82 _ 
BiC(E)/BuC(N) 
5/21/46 Clear 33,630 19.1 18.47 9.14 6.45 1.59 322 4.06 2.02 
*29,770 22.9 15.53 6.74 4.35 1.23 217.3 3.54 2.31 
24,830 28.3 11.25 4.84 2.39 0.46 159 5.18 2.32 
*19.940 35.0 6.32 2.84 1.59 0.39 87 4.07 2.22 
*15,000 43.0 2.97 1.13 0.72 0.32 _ 2.25 2.63 
14,840 43.2 2.54 1.02 —0.005 0.01 42.6 _— 2.49 
*9,980 52.5 0.57 0.7 0.05 —0.37 16 _ 0.81 
4/11/46 Clear 29,900 22.8 17.67 7.93 250 2.22 
*26,900 26.0 12.47 5.96 184 2.09 
24,780 28.5 9.97 4.79 154 2.08 
*22,000 32.0 7.Al 3.65 113 2.03 
19,800 35.2 6.25 2.45 85.8 2.55 
*17,000 39.7 4.01 1.93 59.2 2.08 
*12,000 48.4 1.34 0.77 26.2 1.74 
3/20/46 Clouds at 22,000 29,800 22.8 16,97 8.84 3.74 252 1.92 
24,930 28.3 12.67 6.21 2.74 164 2.04 
14,930 43.2 4.74 2.05 1.02 45 2.31 
4/5/46 Cloud up to 25,000 ft., 10,000 10.07 6.91 7.08 20.6 1.45 
flying 200 ft. under 10,300 3.37 1.97 1.92 25.6 1.71 
toward hole 10,300 4.57 1.77 1.92 24.0 2.58 
10,380 2.97 0.97 1.52 24.0 3.06 
10,420 2.27 1.17 1.52 19.2 1.90 








Note:Counter filling indicated by 
(E) =enriched boron 
(N) =normal boron 


*Coming down. 


Counter shielding indicated by 
A =bare (no shield) 
Cd =cadmium shield 


B.C =boron carbide shield 
P =paraffin shield 
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Fic. 1. Neutron distribution in the atmosphere. 


their counting rates. Thus, this counter was 
suspended from the ceiling of the radar room near 
hydrogeneous material whose geometry did not 
remain constant during a run. However, the 
0.030-inch cadmium cover surrounding the 
paraffin shield rendered this counter insensitive 
to neutrons thermalized in the vicinity. Initially 
it was feared that the gasoline carried by the 
plane (some 6500 gallons) would affect the 
counting rates. However, measurements taken at 
the same altitude at the start and at the end of 
a flight, after the gas load was some 3,500 gallons 
less, gave the same results even for the paraffin- 
shielded counter in the radar room. The radar 
room is much nearer the gasoline tanks than is 
the tail compartment. Consequently, it is be- 
lieved that the counters in the tail compartment 
can be considered to be in free space, except for 
the dural tail structure. 

The recording equipment consisted of pream- 
plifiers placed near their respective counters and 
linear amplifiers, scalers, and mechanical regis- 
ters. The filament supply lines for the pream- 





plifiers and the co-axial signal leads were strung 
from the tail compartment to the recording 
equipment which was placed in the radar room. 
Only four sets of recording equipment were used, 
so that it was not possible to use both the 
enriched and normal counters with all types of 
shielding simultaneously. The normal boron 
counters were used only in replacement for the 
enriched bare counter and the enriched boron- 
shielded counter. 


3. RESULTS AND CONCLUSIONS 


The ratio of the slow neutron efficiencies of 
the enriched and normal boron counters was 
measured in the laboratory. This ratio was com- 
pared with that obtained at various altitudes for 
both the bare and the boron-shielded counters 
and was found to remain constant within the 
statistical counting errors. This meant that the 
counters were counting neutrons plus natural 
background at all times, and that the effects of 
showers were negligible. It is clear that, if one- 
inch thickness of B,C did not produce enough 
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showers to be detectable, a thickness of 0.030 
inch of cadmium would not either. Hence, the 
cadmium-covered enriched boron counter was 
never replaced by a normal boron counter. 

The observed counting rates obtained in five 
flights are given in Table I. In two of these 
flights the recording equipment for one of the 
counters failed, and in one flight counting rates 
were obtained with five counters by using the 
same recording equipment part time at each 
altitude for two different counters. The geom- 
etries used for each counter are indicated in the 
table. Each of the counting rates given in Table I 
is obtained from a run of about 30 minutes. The 
statistical error can be estimated in any par- 
ticular case from the total number of counts 
obtained by multiplying the counting rates given 
in the table by a factor of 30. Figure 1 shows 
curves of counting rate plotted against atmos- 
pheric pressure for the four enriched counters. 
These curves are parallel within experimental 
error. The curve for the boron-shielded counter 
includes data from only two of the flights as 
indicated in Table I. 

The two curves with titles to the left in Fig. 1 
are included for comparison and are plotted to 
an arbitrary scale. The upper curve for ionization 
caused by cosmic rays? is parallel to the curves 
obtained in these experiments. The lower of these 
two curves represents recent data obtained by 
Korff and Hamermesh! by means of an un- 
shielded BF; counter carried by an unmanned 
balloon. The difference in slope between this 
curve and the others shown may arise from the 
relatively poor statistical accuracy, especially at 
the lower altitudes, which is attained in balloon 
flights. 

One flight was made during cloudy weather, 
and the condition of the aircraft kept it below 


*R. A. Millikan, H. V. Neher, and S. K. Haynes, Phys. 
Rev. 50, 992 (1936). 
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the clouds. It was extremely difficult to note 
the plane’s position relative to the main body of 
clouds, as it was about 200 feet beneath them at 
all times. The only data obtained were taken 
under an extremely thick cloud. Observations 
were continued while the plane flew toward ap 
opening in the overcast. Directly under the body 
of the cloud, the counting rate was high, except 
in the paraffin-covered monitor. As the opening 
was approached, the counting rate in the other 
counters decreased. It is hoped that this effect 
can be checked in later experiments incorporating 
a humidity measuring device. Data from this 
flight are not plotted in Fig. 1. 

Within the error of these experiments the 
cadmium ratio is found to be constant, with a 
value about 2.2, at all altitudes from 7000 to 
34,000 feet above sea level. The value of the 
ratio at the lower of these altitudes was obtained 
in the laboratory and is not shown in Table |, 
It is clear from the slopes of the graphs shown 
in Fig. 1 that the neutron counting rate of all 
the counters varied in the same way with the 
altitude as does the total ionization produced by 
the cosmic radiation. Thus, no variation of 
spectral distribution with altitude was observed, 
although it is evident that the experiment was 
not designed to detect variations in the neutron 
energy distribution in the region above 100 ev 
or so. In the pressure range 20 to 50 cm of Hg 
the neutron flux, 7, observed with all counters 
is quite accurately given by J = J e~°-°*?, where 
p is the pressure expressed in cm of Hg. 

The authors wish to express their appreciation 
to the officers and men of the 428th A.AF. 
Base Unit (Special) who were so efficient and 
cooperative in flying a B-29 for this work. 

This paper is based on work performed under 
Contract No. W-7405-eng-36 with the Man- 
hattan District at the Los Alamos Scientific 
Laboratory of the University of California. 
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Experiments are described by means of which the collision cross sections of carbon and 


hydrogen were determined for neutrons of nine different energies between 6 Mev and 22 Mev. 
The experimental method was such that data could be taken on groups of neutrons having 
energies within an interval small compared to the neutron energy, and at any energy produced 
by an energy heterogeneous source. The results are consistent with those of other observers 
at the four neutron energies in this range which have been previously investigated. A com- 
parison with various theories shows that at high energies the observed hydrogen cross sections 
are about 10 percent larger than pure S-wave cross sections computed on the basis of a square 
well interaction between neutron and proton, and that they agree very well with the “sym- 


metrical” theory of Rarita and Schwinger. 








1. INTRODUCTION 


EUTRON sources in the 5- to 25-Mev 

energy range available up to 1943, when 
the present experiments were carried out, 
produced neutrons distributed continuously in 
energy. Therefore, in order to measure cross 
sections for neutrons in energy intervals small 
compared to the neutron energy, it was necessary 
to use a detecting system which could dis- 
criminate against neutrons outside of a particular 
interval. The first method which accomplished 
this made use of a nuclear reaction with a 
definite energy threshold as neutron detector. In 
this case, the most energetic neutrons produced 
by the source must be a reasonable interval 
above the detector threshold, and the method is 
limited to those energies at which there are 
suitable threshold reactions. It has been used 
at 7 Mev by Grahame and Seaborg' with 
Fe**(n, p)Mn**, at 12 Mev by Salant and 
Ramsey? with Cu™(n, 2n)Cu®, and at 21 Mev by 
Sherr? with C”(n, 2n)C®. 

Aside from the present results, the only other 
total cross-section data known to the writer on 
neutrons above 6 Mev are those of Ageno, 
Amaldi, Bocciarelli, and Trabacchi.* They meas- 
ured the cross sections of hydrogen, deuterium, 

* Now at the University of Minnesota. 
ass) C. Grahame and G. T. Seaborg, Phys. Rev. 53, 795 
see O. Salant and N. F. Ramsey, Phys. Rev. 57, 1075A 

*R. Sherr, Phys. Rev. 68, 240 (1945). 

*M. Ageno, E. Amaldi, D. Bocciarelli, and G. C. Tra- 


bacchi, Naturwiss. 31, 231 (1943); Phys. Rev. 71, 20 
(1947). 
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and carbon at 4.1, 12.5, and 13.5 Mev with 
neutrons from Be, B, and Li targets, respectively, 
bombarded with 1-Mev deuterons. Their de- 
tecting method was similar in principle to that 
of the present experiments. However, since only 
1-Mev deuterons were available, they were re- 
stricted to energies near those released in the 
nuclear reactions. 

A very complete study of the cross sections of 
carbon, hydrogen, deuterium, and oxygen be- 
tween 0.35 Mev and 6.0 Mev has been made by 
Bailey et al. Above 2 Mev, these workers also 
used a heterogeneous neutron source and a dif- 
ferential method, somewhat different from that 
described below, to observe the effects of the 
neutrons in small energy intervals. 

The results obtained by these observers on 
carbon and hydrogen for neutron energies 
greater than 5 Mev are shown in Figs. 3 and 4, 
along with the results of the present experiments. 
Theoretical work relating data of this sort to 
the problem of nuclear forces will be referred to, 
and the experimental results compared with 
those of various theories, in the last section of 
the paper. 


2. EXPERIMENTAL METHOD AND APPARATUS 


The fast neutron collision cross section of a 
nucleus can be determined most easily by means 
of a transmission experiment, that is, one in 
which the fraction of a neutron beam is measured 


’ Carl L. Bailey, W. E. Bennett, Thor Bergstrahl, R. G. 
Nuckolls, H. T. Richards, and J. H. Williams, Phys. Rev. 
70, 583 (1946); 70, 805 (1946). 
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which a known amount of the element trans- 
mits unaltered in direction and in energy. In 
general, cross sections obtained in this way will 
be the sum of cross sections for three different 
processes : elastic scattering, inelastic scattering, 
and radiative capture. However, the capture 
cross sections of both hydrogen and carbon for 
neutrons of the energies dealt with in the present 
experiment are negligible.* In the case of hy- 
drogen, moreover, inelastic scattering does not 
take place, so that a transmission measurement 
will give just the total scattering cross section 
appropriate for comparison with theory. The 
result of this kind of experiment will not make 
possible the separate evaluation of the two cross 
sections for carbon unless some other information 
can also be used. (This question is discussed 
below in Section 4.) Transmission measurements 
are particularly advantageous when the source 
is such that the detector must define the energy 
interval, since in a transmission experiment the 
detector need not be altered in any way between 
the beam strength measurements to be com- 
pared, and its energy response will necessarily be 
the same. Thus the method is free of the source 
of error which has been most serious for angular 
distribution measurements.** 

Neutrons for the present experiments were ob- 
tained by bombarding beryllium and lithium 
with 10-Mev deuterons produced by the Michi- 
gan cyclotron. In the first case the reaction is 
Be*+D?—B"+n!, and the most energetic neu- 
trons had about 14-Mev energy. Neutrons from 
this reaction were used to study three energy 
intervals between 6 and 11 Mev. They were not 
used at higher energies because of the relatively 
small number near the upper energy limit. The 


*H. A. Bethe, Rev. Mod. Phys. 9, 160 (1937). 

** For example, in the latter experiments on hydrogen 
(references 7 and 8), the recoil proton intensity is measured 
at two different angles with respect to the neutron beam, 
and the angular distribution is then inferred from a com- 
parison of these two intensities. If the comparison is to be 
significant, the neutron energy intervals used at the two 
angles must be exactly the same because the neutron 
intensity in the interval depends very strongly on where 
the interval is cut off at the low energy end. This means 
that the proton cut-off energy will have to be different 
at the two angles by exactly the right amount. This adjust- 
ment can cause serious errors because, for example, it 
depends quantitatively on the somewhat uncertain range- 
energy cb a of protons in aluminum. 

7 Howard Tatel, Phys. Rev. 61, 450 (1942). 

8E. Amaldi, D. Bocciarelli, B. Ferretti, and G. C. 
Trabacchi, Naturwiss. 30, 582 (1942). 


lithium reaction, Li7+D*—-Be‘+n!', produced 
neutrons having energies up to about 25 Mey 
with which 7 energy intervals between 10 and 
23 Mev were investigated. 

The detecting method used was one which wil] 
determine the number of neutrons in an energy 
interval of almost arbitrary size at any energy, 
The essence of the method is this: the protons 
in any energy interval, which are knocked out of 
a thin paraffin layer at a particular angle, are 
all due to neutrons in a corresponding but slightly 
larger energy interval. (The neutron energy 
interval is larger by a small factor depending on 
the scattering angle and also by an additive 
amount equal to the energy which a proton can 
lose in going through the paraffin layer.) The 
intensity of the neutron beam is measured, in 
units of unknown but constant absolute mag. 
nitude, by counting these recoil protons. To 
determine the neutron intensity in a particular 
energy interval it is necessary to take the dif- 
ference between two sets of counts during which 
the same number of neutrons, distributed in the 
same way with respect to energy, are produced 
by the source. This can be approximated by 
keeping conditions at the source constant, and 
then measuring while two equal charges of 
deuterons fall on the cyclotron target. Between 
the paraffin layer and the ionization chamber are 
placed aluminum absorbers such that one of 
these sets of counts measures the total number of 
neutrons having energies above the lower limit 
of the interval, while the other measures the 
number which have energies above the upper 
limit of the interval. 

The efficiency of this system is necessarily low 
because only a very small fraction of the neutrons 
will collide with protons in a layer of paraffin 
thin enough so that there will be a relatively 
small energy difference between protons pro- 
duced on the two sides of the layer by neutrons 
of the same energy. Also, unless the angles sub- 
tended by the detector at the paraffin layer and 
by the paraffin layer at the detector are small, a 
large energy spread will be introduced by the 
collisions which occur at widely different angles. 
Thus the number of protons observed is small, 
and it is important to have as large a ratio as 
possible between the pairs of counts whose dif- 
ference is to be taken. When the Li target was 
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used, the largest values of this ratio occurred 
when the upper limit of the interval was some- 
where near the energy of the most energetic 
neutrons produced by the source. In general, 
therefore, enough Al was put in front of the Li 
target in the cyclotron to bring about an ap- 
proximate coincidence of these two energies. 
Under these conditions, a thick paraffin layer 
would do as well as a thin one, since only those 
protons losing very little energy in the paraffin 
would be able to get through the Al absorbers. 

However, with neutrons from the Be target 
the best values of this ratio occurred when the 
upper limit of the energy interval was several 
Mev below the peak energy of the source. Also, 
in order to work on the 9-13-Mev region with 
the more abundant Li neutrons, it was necessary 
to cut off the energy interval below the peak 
energy of the source, which could not be made 
as low as the 14 Mev released in the reaction. 
Therefore a thin paraffin layer (one in which the 
protons lost on the average about 4 Mev) was 
used, making it possible for the detector to put 
an upper as well as a lower limit on the working 
energy interval. 

The arrangement of the apparatus is shown 
schematically in Fig. 1 for the low energy runs 
in which beryllium was used as the cyclotron 
target. The set-up used at high energies with the 
lithium target differed from this only in the 
respects noted below. The following lettered 
remarks correspond to the letters in the figure. 
(a) The arrow represents the deuteron beam 
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NOTE: DIMENSIONS INCM 


NOT TO SCALE. 


produced by the cyclotron, with energy about 10 
Mev and intensity 8-15 microamperes. 

(b) This line represents the aluminum ab- 
sorbers of various thicknesses which were put in 
front of the target during most of the runs with 
Li to slow down the deuterons. 

(c) The cyclotron target was a piece of beryl- 
lium or lithium about 3.0 cm. wide, 2.0 cm high, 
and 0.4 cm thick, screwed tightly to the water- 
cooled brass target cup. The dimensions of the 
effective neutron source, i.e., the part of the 
target receiving most of the deuteron beam, were 
about 2.0 cm horizontally and 0.4 cm vertically. 

(d) The scatterers were blocks of graphite or 
paraffin about 15 cm square, set up 282 cm from 
the cyclotron target and 154 cm from the paraffin 
layer. Various thicknesses were used. The deter- 
mination of their composition and density will 
be described in the next section. 

(e) The paraffin layer from which protons were 
knocked by a small fraction of the neutrons was, 
for the runs with the Be target 9 cm high, 3 cm 
wide (perpendicular to the recoil protons), and 
0.01 cm thick. When Li was used as the cyclotron 
target, the layer was 12 cm high, 2 cm wide, and 
0.028 cm thick. 

(f) At low energies an evacuated box, labeled 
“collimator,”’ with an aluminum window to let 
the protons out, contained the paraffin layer and 
extended nearly to the ionization chamber. Its 
angle and length were such as to put the ioniza- 
tion chamber outside of the direct neutron beam 
(in order to cut down the number of recoiling 
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nuclei in the chamber), and to cause the sensitive 
region of the chambers to subtend a sufficiently 
small angle at the paraffin layer so as not to 
widen the energy interval intolerably. At the 
higher energies this box was omitted, and an air 
space of 32.7 cm took the place of some of the 
aluminum absorber. In this case the center of the 
proton beam made an angle of 31.4° with that of 
the neutron beam. The minimum proton energy 
necessary to cause a coincidence is computed for 
the two cases in Table I of Section 3. 

(g) The aluminum absorbers defining the 
limits of the energy intervals were put into the 
proton path at this line. The thicknesses used, 
and the resulting proton and neutron energies, 
are shown in Table II of Section 3. 

(h) The double ionization chamber and circuits 
arranged to record only coincidences were neces- 
sary because in a single chamber the background 
count in the cyclotron room was much too high. 
The chamber and circuits will now be described. 

The coincidence ionization chamber was built 
and used by H. Tatel to measure the angular 
distribution of protons scattered by 11-Mev 
neutrons. Its construction has been described 
in Tatel’s article,’ so that for our purposes 
a statement of the critical dimensions and the 
conditions under which it was used in the present 
experiments will be sufficient. The general shape 
and most of the dimensions can be seen in Fig. 3. 

Each of the two sensitive regions is 4 cm high, 
2 cm wide in the direction of the electric field, 
and 2.3 cm deep along the proton path. Cor- 
responding points in the two chambers are 4.2 
cm apart along the proton path. Between them, 
perpendicular to the proton path, is a thin 
aluminum shield. The protons entered the 
chamber through an aluminum window (1 mil 
thick for the low energy runs and 2 mils thick 
for those at high energy), and then had to travel 
4.5 cm before reaching the sensitive region of the 
first chamber. Throughout the present experi- 
ments the chambers were filled with oxygen at 
1.6-atmospheres pressure, and the ion collecting 
potential was about 5700 volts. A test of the 
dependence of the counting rate on the collecting 
voltage showed that with the electrical sensi- 
tivity at its standard value, and an oxygen 
pressure of 1.6 atmos., the voltage had to drop 
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below 5100 before the counting rate decreased 
observably. } 

The low voltage collecting plate of each 
chamber was connected to the first grid of , 
four stage resistance-capacitance coupled pulg 
amplifier, the gain of which could be varied. The 
amplifiers were connected through discriminator 
to a Rossi coincidence circuit, the output of 
which went to a scaling circuit and mechanicaj 
counter. The coincidence resolving time of this 
system was measured by Tatel with alpha. 
particle sources of various strengths and found 
to be 1.2 X10~ seconds. 

Since the method used to determine trang. 
mission fractions required that all protons which 
passed through the chambers be counted, it was 
necessary to measure the over-all sensitivity of 
the counting systems to discover what proton 
energies were allowable. Accurately reproducible 
settings of the electrical sensitivity, which 
depended on the relation between amplifier gain 
and discriminator bias, were made by means of 
an 1000-cycle oscillator in the manner described 
by Tatel. The operating value was taken to be 
the most sensitive setting at which no noise 
pulses were counted (this test was made with the 
cyclotron running). Then the over-all sensitivity 
of each counting system, i.e., the minimum 
energy loss a particle must have in the sensitive 
region of the chamber in order to produce a 
countable pulse, was determined by tests made 
at reduced pressure with polonium alpha-par- 
ticles. The first chamber and its circuits were 
found to be such that under operating conditions 
(electrical sensitivity as above, 1.6 atmos. Oy, 
and 5700 volts) a particle would be counted if it 
lost more than 78 kev/cm in the sensitive region. 
This is equivalent to 44 kev/cm in standard air, 
so that the most energetic protons which the 
first chamber could count were those which 
had about 11.6 Mev as they entered the sensitive 
region. (For alpha-particles the range-energy 
relation given by Holloway and Livingston® was 
used, and for protons that of Livingston and 
Bethe.) To be detected by the second chamber 
and its circuits, a particle had to lose 73 kev/cm 


*M. G. Holloway and M. S. Livingston, Phys. Rev. 4, 


31 (1938). 
1 M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 


9, 268 (1937). 
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or 168 kev in the 2.3 cm long sensitive region. 
This is the energy loss of a proton in the last 
0.24 air cm of its range, so that the least energetic 
protons which could cause coincidences were 
those whose paths ended about 0.24 air cm 
inside the sensitive region of the second chamber. 
These protons had about 2.5 Mev as they entered 
the sensitive region of the first chamber. 

Thus the system as it was used would count 
any protons which, as they entered the sensitive 
region of chamber 1, had energies between 2.5 
Mev and 11.6 Mev. Because of the curvature of 
the range-energy relation, this will correspond to 
a smaller energy difference in the protons before 
they have passed through the Al absorbers. 
However, the initial energy range of the protons 
which the chambers could count was ample even 
for those cases in which a working interval was 
used, the upper limit of which was several Mev 
below the energy of the fastest neutrons produced 


by the source. 
3. MEASUREMENTS 


The composition of the paraffin scatterers was 
determined by means of melting point measure- 
ments. Two samples of the paraffin changed 
from completely liquid to completely solid 
between 60° and 56°C. This implies an average 
chain length of 26.6 carbons, and gives for the 
ratio, f, of the number of H atoms to the number 
of C atoms, f = 2.075+0.2 percent. The scatterers 
used were two blocks of this paraffin about 6 
inches square and 13 inches thick. The surfaces 
through which the neutrons were to pass were 
scraped flat and parallel (to +0.001 inch) with a 
steel straight edge. Measurements of size and 
mass gave the mass per unit area to within 0.1 
percent. For one block this was 3.163 g/cm?, and 
for the other, 3.139 g/cm?. The graphite used for 
the carbon scatterers was obtained from the 
National Carbon Company of Cleveland, and 
according to their analysis contained less than 
0.1 percent of substances other than carbon. It 
was in the form of pieces } inch thick and 6 inches 
square whose surfaces had been machined flat 
and parallel, and the mass per unit area of each 
was determined as above. 

The calculation and final values of the neutron 
energy intervals used are shown in Tables I and 
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Il. These may be made clearer by the following 
remarks. 

The difference between two counts will give 
the number of recoil protons in an energy 
interval (E,o—E,:) which has been sharply 
cut off at both ends by the aluminum absorbers. 
The lower and upper limits, Eo and E,, of the 
neutron energy interval to which this corresponds 
are given by the expressions 


Eo => Ey0/cos?6o, Ey, = (E51 +AE,) /cos*@2. 


Here 69 is the minimum angle, and 6, the maxi- 
mum angle, between the paths of neutron and 
proton at which the proton can reach the ioniza- 
tion chamber, and AE, is the energy a proton 
loses in going all the way through the paraffin 
layer. The probability that a neutron of a par- 
ticular energy will produce a proton in the 
interval E,o—Ey:, i.e., the efficiency of the 
neutron detector, is not the same for all neutron 
energies between Ey» and Ey, but behaves some- 
what like the solid curve in Fig. 2. The dotted 
curve shows how it would look if #=6) and 
AE,=0. Here E, = E,:/cos*6. 

The actual effective energy distribution of the 
neutrons in an interval is the product of this 
detector efficiency and the distribution function 
of the neutrons produced by the source. Since 
this function was always greater at the low 
energy end of the interval, the effective neutron 
energy is lower than the center of the efficiency 
curve. The available information about the 
neutron distribution was sufficient to justify a 
guess at the effective neutron energy for each 
interval (last column of Table II) but not to 
locate it exactly. 

At each new energy interval, various Al thick- 
nesses were tried in the target cup, and some- 
times the absorbers in front of the chamber were 
also changed, until what appeared to be the 
optimum combination was discovered. The 
criterion used in making this adjustment was not 
the relative position of the upper limit of the 
interval and the peak energy of the source 
(though in the runs on lithium these probably 
coincided fairly closely) but the ratio of the 
number of neutrons (per deuteron) in the interval 
to the number above it. 

When this adjustment was such as to produce 
an adequate ratio for an appropriate energy 
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interval, the fraction of the neutrons trans- 
mitted in this interval by a particular scatterer 
could be determined. In the discussion of this 
process the word “count’’ will mean a deter- 
mination of the number of coincidences which 
occur while a definite number of deuterons, 
measured in arbitrary but constant integrator 
units, fall on the cyclotron target, for a particular 
configuration of scatterer and absorbers (e.g., 
scatterer in the beam and upper limit absorber 
in front of chamber). A “‘run”’ is a determination 
of the transmission fraction of a particular 
scatterer for the neutrons in a particular interval, 
and always consists of several of each of four 
different kinds of counts. The scatterer-absorber 
arrangements for these counts are (a) lower limit 
absorber in front of chamber and scatterer not 
in beam, (b) upper limit absorber in front of 
chamber and scatterer not in beam, (c) and (d), 
similar to (a) and (b), respectively, but with the 
scatterer in the beam. Let A be the total number 
of coincidences observed, and WN, the total 
number of integrator units of deuterons striking 
the cyclotron target during all the counts of 
type (a) in a run, and let B, Ng, C, Ne, D, and 
Np be the corresponding quantities for counts 
of the other three types. Then the observed frac- 
tion 7.4, of the neutrons in the working interval 
transmitted by the scatterer is given by 


C/Nc—D/No 
 A/Ni—B/Neg 





obs 


It would have been desirable, of course, to 
measure the length of a count in terms of the 
number of neutrons in the working energy 
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interval produced by the source, but this would 
have required essentially another complete 
detecting system. Instead, the length of a count 
was measured in terms of the charge of deuterons 
falling on the cyclotron target by means of g 
current integrator attached to the latter. The 
number of “integrator units” of deuterons was 
proportional to the number of neutrons produced 
as long as cyclotron operating conditions were 
kept absolutely constant. The two factors which 
were observed to have the greatest effect on this 
ratio were (1) the strength of the magnetic field, 
which affected the energy of the deuterons, and 
(2) the pressure in the cyclotron tank on which the 
energy distribution of the deuterons depended, 
Numerous tests showed that by careful adjust- 
ment and operation of the cyclotron these things 
could be kept effectively constant during times 
the order of 30 minutes, but seldom for as much 
as an hour. Therefore, it was desirable to take 
at least one of each of the 4 kinds of counts 
about every half hour. Since the number of coin- 
cidences occurring in a count of 5 to 10 minutes 
was in general not nearly large enough to produce 
a sufficiently accurate value of the transmission 
fraction, at least 5, and usually 8 or 10, counts 
of each kind were taken during each run. The 
value of 7... from each run is given in Tables 
III and IV, Section 4. The relative standard 
deviation in 7.»,, 47, given in these tables was 
computed on the assumptions that the standard 
deviation in each of the numbers A and B was 
equal to its square root, that the uncertainty in 
Na and Nz, was negligible, and that Ng = N4 and 
Np=Nce. (This latter condition was nearly 
always and always nearly fulfilled.) The ex- 


TABLE I. Range of countable protons of least energy. 














Be target (as in Fig. 1) Li target 

Aluminum window, collimator: 0.001 inch or 7.35 mgm/cm? 0 mgm/cm* 
Aluminum window, detector: 0.001 inch or 7.35 mgm/cm? 13.6 mgm/cm* 
Aluminum internal shield detector: 0.0007 inch or 4.30 mgm/cm? 4.3 mgm/cm* 
Total aluminum 19.00 mgm/cm? or 12.5 air cm. 11.8 air cm 
Correction for Al 0.1 air cm 0.1 air cm 
Air distance, collimator to detector 3.5 air cm 
Equivalent range, paraffin layer to detector 32.7 cm X 741/760 (distance X mean 

air press. in atmos.) 31.9 air cm 
Equivalent range in Ist chamber: 8.7 cm X 1.6 X 1.1 (distance XO: press. X Oz st. 

power) 15.3 air cm 15.3 air cm 
Equiv. range to register in 2nd chamber: 0.2 air cm 0.2 air cm 
Total: Minimum range necessary to cause coincidence: 31.6 air cm 59.3 air cm 
Corresponding proton energy: 4.80 Mev 6.88 Mev 
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TABLE II. Calculation of neutron energies. Eo= Eyo/cos*@o, Ey = (Ep: +AE,) /cos*és. 











——— 
Al Al Proton energy (Mev) Neutron energy (Mev) 
thickness, thickness, Loss in Min. Max. Effective 
Run Cyclotron tar; chamber Min. Max. paraffin Eo Ex value 
ousabere target (mils) (mils) Epo En SE» (@=23°) (0:=28.5°) E 
0 ea 
7, 10, 13 Be 0, 6 0-3 4.80 5.91 0.40 5.68 816 6.5 
89 Be 0 3-6 5.91 6.96 0.34 699 944 7.8 
11, 12, 14, 15 Be 0 6-12 6.96 8.73 0.29 8.24 11.70 9.3 
(09 =28.7°) (@2=34.1°) 
1. 16 Li 13 0-6 6.88 8.62 0.90 8.94 13.9 10.6 
17 Li 13 6-12 8.62 10.16 0.78 11,2 15.9 12.8 
2. 26 Li 12.5 6-13 8.62 10.40 0.78 11,2 16.3 12.9 
3, 18 Li 11 12-19 10.16 11.74 0.70 13.2 18.1 14.8 
4 19 Li 8.5 18-25 11.53 12.86 0.63 15.0 19.6 16.5 
5 Li 6 24-30 12.68 13.8 0.56 16.5 20.9 18.0 
20, 23 Li 6 24-31 12.68 14.0 0.56 16.5 21.2 18.1 
6, 25, 21 Li 3 30-38 13.8 15.2 0.50 17.9 22.9 19.6 
22, 24 Li 0 36-46 14.9 16.4 0.45 19.4 24.6 21.1 








pression for AT then becomes 
AT=[(C+D)/(C—D)*+(A+B)/(A —B)*]}}. 


4. CORRECTIONS AND FINAL RESULTS 


The procedure outlined above either eliminates 
or takes account of all significant random or 
accidental errors involved in the determination 
of the transmission fraction. Evidence for this is 
the fact that the differences in values of a cross 
section at a particular energy, obtained from 
runs on different days, and with different scat- 
terer thicknesses, are in general not greater than 
the assigned standard deviations of these values 
(see Table III). 

The method is subject to other errors the 
presence of which cannot be discovered by 
repeating measurements if, as was the case, the 
same geometry is used on all the occasions. The 
following four kinds of systematic errors may 
occur. 

(a) Neutrons which leave the source in such 
a direction as not to hit the scatterers may be 
scattered into the detector by nearby objects 
such as the cyclotron magnet or the water wall. 
Though the water wall could be expected to 
absorb most of these neutrons, a test of the 
magnitude of this effect was made by putting 75 
cm of paraffin between the neutron source and 
the detector. (Judging by the measurements on 
paraffin scatterers this reduced the number of 
neutrons in the beam by a factor of more than 
300.) Counts were then taken at several different 
energy intervals, and at none of them was there 


a significant difference between the number of 
coincidences observed with the lower limit ab- 
sorber in place and that with the upper. 

(b) Since the scatterer and the detector are not 
of negligible size compared to their distances 
from the source and from each other, some 
neutrons scattered by the scatterer will reach 
the detector. The magnitude of this effect was 
calculated as described below, and found to be 
such that the observed values of T are the order 
of 1 percent too large. 

(c) The scatterer will also scatter into the de- 
tector a fraction of those neutrons which have 
already been scattered once by some object 
nearer the source. However, the intensity at the 
scatterer of such neutrons was much less than 
that of the direct beam from the source, because 
on only one side of the source-scatterer line was 
there anything sufficiently close and massive to 
scatter neutrons appreciably (a vertical section 
of the magnet yoke). Since this iron certainly 
had no focusing effect, the “‘apparent source”’ in 
it was much weaker than the actual source. 
Since the direct neutrons which the scatterer 
scattered into the detector made only about a 
1 percent change in 7, the effect of this sub- 
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Fig. 2. Schematic diagram of detector efficiency as a func- 
tion of neutron energy. 
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TABLE III. Calculation of total cross sections for carbon. 










JR. 





@(C) =final value of total neutron-carbon cross section (elastic and inelastic scattering). 


Ao(C) =standard deviation in o(C). 


T =transmission fraction corrected for scattering into detector. 











—_—= 
Mean Scatterer 
neutron thickness a(C) a(C) 
Run energy (xd) A4T (10-* Ao(C) (10-% As(C) 
number (Mev) (g/cm?) T obs i (percent) cm?) (percent) cm?) (percent) 
7 6.5 6.562 0.634 0.630 9.7 1.40 21 
10 6.5 10.936 0.432 0.428 9.6 1.54 11 1.51 10 
8 7.8 6.562 0.596 0.592 10.4 1.59 19 
9 7.8 10.936 0.455 0.450 18 1.45 22 1.51 14 
11 9.3 6.562 0.639 0.635 9.2 1.38 21 
12 9.3 10.936 0.450 0.446 10.1 1.47 12 1.45 10 
1 10.6 13.123 0.395 0.390 9.5 1.43 10 
2 12.9 10.936 0.499 0.494 10.9 1.28 15 
26 12.9 13.123 0.484 0.478 9.3 1.12 13 1.18 10 
3 14.8 10.936 0.512 0.507 8.6 1.24 13 
4 16.5 10.936 0.497 0.492 9.9 1.29 14 
5 18.0 8.749 0.614 0.609 9.9 1.13 20 
6 19.6 8.749 0.561 0.556 8.3 1.34 14 
25 19.6 15.311 0.376 0.370 9.8 1.30 10 8 
24 21.1 13.123 0.468 0.462 7.4 7 10 





— 





stantially smaller number of indirect neutrons 
may be safely ignored. 

(d) In those cases in which the upper limit of 
the energy interval was several Mev below the 
energy of the fastest neutrons produced by the 
source, the question arose whether some high 
energy neutrons might be scattered inelastically 
into the detector with an energy loss which would 
bring them into the working energy interval. 
Three major factors control the magnitude of 
this effect. In the first place, the total number 
of neutrons above the working interval was 
always smaller than the number in the interval, 
and most of these were within a few Mev of the 
upper limit of the interval. Secondly, Weisskopf"' 
has shown that the mean energy after impact of 
neutrons in the 10—20-Mev range which suffer 
inelastic collisions is the order of one-third of 
their initial energy, so that only a small fraction 
lose less than 3 Mev. And finally, even if the 
angular distribution of the inelastically scattered 
neutrons had a pronounced maximum in the 
forward direction, the calculations described 
below show that only a small fraction of them 
would actually hit the detector. Thus it is clear 
that the contribution of inelastically scattered 
neutrons to the observed transmission fractions 
may be ignored. 

To calculate the effect on the transmission 


"V. F. Weisskopf, Phys. Rev. 52, 295 (1937). 


fractions of the neutrons scattered into the 
detector by the scatterer, it is necessary to know 
something about the angular distribution of the 
scattered neutrons, at least at small angles. For 
carbon (and larger) nuclei enough is known 
about this to make a calculation worth while. 
Placzek and Bethe” have suggested that the 
elastic scattering of fast neutrons by compound 
nuclei may be regarded as essentially a diffraction 
of the neutron beam by the nuclei, which cor- 
respond to absorbing disks or spheres in a beam 
of light. On this assumption the angular dis- 
tribution of the elastically scattered neutrons is 
given by 


K (0) =R°LJi(kRO) /0}?, 


where K(@) is the scattering cross section per 
unit solid angle at the angle 0, J; is a Bessel 
function of the first kind, k is the wave number 
of the incident neutrons, and R is the nuclear 
radius. The experimental results of Kikuchi, 
Aoki, and Wakatuki,'* who observed the angular 
distribution of 15-Mev neutrons scattered by 
various elements, agree very well with this 
expression (except for their measurements at 
14°, the smallest angle at which they could 
work, which are probably in error). 


2G. Placzek and H. A Bethe, Phys. Rev. 57, 1075A 
(1940). 

13S. Kikuchi, H. Aoki, and T. Wakatuki, Proc. Phys. 
Math. Soc. Japan 21, 410 (1939); 22, 430 (1940). 
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COLLISION CROSS SECTIONS 


TABLE IV. Calculation of total cross sections for protons. 


_@(H) =final value of total neutron-proton scattering cross section. 
Vo(H) =standard deviation in #(H). 


T =transmission fraction corrected for scattering into detector. 








@(C) 





Mean Scatterer from 
neutron thickness Fig. 3 o(H) a(H) - 
Run energy (xd) AT (10-™ (10-™* 4o(H) (10-* 4o(H) 
number (Mev) (g/cm?) Tobs T (percent) cm?) cm*) (percent) cm?) (percent) 

13 6.5 6.302 0.302 0.300 6.3 1.56 1.40 8 
14 9.3 6.302 0.414 0.410 5.2 1.44 0.90 10 

15 9.3 3.139 0.628 0.625 4.7 1.44 0.99 18 0.92 9 
16 10.6 6.302 0.447 0.443 4.8 1.39 0.78 11 
17 12.8 6.302 0.444 0.440 4.9 1.34 0.83 11 
18 14.8 6.302 0.507 0.502 a2 1.29 0.61 15 
19 16.5 6.302 0.494 0.489 5.2 1.28 0.66 15 
20 18.1 6.302 0.516 0.511 5.0 1.26 0.59 16 

23 18.1 3.163 0.746 0.743 4.5 1.26 0.45 33 0.55 14 
21 19.6 6.302 0.539 0.533 5.3 1.25 0.52 18 
22 21.1 6.302 0.576 0.569 5.0 1.25 0.41 21 











The value of R to be used in the expression is 
related to the cross section ¢, for this elastic 
scattering by o,=7R®. Since inelastic collisions 
also occur, ¢, is less than the total collision cross 
section, o, calculable from the observed trans- 
mission fractions of the carbon scatterers. Sherr® 
has compared the results of his experiments at 
24 Mev, in which he also measured total cross 
sections, with those of Grahame and Seaborg! 
whose method was such that they measured 
cross sections for inelastic scattering only. Sherr 
finds strikingly good agreement between the 
results of the two kinds of experiments if he 
assumes o,=o¢;=}$0. It is reasonable from the 
point of view of the optical analogy, also, that 
the cross section for the diffraction process should 
be equal to the actual area of the disk. 

On the basis of these assumptions the trans- 
mitted intensity was calculated (neutrons scat- 
tered once and twice by the scatterer were taken 
into account). The resulting expression for the 
observed transmission fraction is 


A,(l+h)? K(6) 
Taam e¥4 + (Nee 





Ge 





(No.x)* L(8) 
4 foe ) 


2 oe? 


Here o is the total and o, the elastic scattering 
cross section of carbon, Nx the number of atoms 
per cm? of scatterer, A, the scatterer area (per- 
pendicular to the beam), / and h the distances 


from source to scatterer and from scatterer to 
detector, K(@) the average value of K(@) as @ goes 
from 0° to 6° (6° is the largest value @ can have 
for particles which reach the detector, and K(@) 
is nearly constant in this region), and L(@) the 
corresponding average value of the function L(@) 
which gives the angular distribution of neutrons 
which have been scattered exactly twice. 

For the three inch thick carbon scatterer and 
a neutron energy of 21 Mev, the second term in 
the bracket is 1.2 percent of the first term and 
the third is 0.17 percent of the first. Thus even 
at high energies, where the anisotropy of the 
scattering is greatest, this correction is small 
compared to the statistical fluctuations. How- 
ever, since it changes all the cross sections in the 
same sense, it was computed roughly for each 
of the energies and scatterer thicknesses used, 
and applied to 7T.», to get T as shown in Table 
III. The corrections to the paraffin transmission 
fractions indicated in Table IV were made in the 
same way with the additional assumption that 
scattering by the protons is isotropic in the 
center of gravity system. (Though this is very 
likely contrary to fact, there does not seem to 
be information now available which would make 
any other specific assumption more probable.) 

The sixth column of Table III gives AT, the 
relative standard deviation in 7, computed as 
described above from the observed numbers of 
counts. The column headed @(C) gives the total 
collision cross section for carbon at each different 
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energy. The only significant contribution to 
Ao(C), the relative standard deviation of &(C), 
came from AT, since the standard deviation in 
Nx was less than 0.5 percent. In those cases in 
which two runs were made with the same energy 
interval, the columns headed o(C) and Ao(C) 
give the values for each run. ¢(C) is the average 
of the two o(C)s, weighted according to their 
reliability. 

The final values of the neutron-carbon collision 
cross sections, represented by solid dots, are 
plotted against energy in Fig. 3. The open circles 
show the results of other observers. (All available 
data above 5 Mev are plotted.) The total length 
of the vertical line is twice the standard deviation 
(resulting from statistical fluctuations) as com- 
puted above, in the case of the present experi- 
ments. For the other observers it is twice the 
stated uncertainty. The length of the horizontal 
line represents half of the total energy interval 
(in those cases where the information was avail- 
able), and it is centered approximately at the 
effective value of the neutron energy. This part 
of the interval contained at least 0.8 of the total 
number of neutrons in the whole interval. Points 
marked ‘“‘B” in Fig. 4 were obtained by Bailey 
et al.,> those marked ‘‘A’’ by Ageno et al,‘ those 
marked “SR” by Salant and Ramsey,? and that 
marked ‘‘Sh’’ by Sherr.’ 

A value of o(C) for each energy interval must 
be used in the calculation of the proton cross 
section o(H) from measurements with paraffin 


' 





Fic. 3. Total neutron 
collision cross section as a func. 
tion of energy. The length of 
the vertical line at each poi 
is twice the standard on 
ation. The horizontal fine 
indicates an energy range 
which contains about 0.8 of 
the neutrons in the interval 
“B” represents Bailey et gj. 
“A” Ageno et al.,4 “sp 
Salant and Ramsey,? and “S}” 
Sherr.? 
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scatterers, and the question arises whether the 
best value at any energy is the actual result of 
the measurement at that energy, or the value 
lying on a smooth curve determined by all the 
points. The results of Bailey et al.° suggest that 
the spacing of the virtual levels of the compound 
nucleus (about 0.3 Mev between maximum and 
minimum) is such that they would not be 
separated by the present experiments because of 
the large energy intervals which were used. 
Moreover, as the energy increases, one would 
expect the width of these levels to become greater 
compared to their separation in energy. Thus it 
is probable that the apparent fluctuations in 
a(C) are mainly due to the uncertainties in the 
measurement of the transmission fractions, and 
the value of o(C) used in the calculations of 
o(H) below were obtained from the smooth curve 
shown in Fig. 3. 

For the paraffin scatterers the relation con- 
necting the cross section with the corrected 
transmission fractions is 


(o(C)+fo(H) |Ncx =log(1/T), 


where f is the ratio of H atoms to C atoms (as 
evaluated in the previous section) and Nex is 
the number of C atoms per unit area of scatterer. 
The values of «(C) used at each energy are given 
in Table IV as well as the final values of the 
proton cross sections é(H). In estimating the 
standard deviation in ¢(H), the contributions of 
f and Nex have been ignored, and the standard 
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Fic. 4. Total neutron-proton 
collision cross section as a 
function of ener te of 

bols as in Fig. 3. “Ao” 
voresents Aoki" and “ZSC” 
Zinn, Seely, and Cohen."* 


o x 1024 in em? 


% 2 4 


deviation in the values of o(C) read from the 
curve was assumed to be about 5 percent. In 
Fig. 4 the final values of ¢(H) and also the results 
of other observers (complete above 3 Mev) are 
plotted against energy. The symbols here, in- 
cluding the initials representing other observers, 
have the same meanings as in Fig. 4. “‘Ao’’ and 
“ZSC” (which did not appear before) indicate, 
respectively, Aoki’ and Zinn, Seely, and Cohen," 
who both used D—D neutrons of sharply defined 


energy. 
5. DISCUSSION 


Because of the uncertain state of theories of 
nuclear radii, a comparison with them of meas- 
ured cross sections of a single nucleus is of little 
significance. For this purpose measurements at a 
single energy on several nuclei differing widely 
in mass, such as those carried out by Sherr,* are 
much more useful. Thus the actual values of the 
carbon cross sections resulting from the present 
experiments are of interest principally because 
of the comparisons they allow, at those energies 
which have been previously used, with the 
experimental results obtained by other observers 
using different methods. As can be seen from 
Fig. 3, in which all the available data above 5 
Mev have been plotted, the results of the present 
experiments agree with those of other observers 
above 6 Mev within experimental error. 

“H. Aoki, Phys. Rev. 55, 795 (1939). 


“W. H. Zinn, S. Seely, and V. W. Cohen, Phys. Rev. 
56, 260 (1939). 
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In the previous section a theory was men- 
tioned according to which the elastic scattering 
of fast neutrons by compound nuclei may be 
regarded as essentially a diffraction of the 
neutron beam by the nuclei. This picture should 
begin to be valid at energies such that the wave- 
length of the neutron beam is small compared to 
nuclear dimensions. At 15 Mev, for example, 
X=1.2X10-" cm, and R=4.5X10-" cm. R is 
the “‘geometrical radius’ of a carbon nucleus 
computed on the assumption discussed above, 
that the cross sections for elastic and inelastic 
scattering are equal so that o(C) =2xR®. Under 
these conditions, resonances will have a rela- 
tively small effect on the total cross section and 


one would expect the latter to be essentially 


independent of energy. This expectation is borne 
out by the experimental results shown in Fig. 3, 
which indicate that the average cross section 
over 2-Mev intervals changes very little between 
15 and 25 Mev. 

The experimental results for the total neutron- 
proton cross section, and also the predictions of 
several theories, are summarized in Table V in 
such a way as to facilitate comparisons. 

Column 1 of Table V indicates the energies at 
which the cross sections have been evaluated. 
Column 2 contains values of ox, read from a 
smooth monotonically decreasing curve through 
the center of the region covered by the experi- 
mental points shown in Fig. 4. All the available 
data above 2 Mev have been considered. None 
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TABLE V. Comparison of measured proton cross sections 
with the results of various theories. 











(1) (2) (3) (4) (S) (6) (7) (8) 
exp Theoretical cross sections (10-* cm?) 

(hev) en ow TWa (sexe) TRI FRSIL RSI 
2.5 2.60 2.18 2.57 2.74 2.69 2.69 2.69 
40 191 1.56 194 2.02 2.05 2.05 2.05 
60 1.39 1.18 149 148 1.52 1.52 1.55 
80 1.08 0.95 1.22 1.16 1.19 1.20 1.32 
10.0 090 0.80 103 094 0.96 0.98 1.18 
12.0 0.78 0.69 089 0.77 080 0.83 1.08 
14.0 069 0.61 0.79 0.66 0.68 0.72 1.01 
16.0 061 O55 O71 O57 O59 0.64 0.97 
18.0 055 049 064 0.50 053 0.58 0.94 
20.0 0.50 045 0.59 044 048 0.53 0.92 
22.0 045 O41 054 040 044 048 0.91 
24.0 041 038 O.50 037 O41 O45 0.91 








of the experimental points lies off the curve by 
more than 11 percent of its ordinate. In general, 
it would be safer to compare each of the theo- 
retical curves directly with the data. However, 
since all the theoretical curves have the same 
character as the one drawn through the data, 
and since it is not necessary to distinguish 
between two theoretical curves which fit well 
over the whole range, the present method of 
comparison is probably adequate. 

In column 3 of Table V are given values of the 
cross section computed from the Wigner formula 


rh? 3 1 
| + | (1) 
ML3E+|e| 4£+]€’| 


The energies taken for the triplet and singlet 
levels of the neutron-proton system are «=2.18 
Mev, e’=0.072 Mev. The value of «’ was com- 
puted by substituting into Eq. (1) the experi- 
mental value of o for very slow neutrons (1-10 
ev) obtained by Hanstein,"® i.e., oo=2010~-* 
cm?. The values of ¢ predicted by Eq. (1) are too 
small throughout the energy range. However, 
beyond 4 Mev they come closer to the experi- 
mental values with increasing energy, and at the 
high energy end agree within 7 percent. There 
may be some significance in the fact that above 
18 Mev this expression, based on the simple 
assumption that the range of the forces is zero, 
fits the data better than any of the other theories 
except the “symmetrical” one of Rarita and 
Schwinger (column 6). Column 4 of Table V 





16H, B. Hanstein, Phys. Rev. 57, 1045 (1940). 


contains values of o obtained from the Wigner 
formula, modified by the inclusion of the first 
power of the range a of the forces. The valye 
which was assumed here is a = 2.0 10-" em. 


lf 
“ML 4E+\/e| 








en 
_ a 
sE+ |’ | 


This expression, and the experimental value of ¢ 
for very slow neutrons given in the previous 
paragraph, lead to a value of 0.069 Mev for ¢. 
The resulting values of ow. agree well with 
experiment up to about 6 Mev, but at higher 
energies the agreement is not so good as that of 
the simple Eq. (1). If in the above formula, a is 
set equal to 2.8X10-" cm, the range of force 
suggested by proton-proton scattering, the re- 
sulting values of o are about 10 percent larger 
than those in column 4, and differ from the data 
correspondingly more. 

The results of Kittel and Breit!’ are given in 
column 5 of Table V. These authors assumed a 
square potential well of radius a=2.8 X10-" cm, 
and, retaining the first four powers of this radius, 
obtained expressions for the total S-wave scat- 
tering cross section. These values are in fairly 
good agreement with the experimental results in 
the low and middle energy range, but above 14 
Mev are too small. On the assumption that the 
S-wave potential well is independent of energy, 
this difference gives an indication of the mag- 
nitude of the P-wave scattering cross section, 
i.e., about 0.03X10-*%* cm? at 16 Mev, and 
0.04 X 10-** cm? at 24 Mev. 

Kittel and Breit have calculated the P-wave 
cross section on the basis of Bethe’s neutral form 
of meson theory and obtain op wave=0.12 X10™ 
cm? for a meson mass n»=177m, and 0.01 X10™ 
cm? for 1. = 330m, as suggested by proton-proton 
scattering. The magnitude of the P-scattering 
has also been calculated by Smorodinsky'* from 
the angular distribution data of Amaldi and 
associates* by a direct method involving no 
meson theory and no assumption about the 
shape of the interaction potential. For a neutron 


17C, Kittel and G. Breit, Phys. Rev. 56, 744 (1939). 
18 J. Smorodinsky, J. Phys. U.S.S.R. 8, 219 (1944). 
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energy of 16 Mev he obtains op wave =0.01 X 1074 
cm’, in fair agreement with the observed value 
assuming Kittel and Breit’s S-wave cross section. 

The last three columns of Table V contain 
estimates accurate to about 1 percent of the 
total cross sections predicted by the three 
theories of Rarita and Schwinger” (I ‘‘symmet- 
rical,”” Il “charged,” and III “neutral’’). The 
tabulated values in each column were read from 
a smooth curve drawn through five computed 
values (at 2.8, 6, 10, 15.3, and 25 Mev) of the 
cross section predicted by the corresponding 
theory.” The experimental values agree fairly 
well with the results of the ‘‘symmetrical’”’ theory 
throughout the energy range, and very well 


19 W. Rarita and J. Schwinger, Phys. Rev. 59, 556 (1941). 
® The author is indebted to Dr. C. L. Critchfeld for 
suggestions concerning these calculations. 
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above 12 Mev. Of the theories considered, this 
one gives the best over-all fit. The ‘‘neutral”’ 
theory, on the other hand, fits so poorly as to be 
unequivocally eliminated by the data. The 
“charged” theory differs consistently and sig- 
nificantly from the experimental results, and is 
at best an unlikely possibility. 

The author wishes to express his gratitude to 
Professors H. R. Crane and G. E. Uhlenbeck for 
their interest and advice throughout the course 
of this work, and to Professor J. M. Cork, whose 
direction of the Michigan Cyclotron Laboratory 
made these experiments possible. He is also in- 
debted to Mr. Wayne Middleton and to the other 
members of the cyclotron crew for their unfailing 
cooperation. The cyclotron operating expenses 
were met by a grant from the Horace H. Rack- 
ham Fund. ° 
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The wave equations of Proca for a vector mesotron (spin 1) are specialized first to the case 
in which there is no magnetic field, and then to that in which the field is central. The wave 
functions satisfy fourth-order homogeneous differential equations. For the states /=j+1, it is 
simpler to consider two simultaneous second-order differential equations in two unknowns (see 
Eqs. (19) and (21) of this article). For a coulomb field, the form of the radial functions at r=0 
is established as u=exp[—aL/r*]r*f, where a=2(sa)i{j(j+1)]}', L-=h/Mc, M=rest mass of 
mesotron, f is an ascending power series in 7, ze is the charge on the coulomb source, a is the 
fine-structure constant e?/fc, and 8 is a constant which can be determined. 


INTRODUCTION 


ORBEN and Schwinger! and, independently, Tamm*?* have discussed the stationary states 

of a vector mesotron in a coulomb field. Corben and Schwinger have written down the explicit 

radial equations when /= j+1, and Tamm? has indicated the general nature of the solution. It is 

the purpose of the present article to give further details of the solution, together with its derivation. 
The notation of Corben and Schwinger will be adopted. 


WAVE EQUATION 


The wave equations of Proca‘ for a vector mesotron in an external field with a four-potential A. 


'H. C. Corben and J. Schwinger, Phys. Rev. 58, 953 (1940). 
?1. Tamm, Phys. Rev. 58, 952 (1940). 

*I. Tamm, Comptes Rendus U.S.S.R. 29, 551 (1940). 

*See N. Kemmer, Proc. Roy. Soc. A166, 133 (1938). 
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— ~<a Y - "4 ) a = KXaB, 
“a 4.) pet =" 8 ) Pa=KXap (1) 


——<4 40) x%= Ky’, (2) 


OX 


(xo=ct, dao = —a®, a;=a‘ (¢=1, 2, 3), Greek suffixes run from 0 to 3, and Latin suffixes from 1 to 3. 
Unless they are summation variables, 7, k, and / will denote 1, 2, and 3 in cyclic permutation.) 
Here x= Mc/h, where M is the rest mass of the mesotron. 
If we deal with stationary states, then We= —ihd¢/dt. Let us write (W+eAo)/hc=k, and assume 
that Ai=A2=A;=0. Then (1) and (2) become, if ® denotes the vector with components (¢1, ¢2, ¢) 
and g that with components g; = xix, 


—tk® —Vgo=kxo, (3) 

ikxo—curlg=x®, (4) 
curl® = xg, (5) 
divxo= —K¢go. (6) 


Combining (4) with (5) and (3), and (3) with (6), 


curl curl® =ikxx—«*®, 


=k?@ —n°@—ikV go, (7) 
V2 go t+div(ik®) = x? go. (8) 

If, now, we set igo= y, the equations become 
curl curl® — k*@+«*@ = — kV, (7a) 
™ curl curl®+x°*®=k(k®—Vo), (7b) 
V2o—K*g=divk®, (8a) 
" —xy=div(k®—Vy). (8b) 


Taking the divergence of Eq. (7a), we have 
x? div®=k div(k®—Vy)+(k®—Voy)-Vk 
= —K«*ko+(k®—Voy)-Vk, 


wi div@+ke=(1/x?)(k®—Vy) -Vk. (9) 


Then Eq. (8a) takes the form 
 W29—ny=k divO@+@-Vk 
= —k*9+(1/x*)k(R®—Ve)-VkR+®-Vk (10) 
and Eq. (7a) can be transformed, by the use of (9), into 
VO + k°®— PO =V div®+kVo 
= — gpVk+(1/x*)V[(k® —Vo) - VR]. (11) 


In the case that the field is central, Eqs. (9), (10), and (11) specialize to Eqs. (38), (37b), and 
(37a) of Corben and Schwinger. Eq. (9) [Eq. (38) of reference 1] is not a supplementary condition, 
but follows directly from the generalized Maxwell equations (3) to (6). 
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It is possible to eliminate ¢ and thus to obtain an equation which @® satisfies. Let O® =curl curl® 
—h’@+x*®. Then 


div[(1/k)O®] = —V2eg= —x2g—div(k®), curl(1/k)O@=0, 
V?[(1/k)O®]=V div[(1/k)O@] = —2Ve—V div(k®) = (x?/k)O@—V div(ke), 
7 (V?—x?)[(1/k)O@]= —V div(k@). (12) 


Equation (12) is a fourth-order homogeneous equation for ®. We shall determine the nature of 


the four solutions near r=0. 
According to Corben and Schwinger, particles with ]=j-++1 have wave functions 


®= (1, Fi(r)+7[r XLIF.A(n}P", e=G(r)P;", (13) 
where 
Fy! + (2/1) Fe!’ —Li(G+1)/r? ]Fet (k? — x?) Fe 
= (2/r?) Fi—(1/r) { Fi’ +(2/r) Fi t+ L9G +1) /r J F2t+kG}. (14) 


This follows from Eqs. (41b) and (41d) of reference 1. It may be rewritten as 


(1/r)(r Fe)” + (1/r) Fi’ + (Rk? —x?) Fo= —kG/r. (15) 
Now let 
Us=rFs and u,= Fy+u2’. (16) 
Then 
uy’ +(k?—x*)us= —RG. (17) 
Equation (44a) of reference 1 takes the form 
—kG! = (—k?+«?) Fi +[j(G+1)/r? Ju, (18) 
or 
—k(G’ —RF,) =«*F,+-[j(§+1)/r? Ju. (18a) 


Differentiating (17), we obtain 
uy!’ + (k?— Kx?) uo’ +2RR' ue = —kG’—k’G 
= (hk? —x*) (ue! — ui) +L j(G-+1) /r? Jur t+ (R’/k) [er + (kh? — x?) a2 J. 


Hence 


uy" — (k! /k)uy' + {k?—«?—[j(G+1)/r?]} ui = — (h’/k) (RPK?) U2. (19) 
Combining Eqs. (41d) and (43) of reference 1, and using (18a) 

Fi’ +(2/r) Fit CiG+1)/r]F2+kG = —(k’/x*)(G’ —R Fi) = (kh /kx®) {°F +[7(G+1)/r? Jur}. (20) 
If we now use Eqs. (16) and (17) and rearrange Eq. (20), we have 


2 ik’ a) 
———— }. 


ut'+ (==) ut {k?@—«?—[j(j+1)/r?]}u2= (--—- (21) 


k x°r2k 


Equations (19) and (21), being two simultaneous equations of the second order, would yield a 
single differential equation of the fourth order if either u; or u2 were eliminated. 
By definition, k=(W+eV)/hc. If the field is coulombian, V=ze/r, where ze is the charge on the 
source particle. Let W=wMc?. Then, if ro=e?/Mc?, 
k=w(Mc/h) + (ze*/her) = (ze*/her) {1+ (wr/zro)}, (22) 


k’=—ze*/her®; k’/k=—(1/r)[1+(wr/ero) }-. (23) 
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Solution Near the Origin 


The origin is an irregular singular point, which necessitates the introduction of an exponential 
factor in the formal solution. This factor may be found by the procedure indicated by Ince,® which 
will always give a formal series in terms of the variable r/", where N is the order of the differentia] 
equation. Application of this method leads to the factor exp[ —aL/r*], where L?=x~ and a is to 
be determined.® 














Let ui=giexpl—aL/r!] and w2=¢2 exp[—aL/r']. (24) 
Then : 
La 
= ( +=) exp[ —aL/r*], (25a) 
r 
aL 3aL a*L? 
w'=(c i" — as Se aad °) exp[ —aL/r!]. (25b) 
ri 4 ri 4r3 
Now let 
gi=aL/ri—(k’/k), (26a) 
3aL a?*L? , JG+1) | k’ aL 
QO oo fpf Lt ’ (26b) 
4 Sl? Ay r? k 2rd 
Then, from (19), (21), (25), and (26) 
¢1’ +9191’ +G2¢1= — (k’/k) (Rk? +L) g2, (27a) 
2 R j(7+1) LAR’ 
+ (a+ =)! +(o+ )e=[--—- |e. (27b) 
rl k rh 


The formalism is simplified by the substitution r=s*. The differentials transform as follows: 
D,=(1/2s)D, and D,?=(1/4s?)[D,?—(1/s)D,]. 
If w=w/nro, then (k’/k) = —(1/r)(1+r)—. Expanding, 


k’/k = —s*(1 —ws*-+-w?s*+ - - -) (28) 
Eqs. (27) become, if primes now denote differentiation with respect to s, 
gi —S—1py' +2s(aLs +s? —w+ - ++) 91’ +45"q2¢1 =4(k? +L) (1 +s?) 92, (29a) 
go" —s~1po' + 2s(aLs*+ 3s? —wt - - +) ge’ +457(g2+aLs—) oo 


=45?° 35-2? —w+---+7(f7+1)L4s*(1 —ws?+w'st—---)]gi. (29b) 
These are of the form 


¢1'+Pigi +P2¢1= Rigs, (30a) 
¢2'’ +Qige’ +Q2g2= Regi, (30b) 
where 
P,=2aLs~*+s—! —2ws+2w?s?+---, 
P,=a*L*s—aLs*+4s~[n*a*— j(j+1)]+---, (31) 


Q:=P,+(4/s); Q2=P2+(4aL/s*), 
Ri =427%a?s*+-42%a*ws-? ++4L-4(1 —ws?+w*st— +), 


R2.= (4/s*)7(G+1)L4(1 —ws*?+w'st— - - -) +12 —4ws?+--- 


SE. L. Ince, Ordinary Differential Equations (Longmans Green and Company, London, wins pp. 424-428. 
* The scalar quantity L here and below is not to be confused with the angular momentum vector L in Eq. (13). 













( 


( 








nential 
which 
rential 
2 is to 


(24) 


(25a) 


(25b) 


(26a) 


(26b) 


(27a) 


27b) 


OWS: 


(28) 


29a) 


9b) 


}0a) 
Ob) 


31) 








VECTOR MESOTRON 


Now write 


P:= > Pis**, P= > P2.s*, 


n=() n= 


V1 = > 1, ns” “s. O2= > Oe ns” ~ 
n=0 n=0 (32) 


Ri= LV Ris”, Re= DY Res", 


n=0 n=0 


@ 
gi= > oi ns"*™, and ge= ¥ v2.5", 


n=0 n=( 


For Eqs. (30a) and (30b) to be mutually consistent, since Pz and R, (also Q2 and Re) are of the same 
order in s, we must have a,=a:. The recursion formulas are then 


(e:-++2)(ar+n+1) $1, n¢2+ Lim Pi, mos, n—m43(i+n—m +3) 
+Da P», m1, n—m+4 = Lim Ri, me, n—m+4) (33a) 
(a1-+n+2)(a1+n+1) G2, n¢2+ Lim Qi, moe, n—mys(ait+n—m +3) 


+2e Qo, m2, n—m+4 = ) > Ro m1, n—m+4- (33b) 
The values of the first coefficients are 


Pyw=2aL, Pu=1, Pi=0; Po=a*l*, Pau=—aLl, P2=4[2’a?—j(j+1)]=—49°; 

Qin=Pim, MAL and Qu=5; Qom=Pom, m~1 and Q2=3aL; 

Rw=42e?, Ri=42%a°w, Ru=4L; Ro=4j(j+1)L', Re=—4j(j+1)Le, 
Ru=4j(j+1)L4w*+12. 


If n= —4, . 
P2o¢1,0=Riogeo, and Qerogeo= Reogioe, (34) 
P2Q20= RioR20, (35) 
a'=162%a?j(j+1). (36) 


If (33a) is multiplied by Reo, (33b) by P20, and the two added, the terms in ¢i,n44 and ¢ n44 will 
cancel out because of (35). We thus obtain, noting that Q4.=P3:+4 and Q2=P2:+4aL, 


(a: +n+2)(aitn+1)[Reogs, nz2t+Poo¢e, n2l+ d% (ar t+n—m+3)(Re0¢1, n—most+P20¢2, n—m+s) Pim 


+ a { (Reo¢s, n—mast+Pr0¢2, n—m+4 1) P2m— RooR im¢e, n—m44— P20Rom¢, ox~m4} 


m=] 
+4(a1+n+2)Pr0¢2, n42+4cLProge, n43=0. 
The second and third lines combine to give 


LX ¢1, n—mesL (a1 t+n—m+3)RooP imt+ReoP 2, n41—P20Re, m1] 
m=() 
+ > G2, n—mes_ (a1 tn—m+3)P oP imt+Pr0P 2, m-i— RoR, m+1 | 


m=( 


=[(a1+n+3)Pi0t+Pa JC ¢1, n+sR20+ ¢2, n+3P 20] 


+ ¥ [(ar+n—m+3)PimtPo2, m1 JL ¢1, n—m+sR20+ 2, n—m4sP20 | 


m=1 


ice , (P2oRe, m41¢1, n—m43t ReoRs, m41¢2, n—m+8)- 
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If n= —3, then with (34) we find 


4aLPog20+2Po¢20oleiPi+P2)]=0, or o=—}F, 
If n=-2, 


i(Reogiot+P20g20) —2P20¢20+4aLPo0g21+2[ —$Piit+Pe2e ]y20P20—P20R22¢10— RooRi2¢20=0. (37) 
But Re= —wR2, Ri2=R0, ReoeP 20¢10 = Ro2Rio¢20= —wReoRio¢geo and so the sum of the last two 


terms in (37) equals zero. 


Consequently 
gn /¢20= 3 —2[ —$Pirt+Pe2]=(1/4aL)(8p?+3). (38) 
Ifn=-—1, 
—$Progi0t+ Pr¢11+P21¢10 = Rioge1 = Riogeo— (1/4aL) (8p? +3), 
oo (6949) 8 
aig onsen — 39 
oe 4cL . a om 


One can thus obtain all the gen and gim in terms of g29 (or ¢10). 
Summarizing, 


1 2 
imsheu| 1+|—(6o'+8)+—] ++ ven } (40a) 
4aL cL 


1 i 
em shen| 1-+——(8p?+3)s+::: | (40b) 
4aL 


¢20 = (a°L*/42*a") g10, where a is given by (36). 
From (16) and (24) 


Fy=(1/s*)e*”'*s4{ poot+gaist+---}, (41) 
1 duz 1 

Fy =u —— —=—e "sf pot gist: -} — | axs-te-sti got past: :} 
2sds_ s? 2s 


peehl*s-4{ FS poo t+ dyes t+: 1 (42) 


That is, F2 is of order r~*/4 exp(—aL/r*) and F; of order r~*/4 exp(—aL/r*) at the origin.’ 

The wave equation then has solutions of the form displayed in Eqs. (13), where F; and F, are 
given by Eqs. (42) and (41), respectively, and s=r!. Equation (17) may be used to compute G. 

For the more general case where the field is non-coulombian, one must solve Eqs. (19) and (21), 
using the appropriate dependence of k on the radius. The procedure is as given in Ince if k = O(1/r‘) 
as r—0, where 7 is an integer, and should not be particularly difficult even if this is not so. 

The solutions which we have found above are expressed as power series near the origin, and it is 
not to be expected that convergence is rapid away from the origin. For practical computations, a 
knowledge of the solutions for all values of r is desirable. It seems that this is to be attained only 
with a machine such as the ENIAC or EDVAC, at least if one is concerned with the scattering 
problem,’ which usually requires phase shifts to be found for many values of j. 


7 The exponential factor agrees with that found by Tamm, references (2) and (3). 
8 The scattering problem will have a solution if, at large r, one can find a stationary state which is a sum of a plane 
wave plus an outgoing wave. This question needs further study. 
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If an “f-interaction” is introduced into the Hamiltonian for two protons, the square well of 
width e?/mc* representing the specific nuclear proton-proton interaction at low energies must 
be modified to compensate for the effect of the “f-interaction.” It is found that the experi- 
mental dependence of phase shift on energy for low energy proton-proton S-scattering would 
require the assumption of different well depths for different energies of the incident protons. 
The necessary changes in the potential would be sufficiently large to destroy the symmetry 
between the specific nuclear proton-proton and proton-neutron 'S interactions. 





1, INTRODUCTION 


HE ‘“‘f-interaction’’ proposed by Pais! in 
order to overcome some of the divergence 
difficulties of quantum electrodynamics would 
add to the coulomb interaction potential between 
two protons, a term V; giving: 


Veoulomb + V;=(e/r)[1—2 exp(—x«r/ro) ], (1) 


where ro=e?/mc*, and the range 1/x must be 
about 2 to account for the mass difference 
between neutron and proton. A more exact 
evaluation gives? 1/x= 1.82. 

It has been shown by Breit and collaborators*~* 
that the experimental phase shift vs. energy 
curve for proton-proton scattering at energies 
between 0.2 and 2.6 Mev is best accounted for by 
taking the specifically nuclear interaction poten- 
tial, which must be added to the coulomb repul- 
sion as a square well of width ro and depth 

¢=11.35 Mev, and that this corresponds almost 
exactly to the square well potential for the 'S 
proton-neutron interaction in the same energy 
range. A preliminary estimate by Breit? indicated 
that the change 6D/ in the well depth necessary 


* Assisted by Contract N6ori-97, Task I and Contract 
N6ori-44 Task XVI of the Office of Naval Research. 

1A. Pais, Phys. Rev. 68, 227 (1945). 

* Unpublished calculation by G. Breit. 

*G. Breit, H. M. Thaxton, and L. Eisenbud, Phys. Rev. 
55, 1018 (1939). Referred to as BTE in the text. 

‘G. Breit, E. U. Condon, and R. D. Present, Phys. Rev. 
50, 825 (1936). 

'G. Breit and J. R. Stehn, Phys. Rev. 52, 396 (1937). 
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to compensate for the effect of the ‘‘f-interaction”’ 
on the proton-proton S-wave phase shift at low 
energies would be sufficiently large to destroy the 
symmetry between the specifically nuclear 
proton-proton and proton-neutron 4S interac- 
tions. The present work investigates the mag- 
nitude of 6D’ as a function of E and 1/« for 
0.2 Mev=E=2.6 Mev, and 1=1/«35S3.5. 


2. PROCEDURE AND RESULTS 
(a) First-Order Perturbation Estimate of 5D/ 


A first order perturbation calculation for the 
phase shift Ko produced by a perturbing central 
field potential 5V on a radial wave function of 
relative motion Fo/r gives® 


Kos— f (6V/E’) Fo*dp, (2) 


where E’ =energy of relative motion; p=2z7r/A; 
A=wave-length of relative motion. To the ap- 
proximation of Eq. (2) the phase shift produced 
by a square well perturbation of a given width is 
proportional to the depth 5Do of the perturba- 
tion. 

For proton energies up to ~3 Mev the effect 
on the phase shift produced by the omission of 
the coulomb potential at distances smaller than 
ro can be compensated* by changing the well 
depth by 6D*= —0.83 Mev. Hence the first-order 


*H. M. Taylor, Proc. Roy. Soc. A134, 103 (1931). 








TABLE I. Values of 5Do/ in Mev. 
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1/x E=0.2 Mev E=14 Mev E=2.6 Mev 
1.0 —1,70+0.1 —1.51+0.02 — 1.46+0.01 
1.5 —2.7 +0.3 —2.1 +0.1 —2.0 +0.1 
2.0 —3.4 +0.3 —2.6 +0.2 —2.4 +0.2 
as —4.0 +0.4 —3.0 +0.2 —2.7 +0.2 
3.0 —4.6 +0.3 —3.2 +0.2 —2.9 +0.2 








estimate 6D,’ of 5D/ is 


5D = —0.89| ff 2Fe exp(—ar/rolr/r| / 


J “Pedr/r] (3) 


Here Fo/r is the regular S-wave function for the 
unperturbed potential 
Vo(r) = —[D*+é6D*]=—10.5 Mev; 
=e /r; 


r<T1o 
r>To. 


For r<ro, Fo is a sine wave. For r>ro two ap- 


proximations, 
Fo ~at + bp, ' (a) 


Fo~A sin(p+6), (8) 


were used in the evaluation of (3). Except for 
E<1 Mev, when d?F)/dr?>0 over a considerable 
region of r>ro, the actual value of Fo (coulomb 
function) lies between (a) and (8), which thus 
give an upper and lower bound for 6D ’. For low 
energies (a), and for high energies (8) is the 
better approximation. Because of the rapid 
decrease of exp(—x«r/ro) with r for r>7ro, at least 
for small 1/x, the inaccuracy introduced by the 
use of (a) and (8) is small compared with the 
inaccuracy of the perturbation method itself. 
Equation (3) gives the following values of 5Do’: 
The uncertainties indicated in the table arise 
because of the use of approximations (a) and 
(8) instead of the exact coulomb function at 
r>ro (cf. appendix). 


(b) Phase Shift Calculation for 
V= (e?/r)[1 —3 exp(—x«r/To) | aac i : 
and Improved Estimate of 5D’ 


In order to improve the estimate of 6D/ 
obtained in (a) for 1/x=2, the most likely value, 
phase shifts caused by 6 V were calculated for the 





regular solutions of the equation 


PY /dr’+(M/h*)(E’ —e/r—6 VF =0, (4) 
with 
6V=—2 exp(—2r/ro)—D!**spprox, (4) 


for E=0.4 and 2.0 Mev. The well D/*+,,,... is 
supposed to represent approximately the specific 
nuclear potential which must be added to the 
coulomb and ‘‘f-interactions’”’ to give the experi- 
mental proton-proton scattering phase shifts, 
It was chosen as 


D!*¢,jprox = 7.8 Mev; r<ro 
=0; r>To. 


This corresponds to an assumed value of §D/ 
=7.8—11.35 = —3.55 Mev, which is larger than 
the values given in Table I, since 6Do/ will, in 
general, underestimate 45D’ in the present case, 

The values of § were obtained by numerical 
integration up to r=3e?/mc*, at which point the 
phase shifts Ko(r=3e?/mc*) were calculated 
from Eq. (7.8) of Breit, Condon, and Present,‘ 
the values of the coulomb functions occurring in 
this equation having been obtained from a table 
in BTE.* 

At r=3e?/mc* the values of § were normalized 
to FocosKo+GosinKo (Fo, Go=regular and 
irregular coulomb function, respectively), and 
the small additional phase shifts 6Ko which were 
due to the ‘“‘tail” of the f-field potential beyond 
r=3e?/mc? were estimated by Eq. (2) with 
Fy=%, 6V = —(2e/r) exp(—2r/ro), the integral 
having been taken from r=3e?/mc? to ~; two 
approximations were used for §, giving upper 
and lower bounds for 6Ko. The perturbation 
method was used because the quantities entering 
into the exact equation for the phase shift are 
not readily available for r>3e?/mc?. Its use is 
justified by the smallness of the values of 6K» 
obtained (cf. appendix). The additional cor- 
rections 6D in the well depth required to bring 
the phase shift calculated by this procedure into 
agreement with the experimental values were 
estimated by the first order perturbation for- 
mula: 


Kalcale.) —Kolexp)— f (6D/E')Be. (5) 


One thus obtains corrected values for the well 
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depth, Dift+¢= D!**,5r0x + 6D, and for 6D/ = D/te 
—D*. The results are summarized in Table II. 


(c) Phase Shift Calculation for 
V=(e’/r)[1—2 exp(—xr/r9) ] 
Without Square Well 


The results so far obtained show that within 
the range of values for 1/« considered physically 
admissible, the contribution to the phase shift 
caused by the ‘“‘f-interaction’’ becomes com- 
parable for low energies to the total observed 
phase shift because of a near resonance of the 
“f-field” wave functions for large 1/x. A rough 
estimate of the value of 1/« which at a low 
energy would give the total observed phase shift 
from the ‘‘f-interaction’’ alone, without any 
square well, can be obtained by putting FoCop 
(C)=normalization constant for regular coulomb 
function) and then estimating the phase shift, 
Ky, due to the “‘f-interaction” by use of Eq. (2), 
with 6V=—2.exp(—x«r/ro). For E=0.2 Mev, 
C?=0.271, and one has, carrying out the inte- 
gration, 

Ko 0.00735x~. 


A factor 4/3 approximately corrects for the 
fact that the normalization constant Co is too 
low, since the attractive ‘“‘f-field” potential 
increases the wave function at small p. For 
E=0.2 Mev the experimental value of the phase 
shift is Ko=6.8°0.12 radians, so that one has 
finally 


x-*&(0.12/0.00735)(3/4);  1/«2¥3.5. 


Calculation, by the method of section (a), of 
the phase shifts produced by an ‘‘f-field’’ poten- 
tial with 1/x=3.5 alone superimposed on the 
proton-proton coulomb potential for E=0.2 Mev 
and E=2.0 Mev gives: 


Kof =0.1225 40.0065 
=7.02°+0.37°; E=0.2 Mev (6) 
Ko =0.2790 +0.0003 
=14.82°+0.01°; E=2:0 Mev. 


The indicated uncertainties are again due to the 
use of approximations for § at r>3e?/mc*. (cf 
appendix). The experimental values are: 


E=0.2 Mev 
E=2.0 Mev. 


1/x=3.5 


Ky=6.8°; 
Ky=45°; 
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TABLE IT. Phase shifts and corrected well depths 











for 1/x=2. 
Ko, experi- 
E mental 6D Dire 6D! 
(Mev) Ko* (from BTE) (Mev) (Mev) (Mev) 
0.4 0.275+0.001 
=15.8° 15,2° —0.17 7.63 —3.72 
2.0 0.638 +0.000 
236.55° 45° 1.12 8.92 —2.43 








* Calculated for the solutions of Eq. (4). 


3. INTERPRETATION OF THE RESULTS AND 


CONCLUSIONS 

The results show that the phase shift con- 
tribution caused by the ‘‘f-field”’ increases much 
more slowly with energy than the observed 
proton-proton scattering phase shifts at low 
energies. For 1/x=3.5 the total observed phase 
shift at E=0.2 Mev is accounted for by the 
“f-field,”” but at E=2.0 Mev only 1/3 of the 
observed phase shift can be attributed to it as 
is seen from Eq. (6). 

If the “f-field” with the lower value of 1/x=2 
is combined with a square well of width e?/mc’, 
the well depth must be reduced from the usual 
value D¢=11.35 Mev by several Mev and must 
be made velocity dependent to account for the 
observed phase shifts (Table II). The first-order 
estimates of 6D’ (Table 1) are approximately 
correct for 1/x=2, but they fail to show the 
large increase of the “‘f-field” phase shift at low 
energies for larger 1/x, which is due to a near 
resonance of the “f-field’” wave function. Even 
for 1/x=2 the first-order perturbation values 
underestimate the “‘f-field’” contribution to the 
phase shift at low energies. This is due to the 
fact that Eq. (3) underestimates the contribution 
of the “tail” of the “‘f-field’’ beyond ro, since the 
unperturbed wave functions Fy are concentrated 
at small p because of the large depth of the 
unperturbed well, D¢= 11.35 Mev. This “‘tail’’ is 
important at low energies. 

In summary, the calculations show that the 
“f-field”’ potential would have to be combined 
with a velocity dependent square well of width 
e?/mc to account for the observed proton-proton 
scattering. The variation of well depth with 
energy might be avoided by reducing the range 
of the well and correspondingly increasing its 
depth. In either case, the symmetry between the 
specific nuclear proton-proton and _ proton- 
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neutron 1S potentials would have to be aban- 
doned if an “‘f-interaction” were introduced into 
the proton-proton potential. 

We would like to acknowledge our indebted- 
ness to Professor Breit whose preliminary 
estimates first indicated the nature of the effects 
discussed in this paper and who helped with 
suggestions and discussion throughout the work. 


APPENDIX 


Remarks on the Use of Approximate Wave 
Functions in the First-Order Estimates by 


Eq. (2) 


In (a), approximations (a) and (8) for Fo at 
r>ro were used in the evaluation of Eq. (3). 
These approximations correspond to the re- 
placing of the coulomb repulsion e?/r in the dif- 
ferential equation for Fo by a constant potential 
Vo. The solution for Vo=£’ is (a); for Vo=0 it 
is (8). A better approximation for the purpose 
of evaluating the right side of Eq. (3) would be 
obtained by letting V» have such a value that 
the effective kinetic energy E’—Vo should 
approximate the average of the actual kinetic 
energy, in the region which matters for the 
integrals which occur in the perturbation formula 
Eq. (3). This Vo will be referred to as (e?/r) 
because E’ —e?/r is the kinetic energy. 

If 5Do/(a@) and 5D,‘(8) are the values obtained 
by the use of approximations (a) and (8), respec- 
tively, a better approximation for 5Do/ is then 
obtained by taking their weighted average in the 
form 


1 
éDf =! (e?/r) 6D! (a) + (E’ — (e?/r)) 8D (8)}. 


Since about half of the integral from ry to © in 


Eq. (3) is contributed by the region r,<, 
<ro(i+1/x), e/[ro(it+1i/x)] gives the right 
order of magnitude for (e/r). The values of 
éDo’ given in Table I correspond to using this 
value. Since the differences 5Do/(a) —5D,/(8) are 
no larger than the inaccuracy of the first-order 
perturbation method itself, the exact value of 
(e/r) to be used is not important. The uncer. 
tainties indicated in Table I are the smaller of 
the differences 


| 6D of — 5D f(a) | and | 6D — 5D‘ (B) | 


in each case. 

In (b) and (c) two approximations were used 
for § in the evaluation of 5Ko(r>3ro) by Eq. (2): 
(1) approximate coulomb function, (2) approxi- 
mate coulomb+“‘f-field” function. 

Parabolic approximations were used, with the 
coefficient of p* determined in each case by taking 
an ‘‘average value” of the appropriate potential 
in the sense explained above. The exact nature 
of the approximations used is unimportant in 
this case, since r~'exp(—x«r/ro) is very small 
when r — 379 is sufficiently large for the differences 
between approximations (1) and (2). to become 
important. 

The contributions to the total calculated 
phase shifts Ko* of 5Ko* as compared to K,* 
(r =3e?/mc*) in Table II are: 


E (Mev) Ko*(r = 3e*/me*) 5Ko* 
0.4 0.2068 0.0685 +0.0010 
2.0 0.6069 0.0313 +0.0000 


The values of 5Ko* given are the averages 
of the values obtained by using approximations 
(1) and (2). The indicated uncertainties are the 
differences between the averages and _ these 
values. They are seen to be negligible. 
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The Physical Significance of Birkhoff’s Gravitational Equations 


HERBERT E. Ives 
Bell Telephone Laboratories, New York, New York 


(Received December 11, 1946) 


Birkhoff’s gravitational equations are put in terms of dt in place of the local time ds used 
by him. The transformed equations show that Lorentzian mass has been used, and to the 
Newtonian attractive force is added a force normal to the direction of motion, v*/c* times the 
component of the gravitational force normal to the motion. 





1. INTRODUCTION d*x GMx 


can | dx\? (2 : 
HE theory of gravitation proposed by the ds? 7 r3 rc? (5) si x) 








late G. D. Birkhoff was developed by him 
from a very general mathematical standpoint, 
and the physical significance of the theory is 
difficult to grasp in the form in which it has been 
presented. As indicated by the title of his most 
detailed presentation, ‘“‘El Concepto Matematico 
de Tiempo y la Gravitacion,”! Birkhoff laid 
stress on the use of ds (local time) in place of dt 
(Newtonian time); the “forces” figuring in his 
development are Minkowski forces instead of 
physical (Lorentzian) forces. Yet his solution for 
a planetary orbit (which gives the correct 
advance of perihelion of Mercury) contains no 
terms in ds, raising the question whether this 
result does in fact depend on his adherence to 
this variable. Actuated by the desire to see what 
form his ‘‘force’’ equations would take if ex- 
pressed in terms of d¢ and ordinary physical 
force (time rate of change of momentum), I have 
carried through the necessary transformations, 
with the results given below which exhibit some 
features of interest. 


2. THEORY 


Birkhoff’s gravitational force, which is the 
third term of an expansion which he states? “‘it 
is natural to set”’ as typifying ‘‘all force vectors,”’ 
gives, in his hands, for the motion of a particle 
moving in a plane in the gravitational field of a 
mass M, 





1G. D. Birkhoff, Boletin de la Sociedad Matematica 
Mexicana 1, 1 (1944). 

*A. Barajas, G. D. Birkhoff, C. Graef and M. S. Val- 
larta, Phys. Rev. 66, 142 (1944). 


GM dx dr 


————-—, (1 


’ 
r°c? ds ds 


d*y GMy GMy dx\* /dy\? 
ee ) +z) | 
s 


ds? r3 rc? L ds 





GM dy dr 


——-—, (2 


r¢2 ds ds 


It is these equations that we wish to put in 


-terms of dt. To transform the right-hand side we 


use the relations 
ds =dt(1—v?/c*)t, v?=%*+y’, 


where the dots indicate differentiation with 
respect to ¢. We note that the left-hand side is 
the Minkowski force, 


F,(1 —vy*/c?)-4 


in which F, is the Lorentzian force, that is, the 
time variation of momentum where mass varies 
with velocity according to the relation 


m = mo(1 —v?/c?)-4, 
We then have 
(Fx)2=L—(GMmox/r*)(1+0*/c*) 
+ (GMmot?/r*c*) (1 —v?/c?)-4, (3) 
(Fx)y=C—(GMmoy/r*) (1 +0?/c*) 
+ (GMmoyr/r?c?) \(1—v?/c?)-4. (4) 
Putting r*=xi+yy we get 


GMmox (yx —<y) 


Faden ———] +4 
_ r r 





Jo —v*/c*)-', (5) 
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GMmfy (yx—<azy) 
(Fi),= -——| --a—— 
rf? Lr rc? 


|a-e/er. (6) 


Comparing these with the Newtonian forces 
F,= —GMmyx/r', (7) 
y= —GMmy/r', (8) 


it appears that the Birkhoff equations differ in 
the use of the Lorentzian mass, and the addition 
of terms in (yx—Zy). 
It is instructive to put these equations in 
polar coordinates through the relations 
yx — xy = 76 
and 7 R 


Fr=(F.x+Fyy)/r, Fo=(—Fiyt+F,x)/r 


used: on the right-hand side, using on the left the 
statement of the Lorentzian force as given by 
Eddington,’ yielding 


mo, (#—76*)[1 — (7? +776?) /c? “3 


d 
+1 — (F776) /c2}3 r 
t 


= —(GMmy/r?)(1+176?/c?) 
x[1—(# +1767) /c?', (9) 


me| (276-+r6)[1 — (#2+7°62) /c?}-3 


d 
+O +re)/c}| 
t 


= (GMmy/r?) (7r6/c?) 
X [1 — (7? +776?) /c?}-4. (10) 
These equations are solved by multiplying (9) 
by 7, and (10) by r@, and adding, which gives 
<1 — (i +1776") /e?] 
dt 2GMr 
[1 —(#°+776") /c?] ‘i re 





from which 


[1 — (#2+176) /c?]=exp[—(2GM/rc?) +k]. (12) 


3A. S. Eddington, Phil. Mag. [6] 34, 321 (1917). 
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Using this result in (10) gives the areal con- 
stant 


r°6=h exp[ —(2GM/rc*) +k], (13) 


leading to the solution, replacing r by 1/x 


d? 


u 
ao tan (c*/h*)[exp(k+2GMu/c*) —1] 


Xexp(k+2GMu/c?) (14) 


which is identically Birkhoff’s expression, from 
which the advance of perihelion of Mercury is 
correctly indicated. 

Equations (9) and (10) may be put in simpler, 
although less instructive form, by inserting the 
value of [1—(i+7°@) /c?]! from (12), giving 





GM [ 72 — 726? 
7¥—r@@= — 1— | (15) 
r? L c 
d(r?6) /dt = 2G M7r6/c?. (16) 


These are identically the expressions given by 
Fernandez‘ who also puts his equations in terms 
of dt instead of ds, in discussing the problem of 
two bodies in Birkhoff’s theory. 


3. DISCUSSION 


Examination of (9) and (10) shows that 
Birkhoff’s equations for a planetary orbit are 
the equations which one would obtain by using 
Lorentzian masses throughout in place of the 
invariant masses of the simple Newtonian theory, 
with the addition of terms in 7, 76, and c? on the 
gravitational side of the equations. 

First consider the appearance of the Lorent- 
zian masses. Birkhoff in his presentation lays 
great stress on the use of local time. The idea 
that a planet is controlled in its course by the 
time indicated on a clock carried on it lacks sub- 
stantiality. Physically, forces act on masses. 
Actually the factor by which local time is dis- 
tinguished from absolute time is the “‘contraction 
factor,”’ (1—v?/c?)!. Now this is also the factor, 
appropriately placed, by which local mass is dis- 
tinguished from stationary mass. Our trans- 
formed equations are force equations in terms of 
Lorentzian local mass. Hence although Birkhofl 


4C. G. Fernandez, Boletin de la Sociedad Matematica 
Mexicana 1, 25 (1944), especially p. 36. 
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really stressing the importance of the contraction 
factor,» and his development could just as well, 
but with real physical significance, be founded 
on local mass. 

Next consider the added terms on the right 
in Eqs. (5) and (6). Their characteristic is 
that when multiplied respectively by # and y 
and added they cancel exactly, leaving us with 
the simple Lorentzian Eq. (11) which may be 
written 


—(1—v*/c2) 
dt 





GMm, 
——mo + =U. 
2 (1—v?/c?)! = r%c?(1—v?/c?)! 


This equation states that the sum of the deriva- 
tives of kinetic and potential energies is zero 
which may be accepted as a fundamental re- 
quirement. An infinite number of such additional 
terms could be set down, which cancelling out, 
would leave the energy equation unchanged. 
Any term of this sort has, however, the result of 
changing the value of r*6, which is obtained by 
integration of the second of the equations in 
polar coordinates (10), alone. The whole difference 
between the Birkhoff solution (14) and the 
simple Lorentzian solution comes about from the 
value of r?6, the “areal constant,’’ which occurs 
(squared) in the denominator of (14), the nu- 
merator having always the same fixed value 
determined by (12). The terms added by Birk- 
hoff in his force equations change this from 


r6=h exp[4k—GM/rc*], 
the Lorentzian value, to 
r°6=h exp[k—2GM/rc*], 


thereby changing the predicted advance of 
perihelion from 4, the value obtained from the 
Newton-Lorentz solution, to the full observed 
value. The prediction of these particular added 
terms is, in the light of this analysis, the key 
contribution of Birkhoff’s theory. 

Let us study these added terms more in detail. 
They are completely described by the following 


‘The contraction factors are the coefficients of mass, 
length, and clock rate, which are demanded to insure the 
conservation of ener, and momentum in radiation- 
(1948), interactions. Cf. H. E. Ives, Phil. Mag. [7] 36, 392 
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talks of local time as the significant factor, he is 
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Fic. 1. Resolution of the force on a particle. 


statement: A particle in motion in a gravitational 
field experiences an additional force, normal to the 
direction of motion, of value v*/c* times the com- 
ponent of the gravitational force normal to the 
direction of motion. 

The proof of this theorem is exhibited in the 
accompanying Fig. 1, in which 

fy =gravitational force normal to direction of 
motion— 


frrO(¥? +1°6)-4; 


fu =force due to motion, normal to direction 
of motion— . 


fv? /c? =fp(r6/c*) (#2 +7°6*)!. 


The component of added force in the radial 
direction is 


fur6(i?+126?)-3 = fer26?/c?, 


and the component of added force in the direc- 
tion normal to the radius is 


— fui (7? +r°6?)-4 = —frir6/c?. 


Putting —(GMmp)/r)[1—(#+r6@)/e}-? for 
fe we get the added terms in (9) and (10). 

Denoting the Lorentzian forces (time vari- 
ation of momentum, mass varying with velocity) 
by F,, and the Newton-Lorentz gravitational 
force —(GMmp)/r)[1—(F#+r@)/2}* by Gz, 
the Birkhoff gravitational equations for a 
planetary orbit may be concisely expressed as 


follows 
(Fr)r=(Git+(Gr)nv*/c* Je, (17) 
(Fr)e=(Gi+(Gr)nv*/c*)s, (18) 
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where the subscript N denotes components 
normal to the path of the particle, and the sub- 
scripts R and @ denote components in the radial 
direction and perpendicular to the radius, respec- 
tively. 

As to the origin or nature of this transverse 
force produced by motion, it is of interest to 
observe that it is similar to the ‘‘fundamental 
law”’ proposed in a posthumous note by Gauss® 
for the mutual action of two elements of elec- 
tricity in relative motion. The occurrence of ¢* 
stems from the idea of the attraction being 
transmitted with the speed of light. 

The Birkhoff force equations for a planetary 
orbit can be summarized, according to this anal- 
ysis, as follows: They are the equations one 


°C. F. Gauss, Werke (Gottingen, 1863-74), Vol. 5, p. 616. 
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would obtain from a Newtonian attractive force 
acting on the Lorentzian local mass of the planet, 
with the addition of forces caused by motion, 
transverse to the path of the planet, which do 
not affect the conservation of energy but alter 
the areal constant. The precise value of these 
added forces and the method of obtaining them 
is thus of crucial importance. Apparently these 
forces are not introduced by Birkhoff with con. 
scious resort to physical concepts, but they are 
present because of his choice for gravity of the 
third term? of a formal expansion in rational and 
integral components of a typical force function, 
in which successive terms are of increasing com- 
plexity and hence provide for additional force 
components. An independent physical derivation 
of these transverse forces would be welcome. 
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Refraction of Plane Non-Uniform Electromagnetic Waves between Absorbing Media 


L. PINCHERLE 
King's College, University of London, London, England 
(Received August 6, 1946) 


When a plane non-uniform electromagnetic wave is refracted between two conducting 
media, there are two possible positions for the propagation vector in the second medium. 
Consideration of the energy flow shows that each solution holds within a certain range of 
values of the (complex) angle of incidence, the transition from one to the other occurring in a 
discontinuous way. The two cases of the electric vector perpendicular and parallel to the plane 


of incidence are discussed. 


(1) INTRODUCTION 


HE problem of the refraction of a plane 
non-uniform electromagnetic wave at the 
plane boundary between two conducting media 
is not generally fully discussed in textbooks where 
it is pointed out that, with the use of complex 
angles of incidence and refraction and of complex 
propagation vectors, the problem is formally 
identical to the usual one in which perfect dielec- 
trics are involved.! 

It is the purpose of this paper to complete this 
treatment discussing the new physical features 
which appear when both media are conducting, 

1 See, for instance, J. A. Stratton, Electromagnetic Theory 


(McGraw-Hill Book Company, Inc., New York, 1941), 
pp. 500-524. 


in particular, a discontinuity occurring in the 
(complex) propagation vector in the second 
medium at the (complex) angle of incidence for 
which there is no average flow of energy across 
the boundary. 
(2) THE PROPAGATION VECTOR IN THE SECOND 
MEDIUM 

Let the boundary be the plane y—z, and the 
plane of incidence the plane x —2, the x-axis being 
directed from the first medium into the second. 


Let any field component be represented by 
Ee-ki:ttiot (1) 
with ; 
ki=a+ib; k.=A+7B, (2) 


E being a complex amplitude. 
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Given k, and the constants of the second 
medium, k, is determined by: 

(a) The relations holding generally for non- 
yniform waves? 


ko? =A?—B?+2iA-B=wpe.t+iwuor, (3) 


where all letters have their usual meaning, and yu 
will be assumed to be the same for both media; 

(b) The boundary conditions, which give 
(Snell’s law) : 


k; sin8, = ke sin, (4) 


#;, 0: being the (complex) angles of incidence and 
refraction. Equation (4) is equivalent to the two 
relations between real quantities: 


a,=A,; b,=B,: (5) 


a, and b, can be expressed as functions of the 
constants of the first medium and of 6; it is 
convenient in the following to consider them as 
the independent variables in terms of which to 
express the propagation vector in‘ the second 
medium. 

Conditions (3) and (5) are not sufficient to 
determine this vector in a unique way. In fact, if 


2A ,B,=2a,b,>wyoe, (6) 


the product A,B, is negative, but, without con- 
sideration of the energy flow, we cannot imme- 
diately decide which of the two factors is nega- 
tive. Only when A,B, is positive, physical 
reasons will immediately rule out the solution 
with both quantities negative. Even in the case 
of perfect dielectrics it is possible formally to 
construct a second solution for the field in the 
second medium, satisfying Maxwell’s equations 
and the boundary conditions and corresponding 
to a “refracted wave’’ moving towards the 
boundary and on the same side of the normal as 
the incident wave: this solution is absurd from 
a physical point of view. It is an example of the 
striking way in which the presence of a conduc- 
tivity in both media modifies these phenomena, 
that in our case a decision between the two 
possible solutions is not so immediate, and that, 
as we shall see, the requirement that the energy 
should flow from medium 1 to medium 2 leads 
to the necessity of choosing one solution within a 





* See, for instance, J. A. Stratton, reference 1, p. 341. 
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certain range of the complex angle of incidence 
and the other within another range, the transi- 
tion occurring in a discontinuous way. 


(3) THE ENERGY FLOW 


(a) Electric vector Perpendicular to the Plane 
of Incidence (E,) 


If E, E’, E” are the complex amplitudes of the 
electric vector in the incident, reflected, and 
transmitted waves, respectively, the average 
flow of energy perpendicular to the boundary in 
the first medium is given by: 


(Sz )w = Re($E,H,*) 
= (e-%2* /2uw)[ | E| 2b,e-%=* — | E’ |*b,e%=* 
—2a,|E| | E’| sin(2b.«+~7)], 
(x<0) (7) 
y =arg (E'/E). (8) 


with 


The first two terms may be interpreted as the 
energy flow associated with the incident and 
reflected waves, respectively. The third term,’ 
which vanishes for ¢,=0, is caused by inter- 
ference between the incident and reflected waves : 
there is a local radiation from the maxima to the 
minima of the standing waves system. 

The average flow of energy perpendicular to 
the boundary in the second medium is: 


(Se m= (1/2ya)er*aartee)| 2” |B,, 
(x>0) (9) 


so that the condition that S, be continuous at 
x =0 is expressed by: 


6,1 —|E’/E|*]—2a,Im(E’/E) 
=B,|E"/E\|*. (10) 


This condition is automatically satisfied by the 
expressions for E’ and E”’, whatever the sign of 
B,; however the physical conditions of our 
problem impose that the energy flow be directed 
from medium 1 to medium 2, therefore the 
additional relation for the unique determination 
of ke is in this case 


B,20. (11) 
*M. Born and E. Ladenburg, Physik Zeits. 12, 198 
(1911). 
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%a, 


Fig. 1. E perpendicular to the plane of incidence. Discon- 
tinuity along AO. 


(b) Electric Vector in the Plane of Incidence (£,,) 


In this case the condition of continuity of S, 
at the boundary gives, if we indicate by H, H’, 
H" , the amplitudes of the magnetic vector in the 
incident, reflected, and transmitted waves, re- 


spectively : 


a,01+ brew 
—__—_—(|H|*+ |’ |*) 
34- 


01 €1 


b.01 — Az€1W 
o+ €1:°w? 
A,o2+ Brew 
mailed 
o2?+ €27w? 


with 


| H| | H’|siny 


y=arg(H’/H), (12) 


and the condition that the energy must flow in 
the right direction is expressed by 


A,o2+B,we:> 0. (13) 
(4) THE DISCONTINUITY FOR THE CASE E, 


Writing explicitly the equations determining 
A, and B, (see (3) and (5)): 


B?Z—A,? = eq4uw’? —b,?+4,’, (14) 
2A,B,=o2w —2b4,, (15) 


we see that either A, or B, is zero when 


2a,b, = 72, (16) 
and that, when (16) is satisfied : 
B, is zero if 6,22 e2uw*+a,’, (17) 


Az is zero if b,2& é2uw*+a,’. 
In the case of small conductivities (¢<ew) 


b,—(e14) 'w sin gi, (18) 
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¢1 being the real angle of incidence, and a,—0 
so that Eq. (17) becomes 


singi=n, 


with n=(é2/e,)}. (19) 


Although the notion of critical angle loses all 
meaning in the case of conducting media, we 
shall for brevity say that, when (16) is satisfied, 
B, vanishes in the region of total reflection, and 
A, in the region of ordinary refraction. 

If now we plot the curves (Fig. 1) 


b,?—a,? = equw* (line a), (20) 
2a.b,=uwo2 (line y), (21) 


and divide the plane a,v, into four regions (no 
loss of generality is entailed in excluding negative 
values of 6,), Bz vanishes on the boundary 
between regions II and IV, and A, on the 
boundary between | and III. The product A,B, 
is negative in regions III and IV. As B, must 
always be positive (Eq. 11), and no change of 
sign of either factor occurs in crossing line « 
(Eq. 14), it follows that A, is positive every- 
where left and negative everywhere right of line 
y. In the region of ordinary refraction the transi- 
tion from positive to negative A, occurs continu- 
ously, passing through A ,=0, but in the region of 
total reflection A, does not vanish on line y; 
therefore it must be there discontinuous. Indeed 
it jumps from the value +(b,—a,?— ey}! on 
the left to the value —[b,?—a,? —e2yw* }' on the 
right—namely, the planes of equal amplitude in 
the second medium, perpendicular to A, have, 
on the two sides of the discontinuity, symmetrical 
positions with respect to the boundary. 

To the two different orientations of k» cor- 
respond different expressions for the field com- 
ponents, and also for these the transition from 
one set of values to the other occurs discon- 
tinuously when B,=0. Experimental verification 
of these considerations could be better obtained 
by investigation of this last phenomenon. 

A negative value of A, means that the am- 
plitude increases exponentially for x— ~ ; but 
in an actual case the cross section of the beam 
is not infinite, and then the amplitude decreases 
in penetrating into the second medium, as follows 
from A-B>0. 

A difficulty arises from the fact that at the 
discontinuity (B,=0) there is no average flow 
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of energy across the boundary, and it is not 
apparent how the field in the second (dissipative) 
medium can be maintained. This difficulty is 
removed when a beam of finite cross section is 
considered, for then edge effects allow for some 
energy crossing the boundary even when B,=0.‘ 


(5) THE DISCONTINUITY FOR THE CASE E,, 
In this case the energy condition is: 
(o2/we2)Az+B,20. (22) 


The presence of A, in this expression is caused by 
the fact that, when the electric vector is in the 
plane of incidence, the Poynting vector, for a 
non-uniform wave, is not parallel to B (is not 
perpendicular to the planes of equal phase). In 
this case we have thus the possibility of B, being 
negative, provided (S.), 2 0. 

Taking the equality sign in (22), and elimi- 
nating A:z, Bz from (14), (15), (22), we obtain 
the pairs of values of a, and 6, for which the flow 
of energy across the boundary vanishes,’and for 
which we can consequently expect’ a transition 
from one solution for Az, B, to the other with 
reversed signs. We obtain: 

02 WE 


batt | ob, heawa'| 1+ in | (23) 


WE 02 €27w? 





If o2<e2w, Eq. (23), is represented by a line 
such as line 8 in Fig. 2, where lines a and y are 
the same as in Fig. 1. As before, the product 
A.B, is negative only right of the line y; thus, 
since the three hyperbolae meet at 0, the impor- 
tant part of curve @ is only that in the region of 
total reflection ; the remaining part corresponds 
to a spurious solution. A,, as before, vanishes 
along OP, and is negative right of BOP. B, 
vanishes along AO and now becomes negative in 
the region AOB. Along curve 8, where the 
average flow of energy across the boundary 
vanishes, we must pass from one solution to the 
other; there is thus a discontinuous transition 
from (substituting (22) in (15)): 


B,=([(o2/2wee)(2a.b,— woe) |}, 


24 
A,= aan [ (we2/2e2) (2a,b, — ywwo) |! ( ) 


to the solution with opposite signs. 


‘J. Picht, Ann. d. Physik 3, 433 (1929); F. Noether, 
Ann. d. Physik 11, 141 (1931). 
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b, 


ad 











a, 


Fic. 2. E parallel to the plane of incidence. Dis- 
continuity along BO. 


Particular cases, and the relative position of 
the three hyperbolae for different values of the 
parameters, are easily discussed. 


(6) USE OF THE REAL ANGLES OF INCIDENCE AS 
INDEPENDENT PARAMETERS 


An alternative formulation of the above treat- 
ment may be given assuming as independent 
variables, instead of a, and 6,, the real angles 
¢1, ¥1 that the normals to the planes of equal 
phase and of equal amplitude, respectively, make 
with the positive x-axis in the first medium. 
Formulae become then more expressive. Thus 
Eq. (21) representing the condition for the 
changing of sign of the product A,B, becomes 


cot¢: coty; = (o:/e2) —1, (25) 


and the condition for the polarizing angle, which 
makes H’=0 in the case E,, may be expressed as 


cot yg, coty2.+cot¢:2 coty, =2, (26) 


¢2, ¥2 being the angles corresponding to g; and 
¥1 in the second medium. 

An experimental arrangement for the veri- 
fication of the above theory would consist of a 
prism (large with respect to the wave-length) 
made of the material of medium 1, and bounded 
on the second face by medium 2. By varying 
the angle of incidence on the first face, and the 
angle of the prism, all values of g; and y; can be 
obtained. Such an experiment may provide a 
new method for the determination of con- 
ductivities. 
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The Relativistic Clock Problem 


E. L. Hitt 
Department of Physics, University of Minnesota, Minneapolis, Minnesota 


(Received April 18, 1947) 


The theory of uniformly accelerated motions based on the conformal group of transforma- 
tions in space-time is applied to the clock problem of relativity theory. Two solutions are 
found, both of which are at variance with the usual theory. The bearing of the problem on 
the relation between mechanics and electromagnetic theory is discussed briefly. 


1, INTRODUCTION 


NE of the most interesting problems fur- 
nishing a link between the special and the 
general theory of relativity is the so-called clock 
paradox. This was first introduced into relativity 
theory by Einstein! and has since been discussed 
many times in the literature; for present pur- 
poses we may take the easily accessible treatment 
of Tolman? as standard. The interest in the 
problem centers around the fact that it leads to 
Einstein’s important relation between gravita- 
tional potential and the rate of an ideal clock. 
In the problem as ordinarily stated,? two 
identical clocks are initially in coincidence and 
at rest in the laboratory system, and it is assumed 
that their readings have been adjusted and their 
rates synchronized. One of the clocks, B, is then 
given a strong acceleration for a short time, 
which brings it quickly up to a velocity u with 
respect to clock A, while the latter remains at rest 
in the laboratory system. Clock B continues to 
move with this velocity for a long time, after 
which it is given a further strong acceleration 
which reverses its motion and returns it towards 
A with the velocity —u. Just before reaching 
the latter, B is decelerated in such a manner as 
to bring it into coincidence and relative rest with 
respect to A. The readings of the two clocks are 
then compared. According to the usual theory, 
the elapsed time intervals, as measured by the 
two clocks, are related by the equation 


At, = Atg(1+4u?/2c’), (1) 


to terms in (u/c)*. The elapsed time, as read by 
the moving clock is thus /ess than that recorded 


1A, Einstein, Ann. d. Physik [4] 35, 898 (1911). 
2R. C. Tolman, Relatinty, Thermodynamics and Cos- 
mology (Oxford University Press, New York, 1934), p. 192. 
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by the stationary clock. The apparent paradox 
arising from taking the local-system of clock B 
as a standard of rest, in which A is then supposed 
to perform the contrary motion to that sketched 
above, is explained as an effect equivalent to 
that of the gravitational field which can be used 
to provide the apparent acceleration of clock 4 
in the local-system of B. 

A variant of the problem, in which the clocks 
are not required to be at relative rest when in 
coincidence, has been discussed by Méller! 
Apart from the clock problem itself, for which he 
appears to concur in the result of the usual 
theory, Mgller takes the occasion to discuss the 
problem of determining transformations between 
reference systems in uniform acceleration and, 
in fact, for even more general motions. We shall 
return to this aspect of Mgller’s discussion later 
in this paper. 

Despite the fact that in its kinematical aspects 
the problem appears to be entirely symmetrical 
between the two clocks, Einstein’s theory intro- 
duces a dissimilarity between them by the 
assumption that one of them is in an inertial 
system ; the local-system of the other clock then 
cannot be an inertial system. However, the per- 
sistent difficulty is that there is no @ priort way 
of making a decision as to which, if either, of the 
two clocks is in an inertial system. Its resolution 
in the relativistic argument appears to be quite 
as arbitrary as is the corresponding assumption 
in Newtonian mechanics. 

It has been shown in previous work‘ that one 
can introduce uniformly accelerated reference 
systems into relativity theory by employing the 
conformal group of transformations C, in space- 


3C. Moller, Det. Kgl. Danske Videnskabernes Selskab, 
Matematisk-F ysiske Meddelelser, Bind XX, Nr. 19 (1943). 
‘*E. L. Hill, Phys. Rev. 67, 358 (1945); 72, 143 (1947). 





AUGUST 1, 1047 











1, 1947 


aradox 
‘lock B 
Pposed 
etched 
lent to 
€ used 
lock A 


clocks 
hen in 
Igller.? 
lich he 

usual 
iss the 
tween 
1 and, 
e shall 
1 later 


spects 
trical 
intro- 
y the 
ertial 
: then 
e per- 
i way 
of the 
ution 
quite 
ption 


t one 
rence 
g the 
pace- 


Iskab, 
1943). 
947). 











CLOCK 


time, as a generalization of the Lorentz group 
L,, which establishes transformations between 
systems in uniform relative motion. It is the 
purpose in the present paper to discuss the 
application of this theory to the clock problem; 
in this way we shall be able to lift one of the 
deficiencies of the current theory by establishing 
explicitly the transformations between the local- 
systems of the two clocks and so can make the 
kinematical treatment symmetrical between 
them. However, as the theory is kinematical in 
character, and makes no direct appeal to 
mechanical principles, we can draw no particular 
conclusions concerning the causes of the apparent 
motions of the clocks, nor can we infer any rela- 


- tionship to any particular type of force field such 


as gravitation.° 

Our treatment will be limited to the case in 
which either of the clocks appears to perform a 
linear uniformly accelerated motion in the local- 
system of the other. If the clocks have two coin- 
cidences, they are then not at relative rest at 
either coincidence. An important point in the 
analysis must then be the method of establishing 
a comparison of the natural rates of the two 
clocks. We shall require that the rate at which 
either of the clocks appears to run in the local- 
system of the other, at the instant at which it 
comes to rest in that system, agrees with the 
rates of other neighboring clocks which are 
permanently at rest and synchronized in that 
system. The possibility of establishing extended 
covering times in the local-systems of both 
clocks, despite the inherent uncertainty as to 
whether either system is inertial in a mechanical 
sense, is founded on an appeal to the known form 
invariance of the electromagnetic field equations 
under C,; if such a covering time can be estab- 
lished in either system by purely optical means, 
it can be established in both. 


2. THE SUBGROUP OF THE CLOCK PROBLEM 


Our first task is to abstract from the full 
group, Cy, that subgroup which describes the 
type of motion appropriate to our version of the 


* Even should one bring in considerations concerning an 
external force field, it does not appear mandatory that it 
be interpreted as a gravitational field; the invariance of 
the Maxwell-Lorentz field equations suggests a formal 
interpretation in terms of electromagnetic fields. 
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problem. This can be built up from a translation, 
a uniform velocity, and a uniform acceleration 
along the x-axis. A study of the commutator 
tablet of C, shows that we must adjoin three 
further transformations, representing (a) a 
translation of the time origin, (b) a dilatation in 
space-time, and (c) a ‘“‘red-shift’’ transforma- 
tion. This gives us, in total, a 6-parameter 
subgroup of Cy. However, the two translations 
can be eliminated at once by taking as our funda- 
mental particle representing a clock, that one 
which corresponds to the origin of coordinates 
in space-time; the moving clock will then cross 
the origin of coordinates at the zero of local-time, 
moving along the x-axis, and the clocks will be at 
the origins of space-coordinates in their respec- 
tive local-systems. By including translations we 
should merely shift to other points representing 
the clocks. 

The general transformation of this type has the 
subsidiary differential equations’ 


dx /du= —tp\—}(r?+x* — y?—2*) pe—xrpstx, 
dy/du=y(—xp2—Tps+1), (2) 
dz/du=2(—xp2—rps+1), 

dr/du= —xpi—xtp2—3(7P?+2°+y? +2") pst r. 


For notation, we let (xo, Yo, Zo, and ro) be the 
space-time coordinates in the local-system Sz of 
clock B, while (x, y, z, and r) are those for the 
local-system S,4 of A. We have now to establish 
the relationship between these sets of coordinates 
by integration of Eqs. (2), with adjustment of 
the constants to suit the conditions of the pre- 
scribed motion. 

We start by considering only those points 
which move along the common x-axis of the two 
systems. On setting y=z=0 in Eqs. (2) we 
obtain 


d(r+x) /dp =(r+x)(1—p1) —$(7+x)?(p2+ ps), 
d(r—x)/dp = (r—x)(1+ 1) +3(7—x)*(p2— ps), 





SE. L. Hill, Phys. Rev. 68, 232L (1945). 
7For comparison with the notation of reference 4, we 
have the correspondence 


pi ars/C, pyran /c*, pye—ra/c, pons. 


The theory of the differential equations involved in the 
integration of a continuous group can be found in the 
books by J. E. Campbell, Theory of Continuous me 
(Oxford University Press, New York, 1903), p. 47; L. P. 
Eisenhart, Continuous Groups of Transformations (Prince- 
ton University Press, Princeton, New Jersey, 1933), 
Chapter 1. 
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and by integration we obtain the transformation equations for this class of points in the form 


(1 


y=2=0, 


— p1)(To+x0) 





tT+x= 


(1 — prem -9)# +3 (po+ ps)(To+%0)[1 —e-(-ee] 


(1+ p1)(to—xo) 





T—-X= 


To find the trajectory of the point representing 
the moving clock, B, in the local-system of 
clock A, we set x»=0 and eliminate ro. This 
yields the equation 


(x-+-w~1e* cosh pip)? 


—(r—w'e sinhp,u)? = (e*/w)?, 


where 
pP2—p pot+ps 
aory —e(itene)] + ——__ [1 —e(l-pvde], 
~ 2(1+ pr) 2(1—p1) 


In order to synchronize the clocks we now 
require that for x» =0 


(dro/dr)=1 at (x, r)=(x«, Tx), 


where x» and 7» are, respectively, the point and 
the instant at which the moving clock, B, comes 
to rest in the local-system of clock A. They are 
given by the relations 


Xy,=w'e*(1—coshpy), 7. =w~'e* sinhpy. 


On writing out this condition in full, we find 
the relation 


(p2— ps)(1—pi)[—o+oe tne 


+el(l-eve — el@tevde] 


+(p2+ps)(1+p1)[—o+oe%-9 


+elltevs — eh@-eve] =0 


with ¢= +1. 

Now when p; is considered to have an indeter- 
minate numerical value, the exponentials in this 
expression are linearly independent functions of 
u, so that this condition can be satisfied only by 
assuming 

(p2—ps)(1— pi) =0 
and 


(p2+ps)(1+ 1) =0. (3) 


(1 + pie tens a. (pe a ps)(to— Xo) [1 — en (tee) 





SS 


These relations yield just two solutions, which 
we designate as" 


Casea: pi:=+1, 
Case b: pi=-—l, 


pP2>= — pP3= Pa, 
p2= + p3= po. 


In the next two sections we shall examine the 
solutions for these cases separately. 


3. SOLUTION FOR CASE a. 
The differential Eqs. (2) reduce to the form 


d(r+x) = poly’ +2")dyu, 

d(r—x) =(r—x)[2+pa(r—x) ]dy, (4) 
dy =y[1+pa(r—x) ]dp, 
dz=2[ 1+ pa(r—x) du. 


By integration we find the transformation 
equations between the local-systems S4 and Sp, 
of the two clocks to be 


t+x= (to+X0) + rapa Vo? +20") sinhyg, 
T—X =)ee*(to—Xo), (5) 
y= AaVo, 2=),20 
with 
= 1/[e-* — pa(To—Xo) sinhy ]}. 


The trajectory of the moving clock, B, in S, 
is found by setting x9 = yo = 29 =0 and eliminating 
to. This gives us 


(x — pa! cothy)*— (7+ p.-')?=(1/pasinhy)*. (6) 


In the local-system of clock A, the moving 
clock, B, crosses the origin at time r=0, witha 
velocity £1=—tanhy, proceeds to the point 
X = pa! tanh(u/2) on the x-axis, which it reaches 
at time ts = —1/pa, coming to rest at this point. 
On reversing its motion it again reaches the 
origin, at which clock A is fixed, at time 2t, 
after which it continues its motion indefinitely 
along the x-axis, its speed approaching the value 
¢ asymptotically. By adjustment of the param- 
eters one can cause the ultimate motion to take 
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place in either direction. The whole motion is 
symmetrical between the two clocks, since the 
transformation of Eqs. (5) is inverted by inter- 
changing the coordinates and changing the sign 
of the group parameter uy. 

The relation between the local-times of the 
two clocks is obtained directly from Eqs. (5) 
on setting x0=Yo=20=0, which yields 


To= T+X. (7) 


We see from this that to=7 for x=0, so that 
the readings of the two clock agree at both coin- 
cidences, in contrast with the result expressed in 
Eq. (1). 

It follows also from Eq. (7) that dro/dr is 
greater than or less than unity, according as 
dx/dr is greater than, or less than, zero. This 
relation is peculiar in that it does not depend on 
the algebraic signs of p, and of yu. There appears 
to be an absolute distinction between left and 
right directions along the x-axis; this may also 
be interpreted as an absolute distinction between 
past and future which is conditioned by the 
imposed temporal ordering of the coincidences as 
first and second. 

It is of interest to note that we can arrive ata 
transformation corresponding to uniform relative 
motion of the two systems if we take the limiting 
case pa—0, with 8; = —tanhyu=const., with which 
we find 

x = (xo+BiT0)/(1+8:), 


tT =(to+Bix0)/(1+ 6), 

Cl -sal 

bf 1+ 8; Z0 
The transformation of Eqs. (8) is formed by a 
homogeneous Lorentz transformation with the 
velocity parameter, 81, followed by a dilatational 
transformation of all coordinates with a scale 
factor [(1—81)/(1+81)]*. The complete group 
of transformations of this type is obtained by 
adjoining dilatations to the inhomogeneous 
Lorentz group, giving rise to an 11-parameter 
subgroup of Cy. From the point of view of 
physical interpretation this result throws new 
light on the meaning of the usual simple Lorentz 
transformation, which provides no mechanism 


for the synchronization of clocks in uniform 
relative motion. 





CLOCK PROBLEM 





239 





We observe that our transformation does not 
satisfy the a@ priori requirement imposed by 
Mdller*® that it be of the form 


x=f(x0, To), S=Z0, t=h(xXo, To). 


¥=Yo, 


There appears to be no reason to believe that 
transformations of this functional form have any 
wider validity than the Lorentz group. 


4. SOLUTION FOR CASE b 


It will not be necessary to give the calcula- 
tional details for this case, since the work is 
similar to that of Section 3. The transformation 
equations are 


tT+x=)ee*(To+X0), 
T—XxX=(to—X0) —Aops(Yo' +50") sinhy, (9) 
Y=Awyo, 2=Arzo, 
with 
Ay = 1/[Le-*+ po(ro+x0) sinhy ]. 


This transformation can be obtained from that 
of Eqs. (5) by either of the formal substitutions 


(a) (x, y, 8, rT) —>(—x, ¥, 2,7), Pa>— ps, wn, 
(b) (x, ¥, &, T)—>(x, Tae, T), Pa + po, MoM. 


The remainder of the anaylsis proceeds exactly 
as in Section 3. The two clocks again appear 
to register the same time interval between coin- 
cidences. The moving clock appears to gain when 
it moves to the left and to lose when it moves 
to the right, which is just the reverse of the 
behavior found for case a. 


5. DISCUSSION 


The existence of two solutions of the clock 
problem shows that the conceptions of classical 
kinematical theory are not capable of charac- 
terizing the problem uniquely. Mathematically 
this arises from the circumstance that C, is a 
15-parameter group, while the corresponding 
group of classical kinematical theory has but 13 
parameters.‘ In order to see the nature of the 
influence of the enlarged character of C,, let us 
examine the behavior of the particles which move 
along the x-axis. We find the formula for the 
trajectory of such a particle from Eqs. (5) and 
(9) on setting yo=29=0 and eliminating ro. We 
obtain equations of the form 


(x —xo— p~' cothy)* — (r-Fxo+ p)? = (1/p sinhy)’. 












240 





For the apparent velocities of these particles we 
have 


dx /dt = +(rFxotp)/[(rFxotp)? 
+(p sinhy)~*}}. 


If we now consider those particles at great 
distances in either direction along the axis, for 
any fixed time 7, we find that |dx/dr|—+1 as 
|x9|—+-+ ©. In this way we are brought quite 
directly to a kinematical connection between the 
clock problem and the “expanding universe.” 
This is, in fact, simply another apsect of the 
relationship between C, and cosmological theory, 
which has been discussed by Robertson® and by 
Infeld and Schild.° 

To return to the comparison of the present 
analysis with the current relativistic theory, we 
observe that the divergence between the two 
mathematical procedures rests on the interpreta- 
tion assigned to the “‘line-element”’ associated 
with the group C,. The group is, in fact, charac- 
terized by a differential form of the type‘ 


A(x, y, 8, 7) (dr? —dx* —dy*® — dz’). (10) 


In the relativistic theory an arbitrary assignment 
is made of the “‘inertial’’ system in which \=1. 
In the present theory we refrain from making a 
unique assignment of line-element to the local- 


®*H. P. Robertson, Phys. Rev. 49, 755 (1936). 
*L. Infeld and A. Schild, Phys. Rev. 68, 250 (1945); 
ibid., 70, 410 (1946). 


E. L. 
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systems of the clocks; our discussion leaves the 
line-element of any particular coordinate system 
indefinite to the extent indicated by the form 
(10). We are concerned only with relations 
between coordinate systems, but not with the 
absolute specification of either system. From the 
mathematical point of view the present pro. 
cedure is more closely related to Weyl’s theory! 
than to that of Einstein. 

From the point of view of physical interpre. 
tation, this bifurcation is just that existing 
between mechanical and electrical theories, and 
it seems to the writer that by the apposition of 
the two procedures in the clock problem we are 
enabled to see the basic divergence between the 
two types of theory in a particularly elementary, 
but fundamental, light. The equations of me. 
chanics are form invariant under transforma- 
tions to systems moving with constant relative 
velocity, but the electromagnetic equations are 
insensitive to transformations involving uniform 
accelerations. The failure of the electromagnetic 
theory of mass seems to be a real measure of the 
incompatibility between the two theories. We 
seem here to be close to the root mathematical 
difficulties confronting the attempt to correlate 
quantum mechanics and the theory of relativity 
on any wider basis than that provided by the 
Lorentz group. 


1H. Weyl, Space-Time- Matter (Methuen Press, London, 
1922) Section 35. 
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Fine Structure of the Hydrogen Atom by a Microwave Method* ** 


Wuulis E. Lams, Jr. AND Ropert C. RETHERFORD 
Columbia Radiation Laboratory, Department of Physics, Columbia University, New York, New York 


(Received June 18, 1947) 


HE spectrum of the simplest atom, hydro- 
gen, has a fine structure! which according 
to the Dirac wave equation for an electron 
moving in a Coulomb field is due to the combined 
effects of relativistic variation of mass with 
velocity and spin-orbit coupling. It has been con- 
sidered one of the great triumphs of Dirac’s 
theory that it gave the “right’’ fine structure of 
the energy levels. However, the experimental 
attempts to obtain a really detailed confirmation 
through a study of the Balmer lines have been 
frustrated by the large Doppler effect of the lines 
in comparison to the small splitting of the lower 
or n=2 states. The various spectroscopic workers 
have alternated between finding confirmation? of 
the theory and discrepancies’ of as much as eight 
percent. More accurate information-would clearly 
provide a delicate test of the form of the correct 
relativistic wave equation, as well as information 
on the possibility of line shifts due to coupling of 
the atom with the radiation field and clues to the 
nature of any non-Coulombic interaction between 
the elementary particles: electron and proton. 
The calculated separation between the levels 
2?P, and 2?P3;2 is 0.365 cm and corresponds to a 
wave-length of 2.74 cm. The great wartime 
advances in microwave techniques in the vicinity 
of three centimeters wave-length make possible 
the use of new physical tools for a study of the 
n=2 fine structure states of the hydrogen atom. 
A little consideration shows that it would be 
exceedingly difficult to detect the direct absorp- 
tion of radiofrequency radiation by excited H 
atoms in a gas discharge because of their small 


* Publication assisted by the Ernest Kempton Adams 
ee for Physical Research of Columbia University, New 

ork. 

** Work supported by the Signal Corps under contract 
number W 36-039 sc-32003. 

! For a convenient account, see H. E. White, Introduction 
to Atomic Spectra (McGraw-Hill Book Company, New 
York, 1934), Chap. 8. 

*J. W. Drinkwater, O. Richardson, and W. E. Williams, 
Proc. Roy. Soc. 174, 164 (1940). 

*W. V. Houston, Phys. Rev. 51, 446 (1937); R. C. 
Williams, Phys. Rev. 54, 558 (1938); S. Pasternack, Phys. 
Rev. 54, 1113 (1938) has analyzed these results in terms 
of an upward shift of the S level by about 0.03 cm™. 
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population and the high background absorption 
due to electrons. Instead, we have found a 
method depending on a novel property of the 
27S, level. According to the Dirac theory, this 
state exactly coincides in energy with the 2°P, 
state which is the lower of the two P states. The S 
state in the absence of external electric fields is 
metastable. The radiative transition to the 
ground state 1°S, is forbidden by the selection 
rule AL= +1. Calculations of Breit and Teller‘ 
have shown that the most probable decay mecha- 
nism is double quantum emission with a lifetime 
of 1/7 second. This is to be contrasted with a 
lifetime of only 1.610~-* second for the non- 
metastable 2?P states. The metastability is very 
much reduced in the presence of external electric 
fields’ owing to Stark effect mixing of the S and P 
levels with resultant rapid decay of the combined 
state. If for any reason, the 27S, level does not 
exactly coincide with the 2?P, level, the vulnera- 
bility of the state to external fields will be re- 
duced. Such a removal of the accidental de- 
generacy may arise from any defect in the theory 
or may be brought about by the Zeeman splitting 
of the levels in an external magnetic field. 

In brief, the experimental arrangement used is 
the following: Molecular hydrogen is thermally 
dissociated in a tungsten oven, and a jet of atoms 
emerges from a slit to be cross-bombarded by an 
electron stream. About one part in a hundred 
million of the atoms is thereby excited to the 
metastable 27S, state. The metastable atoms 
(with a small recoil deflection) move on out of the 
bombardment region and are detected by the 
process of electron ejection from a metal target. 
The electron current is measured with an FP-54 
electrometer tube and a sensitive galvanometer. 

If the beam of metastable atoms is subjected to 
any perturbing fields which cause a transition to 
any of the 2?P states, the atoms will decay while 
moving through a very small distance. As a re- 
sult, the beam current will decrease, since the 


*H. A. Bethe in Handbuch der Physik, Vol. 24/1, §43. 
5G. Breit and E. Teller, Astrophys. J. 91, 215 (1940). 
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ical plot of galvanometer deflection due to 
interruption of the microwave radiation as a function of 
magnetic field. The magnetic field was calibrated with a 
flip coil and may be subject to some error which can be 
largely eliminated in a more refined apparatus. The width 
of the curves is probably due to the following causes: 
(1) the radiative line width of about 100 Mc/sec. of the ?P 
states, (2) hy so ae ae of the *S state which amounts 
to about 88 (3) the use of an excessive intensity 
of radiation which's ives increased a in the wings 
of the lines, and (4) inhomogeneity of the magnetic field. 
No transitions from the state 2*S;(m=-—4) have been 
observed, but atoms in this state may be quenched by 
stray electric fields because ofthe more nearly exact degen- 
eracy with the Zeeman pattern of the ?P states. 


Fic. 1. A ty 


detector does not respond to atoms in the ground 
state. Such a transition may be induced by the 
application to the beam of a static electric field 
somewhere between source and detector. Transi- 
tions may also be induced by radiofrequency 
radiation for which hy corresponds to the energy 
difference between one of the Zeeman components 
of 27S, and any component of either 2?P, or 2?P 3/2. 
Such measurements provide a precise method for 
the location of the 22S, state relative to the P 
states, as well as the distance between the latter 
states. 

We have observed an electrometer current of 
the order of 10-'4 ampere which must be ascribed 
to metastable hydrogen atoms. The strong 
quenching effect of static electric fields has been 
observed, and the voltage gradient necessary for 
this has a reasonable dependence on magnetic 
field strength. 

We have also observed the decrease in the 
beam of metastable atoms caused by microwaves 
in the wave-length range 2.4 to 18.5 cm in various 
magnetic fields. In the measurements, the fre- 
quency of the r-f is fixed, and the change in the 
galvanometer current due to interruption of the 
r-f is determined as a function of magnetic field 
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Fic. 2. Experimental values for resonance magnetic 
fields for various frequencies are shown by circles, The 
solid curves show three of the theoretically — 
variations, and the broken curves are obtained by shifting 
these down by 1000 Mc/sec. This is done merely for the 
sake of comparison, and it is not implied that this would 
represent a “best fit.” The plot covers only a small range 
of the frequency and magnetic field scale covered by our 
data, but a complete plot would not show up clearly ona 
small scale, and the shift indicated by the remainder of 
the data is quite compatible with a shift of 1000 Mc. 


strength. A typical curve of quenching versus 
magnetic field is shown in Fig. 1. We have 
plotted in Fig. 2 the resonance magnetic fields for 
various frequencies in the vicinity of 10,000 
Mc/sec. The theoretically calculated curves for 
the Zeeman effect are drawn as solid curves, 
while for comparison with the observed points, 
the calculated curves have been shifted down- 
ward by 1000 Mc/sec. (broken curves). The 
results indicate clearly that, contrary to theory 
but in essential agreement with Pasternack’s 
hypothesis,* the 22S, state is higher than the 
2?P, by about 1000 Mc/sec. (0.033 cm- or about 
9 percent of the spin relativity doublet separation. 
The lower frequency transitions *S,(m=})7? 
2P,(m= +4) have also been observed and agree 
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well with such a shift of the 2S, level. With the 
present precision, we have not yet detected any 
discrepancy between the Dirac theory and the 
doublet separation of the P levels. (According to 
most of the imaginable theoretical explanations 
of the shift, the doublet separation would not be 
affected as much as the relative location of the S 
and P states.) With proposed refinements in 
sensitivity, magnetic field homogeneity, and 
calibration, it is hoped to locate the S level with 
respect to each P level to an accuracy of at least 
ten Mc/sec. By addition of these frequencies and 
assumption of the theoretical formula Av= 7ga*R 
for the doublet separation, it should be possible to 
measure the square of the fine structure constant 
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times the Rydberg frequency to an accuracy of 
0.1 percent. 

By aslight extension of the method, it is hoped 
to determine the hyperfine structure of the 2°S, 
state. All of these measurements will be repeated 
for deuterium and other hydrogen-like atoms. 

A paper giving a fuller account of the experi- 
mental and theoretical details of the method is 
being prepared, and this will contain later and 
more accurate data. 

The experiments described here were discussed 
at the Conference on the Foundations of Quantum 
Mechanics held at Shelter Island on June 1-3, 
1947 which was sponsored by the National 
Academy of Sciences. 
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OME considerations relating the spectroscopic 
and photochemical properties of polyatomic 
gases with their use in Geiger-Mueller counters 
have been developed by the writer.! The self- 
quenching counter is able to operate without the 
usual large resistance because of the absence of 
secondary emission of electrons when the positive 
ion sheath reaches the cathode and is neutralized. 
The neutralization of the ions (at about 10-7 cm 
from the wall) results in neutral atoms or mole- 
cules having an excitation energy of (J—¢) 
approximately, where J is the ionization energy 
and gy the photoelectric work function of the 
surface. In the case of a monatomic or diatomic 
(non-quenching) gas this excitation energy is 
transferred to the metal in about 10-” sec. and 
secondary electrons are liberated if J>2¢. In the 
case of most polyatomic gases (with four or more 
atoms), the neutralized molecules will predis- 

* Also consultant at Clinton Laboratories, Oak Ridge, 
Tennessee. 

1S. A. Korff and R. D. Present, Phys. Rev. 65, 274 
(1944). The theoretical considerations in this paper were 


mainly developed by R Present; the experimental 
bai and its interpretation were done mainly by S. A. 


(Received May 23, 1947) 





sociate (in about 10-" sec.) before they can make 
an inelastic collision with the wall, and the 
emission of secondary electrons will be negligible. 
In order to determine whether a given gas will 
show self-quenching action, it is sufficient to 
examine its absorption spectrum at wave-lengths 
corresponding to an excitation energy of (J—). 
Continuous absorption, or a line broadening 
corresponding to a predissociation lifetime of 
10-* sec., indicate that the molecule undergoes 
primary photo-decomposition. Such a gas will 
make a self-quenching counter. If a mixture of 
quenching and non-quenching gases is used, the 
counter will show self-quenching action only if 
the ionization potential of the non-quenching gas 
exceeds that of the quenching gas. Under these 
circumstances electron transfer changes the com- 
position of the positive ion sheath so that, when 
it reaches the cathode, there is a negligible 
number of non-quenching ions. A more detailed 
discussion will be found in reference 1. 

We should like to point out here that the 
mechanism of secondary-emission quenching is 
compatible, under special circumstances, with the 
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quenching gas being diatomic. The halogens Cls, 
Brz, and I; are examples of diatomic gases which 
would be expected to show self-quenching action.” 
The absorption spectra of these gases have been 
investigated into the far ultraviolet. The ioniza- 
tion potentials are respectively 13.2, 12.8, and 9.7 
volts. If the cathode surface is carbon’ (photo- 
electric work function ~5 volts) the values of 
(J—¢) correspond to wave-lengths of 1500, 1600, 
and 2600A for the three gases, respectively. In 
this region of the spectrum both Cl; and Br. 
show continuous absorption corresponding to 
primary photo-decomposition.‘ One should ex- 
pect that these gases would give rise to no 
secondary emission. The ultraviolet absorption 
spectrum of I, shows several bands in the region 
from 2760 to 1750A.° Several of these bands are 
sharp and are accompanied by fluorescence, indi- 
cating the absence of predissociation except as 
induced by collisions. Decomposition of the 
neutralized I; molecules at the cathode surface is 
not indicated as in the case of the other halogens. 
Nevertheless, if the cathode surface be of carbon,* 
secondary emission cannot occur because the 
ionization potential is less than twice the work 
function. The work function of carbon is 4.7 
volts for ordinary contamination; however, 
halogen contamination generally raises the work 
function by several tenths of a volt.’ For all 
surfaces for which J <2¢ the secondary emission 
by positive ions is quenched. 

2S. H. Liebson, Bull. Am. Phys. Soc. 22, 2 (1947), has 
recently reported that counters filled with suitable mix- 


tures of halogen and noble gases show self-quenching (fast) 
counter action. 

3 The inner surfaces of the cathodes used in reference 2 
were coated with mpetes (graphite). 

‘ Both the Cl, and Bre continuous absorption r 
the ultraviolet extend down to 1560A (limit of 
tions) and the maximum of the Cl: absorption ap 
be near this value. H. Cordes and H. Sponer, 
Physik 63, 334 (1930). 

50. Olden , Zeits. f. Physik 25, 136 (1924); P. 
Pringsheim and B. Rosen, Zeits. f. Physik 50, 1 (1928); 
D. T. W. , Phys. Rev. 47, 1 (1935). 

6 Predissociation induced by gas collisions would be 
ineffective in preventing secondary emission since the 
neutralized molecule would make an inelastic collision 
with the wall before it could collide with another molecule. 
The van der Waals field at the surface has a sufficient 
range to induce predissociation of the neutralized molecule 
but the interaction is insufficient to reduce the lifetime to 
less than 10-™ sec. 

7C. Ouellet and E. K. Rideal, J. Chem. Phys. 3, 150 
(1935). 
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The possible effects of electron attachment 


cannot be neglected when the halogens are used 


in appreciable concentrations. After the initial 
ionizing event there will evidently be a smalj 
number of negative ions moving in toward the 
center wire from the outer regions of the counter, 
When a negative ion reaches the high field region 
near the wire, an electron will detach itself and 
initiate a new avalanche unless the field is below 
the counting threshold. As the positive ion sheath 
moves outward, in a fast counter the field near 
the wire increases; the threshold field is reached 
when the positive ions have traveled the critical 
distance (about half the distance to the wall), 
The negative ions in a fast halogen counter must 
all have been neutralized before the positive ion 
sheath reaches the critical distance. The reason 
for this is found partly in the greater mobility of 
the negative compared to the positive ions! 
Furthermore, the negative ions may be neutral- 
ized in crossing the positive ion sheath (the 
electron affinities of the halogen atoms are less 
than the ionization potentials of the molecules), 
If the partial pressure of halogen gas is not too 
large, it appears likely that the performance of 
these counters should be unimpaired by multiple 
pulses and spurious counts arising from negative 
ions. However, the concentration of halogen 
must be large enough to quench metastable rare- 
gas atoms.? Finally, just as in the case of oxygen 
counters,’ attachment of the electrons formed in 
the initial ionizing event can cause a mean time 
lag of as much as 10~ sec. between the passage of 
the ionizing particle and the initiation of the 
discharge (unless the concentration of halogen is 
only a few percent). This would make the halogen 
counter unsuitable for coincidence experiments in 
which the efficiency must be high and _ the 
resolving time low. 


®H. G. Stever, Phys. Rev. 61, 38 (1942). 

* Immediately after electron attachment the halogen 
molecule ion dissociates into an atomic ion and a neutral 
atom. Because of its smaller mass and effective cross section 
the negative atomic ion should have a greater mobility 
than the positive molecular ion. The discharge terminates 
before the mobilities can be affected by impurities. 

10M. E. Rose and W. E. Ramsey, Phys. Rev. 59, 616 


(1941). 
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A Low Temperature Transformation in Lithium* 


CHaRLeEs S. BARRETT 
Institute for the Study of Metals, University of Chicago, Chicago, Illinois 
June 12, 1947 


ITHIUM has been reported as having a body-centered 

cubic structure in the range from ordinary tem- 
peratures to liquid air temperature."* X-ray diffraction at 
low temperatures now reveals, however, that a transforma- 
tion can be induced in the metal if it is plastically deformed 
at temperatures in the vicinity of — 196°C. 

Analysis of the Geiger-counter spectrometer record of 
the diffraction pattern shows that the new phase produced 
by this treatment is face-centered cubic in structure with 
the lattice constant a9=4.41A at —196°C. The body- 
centered cubic phase in the same sample at this tempera- 
ture has dp = 3.50A. There is an expansion of atomic radius 
accompanying the increase of coordination number from 
8 to 12, as would be expected, and because this expansion 
has the value 2.8 percent it follows that the calculated 
densities of the two phases are identical within experi- 
mental error. The calculated density is 0.534. 

The transformation to the face-centered cubic form is 
accompanied by a series of audible clicks, as in the twinning 
of tin or magnesium and the formation of martensite. By 
analogy with these processes it may be concluded that the 
transformation goes by abrupt shear movement in small 
isolated regions. The constraints imposed by the material 
around these transforming regions can account for the 
fact that we have been unable to transform much more 
than half of any of the polycrystalline samples by any 
type of straining we have tried. Presumably then, proper 
straining of a single crystal might lead to complete trans- 
formation. The face-centered phase disappears fairly 
rapidly when heated above about — 117°C. 

We were led to search for low temperature transforma- 
tions in various metals by C. Zener’s theory of the sus- 
ceptibility of body-centered cubic phases to instability at 
low temperatures,* * which may be summarized as follows. 
A homogeneous shear of 0.35 in a body-centered cubic 
structure along the (110) plane and the [110] direction 
will produce an atomic arrangement that is very nearly 
face-centered cubic. Consider the relative free energy of 
the two phases. The free energy of the body-centered cubic 
phase increases as the temperature is lowered, and increases 
more rapidly than would the free energy of a face-centered 
cubic phase, provided the thermal vibrations of the atoms 
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in some direction of the body-centered structure are ab- 
normally large. In metals and alloys of body-centered 
cubic structure that have filled inner shells of electrons, 
the shear constant (Ci:—Ci2)/2 should be small and is, 
in fact, very small; as this constant applies to shear in 
the (110) plane in the [110] direction, the thermal vibra- 
tion amplitudes in this direction should be large, and since 
the entropy is proportional to the logarithm of the ampli- 
tude it follows then that the entropy should be large. The 
free energy of a phase with large entropy rapidly increases 
as the temperature is lowered, since F=U—TS and 
0F/8T =—S, where F is the free energy, U the internal 
energy, J the temperature, and S the entropy; it is not 
unlikely that F for the body-centered phase in such cases 
will exceed that for the face-centered phase and make a 
transformation possible, as in some compositions of beta- 
brass, where (Ci,—(C2)/2 is twice as large in the face- 
centered form as in the body-centered. Such transforma- 
tions should be aided by shear in the direction for which 
the shear modulus is low, both because it reduces the 
sluggishness of the transformation by aiding the atoms in 
surmounting the potential barrier in this direction and 
because this is the direction required for the shear move- 
ment of the transformation. 

A more detailed account of the investigation of this 
metal and others will be published elsewhere. 

* This work was sup by the Office of Research and Inventions, 
U 5. ye (Contract No. N-6ori-20-IV). 

F. Simon and E. Vohsen, Zeits. f. physik. Chemie 133, 165-187 

OT ca J. Phys. Chem. 32, 354-359 (1928). 

*C. Zener, Phys, Rev. 71, 846 (1947). 


4 C. Zener, Elasticity and Anelasticity of Metals (University of Chicago 
Press, Chicago, to be published). 





On the Hyperfine Structure of the Ground 
State of H and D Atoms 


Otto HALPERN 
Columbia University, New York, New York 
June 16, 1947 


ECENTLY'! measurements of the hyperfine structure 

of H and D atoms have been carried out which are 
extremely interesting from a theoretical as well as an 
experimental point of view. 

According to the latest precision measurements, not 
only the absolute values of the hyperfine structure splittings 
vH, vp differ from the theoretical value by about one part 
in 400, but also the ratio of these quantities differs from 
the theoretical value by about 1 part in 1700, the calculated 
value always being smaller than the observed. The theo- 
retical evaluation was carried out with the aid of the 
well-known Fermi formula 


2I+1 
= I 





a unuo|¥(0)|*, (1) 
in which the nuclear spin is denoted by J; uy is the mag- 
netic moment of the nucleus, yo is Bohr’s magneton, and 
¥(0) the Schroedinger wave function at the origin. Through 
|y(0)|?, the Bohr radius of the ground state enters into 
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the expression of the hyperfine structure splitting, |¥(0)|? 
being proportional to a~*. It is customary to eliminate the 
Bohr radius by introducing the Rydberg constant for 
infinite mass, together with the reduced mass of the 
electron. It then turns out that the hyperfine structure 
splitting is given by the expression 


3 
v= (= ‘ (2) 
I m 


Here C contains universal constants and numbers, while 
m, stands for the reduced mass of H or D, respectively. 
me denotes the electronic mass. The ratio vy/vp is then 
given by 





va/vp = (4/3)(un/up)(mu/mp), 


and the numerical value for vy/»p as observed, is 4.3416 
compared with a computed value of 4.3393. 

Equation (1) is derived by neglecting the small com- 
ponents of the Dirac equation and replacing the large 
components by Schroedinger functions. 

The discrepancy observed which, if the accuracy of 
observation is sufficient, in the case of the ratio at least, 
cannot be ascribed to inaccurate values of the universal 
constants, makes it advisable to re-examine the derivation 
of (1). 

We have obtained a value for » by consistently using 
Dirac’s equation, retaining all four components through 
the perturbation calculation and using the rigorous Dirac 
functions in the evaluation of the matrix element for the 
perturbed energy. The hyperfine structure splitting apart 
from numerical factors is now given by 





= Pe a" —(E/me?)*}'=C' €un (2R/me?)'. (3) 


(h/mc)? (h/mc)? 
The second equation in (3) is obtained by the somewhat 
arbitrary insertion of the empirical Rydberg constant in 
place of R., which would, of course, follow from the Dirac 
equation. 

This calculation leads, within the accuracy aimed at, 
to the same value as given by (2), with the one difference 
that the ratio of reduced to electronic mass appears in 
the three-halves rather than in the third power. This 
correction diminishes the discrepancies between the ob- 
served and calculated values of vq and yp, as follows. For 
H, the discrepancy is reduced to one part in 600; for D, 
to one part in 500; both deviations are obviously still 
large if one believes in the presently accepted values of the 
universal constants. 

The ratio »y/»p, on the other hand, now differs from 
its calculated value by only one part in 8000; this is much 
smaller than the accuracy claimed for the earlier determi- 
nations of wH/up, which enters as a factor into (3) and is 
assumed to be known to about 1 part in 3000. 

In interpreting this result, several points must be kept 
in mind: The accuracy of the experimental determination 
of wH/up and the calculated value which contains it as its 
most uncertain element, are not yet sufficiently good 
to exclude a different dependence on the ratio of the 
reduced masses. 


THE EDITOR 


The theory here used is obviously not consistent: we 
have carried out all calculations with the one-body Dirac 
equation and taken into account the two-body nature of 
the problem by the empirical introduction of Ry and Rp 
This point will need further theoretical study. 

We have also investigated the question of how the 
electronic magnetic moment may be expected to depend 
on the nuclear mass and have found different results 
depending on the physical interpretation given to the 
coordinates which enter into the Dirac equation. 

A detailed paper will follow shortly. 

1 J. E. Nafe, E, B. Nelson, and I. I. Rabi, Phys. Rev. 71, 914 (1947), 


I am greatly indebted to the authors for telli ; 
before publication. ng me about their results 





Nuclear Capture of Mesons and the 
Meson Decay 


B. PoNTECORVO 


National Research Council, Chalk River Laboratory, Chalk River 
Ontario, Canada : 


June 21, 1947 


HE experiment of Conversi, Pancini, and Piccionit 
indicates that the probability of capture of a meson 
by nuclei is much smaller than would be expected on the 
basis of the Yukawa theory.** Gamow‘ has suggested that 
the nuclear forces are due exclusively to the exchange of 
neutral mesons, the processes involving charged mesons 
and the 8-processes having probabilities which are smaller 
by a factor of about 10”. 

We notice that the probability (~10° sec.—') of capture 
of a bound negative meson is of the order of the probability 
of ordinary K-capture processes, when allowance is made 
for the difference in the disintegration energy and the differ- 
ence in the volumes of the K-shell and of the meson orbit. 
We assume that this is significant and wish to discuss the 
possibility of a fundamental analogy between 8-processes 
and processes of emission or absorption of charged mesons. 

An immediate consequence of the experiments of the 
Rome group’ is that the usual interpretation of the 8-process 
as a “two-step” process (“‘probable”’ production of virtual 
meson and subsequent §-decay of the meson) completely 
loses its validity, since it would predict too long 8-lifetimes: 
the meson is no longer the particle responsible for nuclear 
8-processes, which are to be described according to the 
original Fermi picture (without mesons). Consequently 
there is no need to assume that charged mesons have 
integral spin, as the Yukawa explanation of 8-processes 
required. Once we believe that the ordinary 8-process is not 
connected in any way with the meson, it is difficult to see 
strong reasons for the usual assumption that the meson 
decays with emission of a §-particle and a neutrino. We 
shall consider then the hypothesis that the meson has spin 
$h and that its instability is not a 8-process, in the sense 
that it does not involve the emission of one neutrino. The 
meson decay must then be described in a different way: it 
might consist of the emission of an electron and a photon or 
of an electron and 2 neutrinos® or some other process. 
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In the hypothesis that the meson decay is not a 8-process 
(meson of spin }) the process of nuclear absorption or pro- 
duction of a single meson would be accompanied by the 
emission of a neutrino. This analogy between 8-particles and 
mesons suggests, in addition, that just as the production of 
single g-particles is extremely unlikely, while the produc- 
tion of electron pairs is a very likely phenomenon, so the 

uction of a single charged meson would be very 
unlikely, while the production of pairs of mesons would be 
quite probable. The experimental evidence is, in fact,* that 
most, if not all, of the meson showers are created in con- 
nection with large Auger showers. 

The assumption that the emission or absorption of one 
meson is accompanied by the emission of a neutrino would 
explain in a natural way a somewhat puzzling experimental 
result. Among the few pictures of a meson stopping in the 
gas of a cloud chamber, no “‘star” has been observed at the 
end of the meson track.’ The absence of a star must be due 
to a process leaving the capturing nucleus in a not too 
excited state: the mechanism proposed here would explain 
that the capture of a negative meson from a nucleus Z 
results in a nucleus Z—1 close to its ground level, since the 
excess energy could be carried away by the neutrino. 
Actually, in such a process we should expect that the 
emission of a ‘neutrino of high energy with consequent 
production of the nucleus Z—1 in a state of low excitation 
would be more likely than the emission of a neutrino of low 
energy with the production of the nucleus Z—1 in a state of 
high excitation (cf. K-capture process). 

The hypothesis that the meson decay is not a 8-process, 
while the meson absorption is a 8-process, does not require 
that hypothetical particles such as neutral mesons are 
invoked to account for nuclear forces. In fact, a heavy 
electron pair theory of nuclear forces was successfully de- 
veloped by Marshak.* Moreover, a pair theory is capable of 
accounting, at least in principle, for the existence of 
processes in which several pairs of mesons are produced in a 
single act, as suggested by Heisenberg in connection with a 
different problem.® 

Returning to the actual decay of the meson, an experi- 
ment suggests itself which might answer the following 
question : Is the electron emitted by the meson with a mean 
life of about 2.2 microseconds accompanied by a photon of 
about 50 Mev? This experiment is being attempted at the 
present time, since it is felt that the available analysis" of 
the soft component in equilibrium with its primary meson 
component is probably insufficient to decide definitely 
whether the meson decays into either an electron plus 
neutral particle(s) or electron plus photon. 

Pe. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 71, 209 tly 
TE. Fermi, B. Teller, and VWelsskopt, Phys. Rev. 91, 314 (1947). 
*J. A. Wheeler, Phys. Rev. 71, 320 (1947). 

‘*G. Gamow, Phys. Rev. 71, 550 (1947). See also G. Gamow and E. 

Teller, Phys. Rev. 51, 289 (1937). 

*W. Nordheim, Phys. Rev. 59, 544 (1941). 

one A. Loveredo, and V. Tongiorgi, Phys. Rev. 70, 852 

? See for a critical survey: T. H. Johnson and R. P. Shutt, Phys. Rev. 

61, 380 (1942). 

*R. E. Marshak, Phys. Rev. 57, 1101 (1940). 
* References can be found in Cosmic Radiation, edited by W. Heisen- 


berg (Dover Publications, New York, 1946), p. 127 
¥ See reference 9, pp. 84-97. 
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An On-Orbit Injector for Betatrons 
and Synchrotrons 


Davip C, DEPACKH AND ANDREW V. HAEFF 
Naval Research Laboratory, Washington, D. C. 
May 29, 1947 


NE of the primary limitations on the yield from 

present electron accelerators of the betatron and 
synchrotron type lies in the injection process. The injector 
must lie off the orbit in which the acceleration is supposed 
to occur in order not to interfere with the acceleration 
process, and this necessarily implies radial oscillations of 
the injected electrons of an amplitude at least equal to 
the distance between the stable orbit and the injector. 
The frequency of these radial oscillations is approximately 
vo(1—m)*, where vo is the rotation frequency of the elec- 
trons, and n is the exponent in the expression H = H,(r/re)* 
law of decrease of the magnetic field with radius. If » has 
a nominal value of 0.5, it is seen that the oscillation fre- 
quency is 0.71 of the electron frequency and that as a 
consequence we cannot expect any appreciable damping 
of this oscillation by increasing magnetic field in a time 
comparable to several cycles of this oscillation. Moreover, 
the beneficial effect of such damping is further reduced if, 
as is the practice, the injection is performed at high voltage, 
since in a given time interval the relative change of field 
is less with a large field than a small one for a field in- 
creasing approximately linearly with the time. The net 
result of these conditions is that although injected electrons 
may clear the injector on the first transit because of the 
difference between vo and vo(1—m), nevertheless a great 
many of them will strike it on subsequent transits and 
thus be lost. 

In order to circumvent this difficulty, it is proposed to 
employ an injector which surrounds the stable orbit so 
that electrons rotating in the stable orbit may pass through 
the injector structure without interference. Thus in sub- 
stance the injector (see Fig. 1) consists of an annular 





Fic. 1. Schematic drawing of on-orbit injector. The magnetic field is 
normal to the paper. 


cathode and a series of electron lenses which narrow and 
collimate the electron beam, so that as it emerges from the 
injector it appears ideally as a cylindrical tube the axis of 
which is the stable orbit. One can design the lens system 
to make the diameter of this tube of such a size that 
oscillation within the beam is of no importance. The part 
of the injector behind the cathode can be designed so 
that electrons, passing through the injector on subsequent 
transits and having been defocused because of space charge 
or for some other reason, can be refocused as they pass 
through the low potential region near the center of the 
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cathode ring. Since the electron velocity during the in- 
jection pulse falls as the electron enters the injector, goes 
through a minimum as it passes through the cathode ring, 
and rises to its initial value as the electron leaves the gun, 
it is necessary to supply a compensating magnetic field in 
the neighborhood of the injector during the injection pulse. 

An injector of this type has been designed at the Naval 
Research Laboratory and is under development. 





The Beta-Ray Spectrum of C™' 


Pau. W. Levy 
Clinton Laboratories, Oak Ridge, Tennessee 
June 9, 1947 


180° beta-ray spectrometer of 10-cm radius, using a 

Geiger-Mueller tube as the detecting device, has 
been used to investigate the continuous beta-ray spectrum 
of the 4700-year C™ activity. It was necessary to use a 
very thick source because of the long half-life and low 
concentration of C“ in the available material. The spectrum 
obtained with this source appears simple, that is, it has 
only one component and has a maximum energy of 0.154 
Mev. As is to be expected from the fact that no gamma-rays 
have been detected accompanying the decay of C™, no 
conversion electron lines appear. 

The maximum energy of the beta-ray spectrum was 
determined from the Kurie plot, Fig. 1, which was con- 
structed from that part of the spectrum high enough in 
energy so that absorption in the window of the detecting 
G. M. tube is negligible. One cannot say, however, that 
this part of the spectrum is entirely free from the effects of 
self-absorption and back scattering due to the thickness 
of the source. These effects might possibly combine in a 
way that introduces an undetermined error in the 0.154- 
Mev end point and, therefore, one cannot claim the 2 
percent inherent accuracy of the spectrometer at this 
energy. Yet it is surprising that the Kurie plot is straight, 
and, if it were not for the unusual thickness of the source, 
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Fic. 1. Kurie plot of C™ beta-ray spectrum, 
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one would interpret this as evidence that the spectrum was 
free from large disturbing effects. 

The author wishes to express his thanks to L. D. Norris 
for preparing the C™ source. 

This work was done under the auspices of Manhattan 
District Contract W 7405-eng. 39. 

' . oe ; : 
tained in Plutonium Project Report CP-3702 winch is tobe patie 


a contribution to the Plutonium Project Record volume on nuclear 
physics. 





On the Change of the Mean Life of 
Negative Mesotrons 


HAROLD K, TICHO AND MARCEL SCHEIN 
Depariment of Physics, University of Chicago, Chicago, Illinois 
June 23, 1947 


ECENT experiments"? have demonstrated that nega- 
tive mesotrons, stopping in solid materials of low 
atomic number, do disintegrate; this result is in contra- 
diction with the calculations of Tomonaga and Araki* who 
predicted that mesotrons carrying a negative charge 
should be captured by the nuclei of solid absorbers, what- 
ever their atomic number. Subsequently, Fermi, Teller, 
and Weisskopf,‘ and also Wheeler® investigated the energy 
loss of mesotrons after they have stopped ionizing and 
concluded that negative mesotrons should fall into the 
(mesotron) K-shell within 10-" second. From the K-orbit, 
mesotrons can either be captured by the atomic nucleus or 
disintegrate spontaneously. Their decay constant should 
therefore be A+A, where 1/A=7o is the mean life for dis- 
integration and A the probability per second for capture. 
The fact that the disintegration of negative mesotrons is 
observed in light absorbers indicates that A is not much 
larger than A; hence, A~5X 10° sec.~', which is in strong 
disagreement with meson theories which yield A~10" 
sec.—!. According to Wheeler® A should increase roughly 
as 2‘. 

If it is true that the capture probability is considerably 
smaller than that calculated on the basis of meson theories 
of nuclear forces, then the disintegration electrons arising 
from the decay of negative mesotrons should also follow a 
disintegration curve with the modified decay constant 
A+A. The number of negative mesotrons dN, disin- 
tegrating in the time interval dt, is equal to the original 
number of negative mesotrons No~ times the probability 
that they have neither disintegrated nor been captured 
before t, exp(——At), times Adt. The integration of this 
expression yields 


Na = No [A/(A+A)] exp[— (A+A)é]. 


The capture constant A could, therefore, be measured by 
determining the mean life 1/(A+A) of negative mesotrons 
alone. In such a determination one should use ro=2.15 
 10-* second for 1/A. This value of ro was obtained by 
Nereson and Rossi® chiefly in lead, where negative meso- 
trons do not disintegrate. 

If the experiment is carried out in the usual way, i.e., 
without differentiating between positive and negative 
mesotrons, a small change of the disintegration constant 
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Fic. 1. The integral-disintegration curves which should be obtained if 

tive mesotrons both disintegrate and get captured. The values 

(1/A) are the “‘mean-life” value calculated from the average slope of 
these curves. 


should still be observable in materials of suitable atomic 
number. This change of the disintegration constant is due 
to the fact that the experimental data obtained really 
represent a superposition of two decay curves, one pro- 
duced by positive mesotrons with a mean life 1/A and one 
produced by negative mesotrons with a mean life 1/(A+-A). 
Hence, the total disintegration curve is represented by 
the following expression : 


Na= Not exp(—)+ No“ [A/(A+A)] exp[— (A+A)é]. 


Disintegration curves of this kind for values of A between 0 
and «©, and assuming an equal number of positive and 
negative mesotrons, are shown in Fig. 1. The quantities 
(1/A)’ are the “changed” values of the mean life obtained 
by determining the average slope of the curves for a time 
interval from 1X10~* to 4.5X10°-% sec. This “change” 
may be as large as 11 percent. From the curves, however, 
it is clear that in such an investigation particular attention 
must be paid to very short decay times, especially if the 
capture constant is large. This consideration applies also 
to experiments where the ratio of disintegration electrons 
to stopped mesotrons is measured. 

1M. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 71, 209 (1947). 

? T. Sigurgeirsson and A. Yamakawa, Phys. Rev. 71, 319 (1947). 

ic —- = G. Araki, Phys. Rev. 58, 90 (1940). ' 

. Fermi, E. Teller, and V. Weisskopf, Phys. Rev. 71, 314 (1947). 


* John A. Wheeler, Phys. Rev. 71, 320 (1947). 
*N. Nereson and B. Rossi, Phys. Rev. 64, 199 (1943). 





Analysis of the Hyperfine Structure in the Micro- 
wave Spectrum of the Symmetric 
Top Molecule CH;I* 


WALTER Gorpy, A. G. SMITH, AND JAMES W. Simmons 
Department of Physics, Duke University, Durham, North Carolina 
June 7, 1947 


HE preliminary measurements! of the J=1 to J=2 
rotational transition of CH;I have been repeated at 
a lower temperature, —70°C. As a result of the increased 
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absorption at the lower temperature, additional lines have 
been observed and a more accurate estimate of relative 
intensities has been made.' The measurements at low 
temperature revealed that none of the observed lines 
arises from molecules in excited vibrational states. 

A complete analysis of the hyperfine structure resulting 
from the nuclear quadrupole moment of iodine has now 
been made. In Fig. 1 the theoretically predicted structure 
is compared with that experimentally observed.* The 
particular formula used to calculate the positions of the 
CHsglI lines is that used by Coles and Good? to calculate 
the quadrupole effects in NHs: 


Tr = [eQ0* V /dz*}-[(3K?/2J(J+1)—4) 
-[3C(C+1)—4J(J +110 +1)/ 
(2J —1)(2J +3)(2I—1)(27+3)], (1) 
where 


Try represents the hyperfine splitting of the rotational 
energy levels, 

C= F(F+1)—I(I+1)—-J(J +1), 

F=J+I, J+I-1, --- J-I (or tol—J when J=/J), 

I=5/2, the nuclear spin of iodine, and 

eQ(d?V/dz*) represents the quadrupole coupling coeffi- 
cient, Q being the quadrupole moment of iodine, 
V, the molecular potential at the iodine nucleus, 
and zs, the coordinate axis along the C—I bond. 


Because of the nuclear effects the K-levels, associated 
with a given J-level, which otherwise would be degenerate 
for low J-values, are well separated. This allows two 
hyperfine patterns for the J=1 to J=2 transition corre- 
sponding to the transitions 0-0 and 1-+1 for K. The 
relative intensities of the lines of a given K-transition can 
be determined from the intensity rules for atomic fine 
structure. To determine the relative intensities of the two 
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Fic. 1. Theoretical and observed hyperfine structure of the CHy,l 
spectrum for the rotational transition J =1 to J =2. 











patterns arising from different X-transitions we have 
applied the theory of Dennison :* 


Intensity of terms for K divisible by 3 


Intensity of terms for K not divisible by 3 
_4P+4I1+3 


4P+4I 





» (2) 


where J is the nuclear spin of the 3 identical nuclei. Here J 
is } the spin of H; and the intensity ratio is 2:1. 

It is apparent from the above analysis that the micro- 
wave spectra of symmetric top molecules can be used in 
two different ways to ascertain nuclear spins. From for- 
mula (1) the spin of iodine is determined as 5/2. Formula 
(2) allows a determination of the spins of the 3 identical 
atoms on the corners of the tetrahedron. 

The quadrupole moment of the iodine nucleus is nega- 
tive, in agreement with the accompanying ICN measure- 
ments. The value of the coupling coefficient eQ(d*V/dz*), 
1520+15™°, for CHsI is different from the ICN value, 
2070+20 mc. Since Q is the same, this means that the 
0? V/d2* factor is significantly different for the two cases 
and that the C—I bonds in the two molecules are not 
equivalent. The difference in bonding is also revealed by 
the widely differing C—I interatomic distances, 2.00A for 
ICN and 2.13A for CHsl. 

The moment of inertia, Jz, for the ground vibrational 
state determined from the hypothetical frequency corre- 
sponding to no quadrupole-moment splitting of the energy 
levels is 111.4 10-*° g cm.? The C—I interatomic distance 
remains 1.13A, as determined in the previous paper. 


* The research described in this report was supported by Contract No. 
W-28-099-ac-125 + the Army Air Forces, Watson Laboratories, Air 


Materiel Comman 
1 Walter Gordy, A. G. Smith, and James W. Simmons, Phys. Rev. 
71, 917(L) (1947), i the preliminary report a rather large error was 
made | in estimated intensity of one of the lines. This is corrected in 
*For references to the development of the theory of quadru —_ 
moment interactions see the accompanying note on BrCN and 
ts + Coles and W. E. Good, Phys. Rev. 70, 979 (1946). 
Dennison, Rev. Mod. Phys. 3, 280 (1931); see also G. 
craberg, . Infrared and Raman — of Polyatomic Molecules, (D. Van 
Now Company, Inc., New York, 1945), p. 28. 





Proton-Proton Scattering at 10 Mev 


R. E, PEreRLS AND M. A. PRESTON 
The University, Birmingham, England 
June 19, 1947 


STUDY has been made of R. R. Wilson’s experi- 

mental results' in order to analyze the proton-proton 
interaction in the *P state. We have found that the 
P-scattering is compatible with the assumption that the 
nuclear force in this state is given by a 10-Mev repulsive 
potential with a range of 2.5 10-" cm. 

We have compared Wilson's results with the general 
formulae for scattering with arbitrary phases of the S- 
and P-waves. The best value for the S-phase is Kg= 52° 30’. 
Although the cross sections found for angles less than 30° 
in the center of mass system are not sufficiently accurate 
to be used for numerical calculations, they definitely 
show that the P-phase is negative, implying a repulsive 
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TABLE I. Relation between P-wave phases and interaction potentials, 











—=——S=S= 
p K 
Potential K, (calculated) . (perturbation 
—30 Mev 5 6° 2.7° 
—20 3.0 18 
—10 1.3 0.9 
0 0.0 0.0 
10 —0.8 ~0.9 
20 —1.3 eal 18 
30 —-18 —27 








force. The value K,= —0.8° fits the experimental points 
best, though with the stated experimental errors this 
figure is an order of magnitude estimate only. Since Wilson 
states that the absolute values of the cross section are less 
reliable than their ratios, we have confirmed that, ad- 
justing both Ko and Ki, a 10 percent change in scale 
changes K, only by about 0.2°. 

From the S-phase it is possible to calculate for each 
energy the “‘relative slope” of the wave function (ratio 
of the derivative to the function itself) at the limit of 
the range of the nuclear forces. The derivative of this 
quantity with respect to energy is directly related to the 
mean range of forces, without any hypothesis about the 
shape of the well, and combining the results of Wilson 
and Creutz at 8 Mev? with those at 10 Mev and with 
previous values at smaller energies, we have obtained for 
the range of forces the value 2.4, 10-" cm. 

Assuming this range, one can calculate the P-phase 
shift as a function of the depth of a rectangular potential 
well (added to the Coulomb potential). The results are 
given in the second column of Table I. 

Since the experimental value of K, is —0.8°, we see 
that the proton-proton interaction in the *P state is 
represented by a repulsive potential of about 10 Mev. The 
measurements of scattering at 14.5 Mev‘ do not disagree 
with this result. 

Since this conclusion is at variance with the theoretical 
curve given in the graph in Wilson’s note,’ which indi- 
cates a P-phase of about —2.1° for a 10.5-Mev potential, 
we have checked our calculation by means of perturbation 
theory using Eq. (8.1) in the paper of Breit, Condon, and 
Present.* This leads to the values in the third column of 
the table, which all seem to agree with the second column 
for small potential. 

This result is quite compatible with current ideas about 
nuclear forces. If we assume, as is customary, that the 
forces are charge-independent and are mainly of the 
“Majorana” and “Heisenberg” type, we can estimate the 
depth of well for these two forces from the *S and the 
virtual 'S level of the deuteron. With the value of the 
range used above, this leads to a repulsive potential of 
about 12 Mev for the *P state, in quite satisfactory agree- 
ment with the present data. However, it is evidently 
quite possible to admit a certain amount of “Wigner” or 
“Bartlett” force and some tensor force. 

1R. R. Wilson, Phys. Rev. 71, 384 (1947). 


2R. R. ho ay ne E. C. Creutz, ay Rev. 71, 339 (1947). 
2G. Breit, E. U. Condon, and R. D. Present, Phys. Rev. 50, 825 


CR R. Wilson, E. J. Lofgren, J. R. Richardson, B. T. Wright, and 
R. ‘S Shankland, Phys. Rev’ 71, 560 (1947). 
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Meson Production 


W. HoRNING AND M. WEINSTEIN 
Department of Physics, University of California, Berkeley, California 
June 9, 1947 


cMILLAN and Teller' have treated the production 
M of mesons by the collision of heavy particles with 
heavy nuclei. With the essential assumption that the 
collision time is small compared to nuclear periods, they 
have used the Fermi-Thomas model to calculate the effect 
of the presence of many nucleons in depressing the thresh- 
old. In obtaining the cross section near threshold, however, 
they have extrapolated from higher energies, and so over- 
estimated the effect because of the rapid departure of the 
cross section from a simple power law. On the same theo- 
retical basis, a better estimate can be obtained by carrying 
through the perturbation calculation near threshold. 

The process is third order in the meson field interaction. 
The heavy particles transfer momentum by exchanging a 
meson; the final meson is emitted before, during, or after 
the exchange. In the symmetric scalar theory the matrix 
elements do not vary appreciably near threshold, and the 
integral over final states reduces to the volume in phase 
space in which energy and momentum are conserved and 
the exclusion principle’ satisfied. The cross section is 


Cecalar = 0.7 (g?/he)* V(uc/h)(AE/uc?)*!? cm’, 


where V is the volume of the nucleus, u the mass of the 
meson, and AE the excess of incident energy over threshold. 
The result depends strongly on the coupling constant, 
which is usually taken to be (g?/hc)= ys. 

In the symmetric scalar theory there is considerable 
contribution to the integral from terms in which the final 
meson is produced during the exchange of a virtual meson. 
If the mechanism for heavy particle momentum transfer 
were a potential interaction, these terms would not be 
present. The importance of these terms arises directly from 
the finiteness of the ratio of the meson mass to the mass of 
the heavy particle. Also, using a potential interaction, the 
contribution from terms in which scattering precedes 
meson emission is, at threshold, equal and opposite in sign 
to that in which scattering follows meson emission. This 
equality is destroyed if exchange of positive, negative, and 
neutral mesons scatters the heavy particles. These proper- 
ties of the meson field distinguish it from the more familiar 
case of the electromagnetic field, for which a potential 
description of the interaction is valid. 

The energies now available on the Berkeley cyclotron 
make especially interesting the case of bombardment with 
neutrons obtained by stripping an energetically homo- 
geneous beam of deuterons. The energy distribution of a 
beam of neutrons so prepared is governed largely by the 
distribution in momentum of the neutron relative to the 
center of gravity of the deuteron and so can be calculated 
by Fourier analysis of the deuteron wave function assuming 
an exponential well. In the symmetric scalar theory, for 
velocities of the center of gravity of the deuteron between 
0.8 and 1 times the threshold velocity for the neutron, the 
total cross section varies from 0.014 to 2.1 10-* cm’. 
When the center of gravity moves with threshold velocity, 
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the cross section is that for neutrons at 16 Mev above 
threshold. 

These results, of course, share the large uncertainty in 
all present meson theories. 

We wish to thank Professor J. R. Oppenheimer for 
guidance in this problem. 


1W. G. McMillan and E. Teller, Phys. Rev. 72, 1 (1947). 





The Disintegration Scheme of Sc*** 


CHARLES PEACOCK AND ROGER G. WILKINSON 
Indiana University, Bloomington, Indiana 
June 21, 1947 


SMALL 180° beta-ray spectrometer, having good 

resolution, has been used to study the radiations from 
Sc**. Samples were produced by bombarding Sc2O; with the 
external cyclotron deuteron beam. For beta-ray measure- 
ments, sources were prepared by precipitating chemically 
purified SczO; on thin mica backing; for gamma-ray de- 
terminations, the photoelectrons ejected from a thin lead 
radiator were studied. Our results indicate that there are 
two groups of beta-rays of maximum energies 0.36+0.01 
Mev and 1.49+0.01 Mev, as well as two gamma-rays of 
energies 0.88+0.02 Mev and 1.12+0.02 Mev. Acomparison 
of the relative number of particles emitted in each beta- 
group shows that about 6 percent of the beta-rays belong to 
the higher energy group. A straight line Fermi plot is ob- 
tained for the 0.36-Mev group but not for the less intense, 
higher energy group, indicating a transition of the for- 
bidden type. 

Coincidence absorption experiments were carried out by 
Mr. Edward T. Jurney and Miss Margaret Ramsey to 
ascertain the mode of decay. Beta-gamma and gamma- 
gamma coincidences show conclusively that each beta-ray 
of the 0.36-Mev group is followed by the cascade emission 
of the two gamma-rays. Definite results have not yet been 
obtained from coincidence measurements for the very weak 
1.49-Mev group, but it is quite likely from energy con- 
siderations that a beta-ray emission of this group is followed 
by the emission of a 0.88-Mev gamma-ray. A study of the 
photoelectron lines bears this out to some extent. By using 
Gray’s' empirical curve for the variation of the photo- 
electric absorption coefficient with energy, it is found that 
the 0.88-Mev gamma is roughly 5 percent more abundant 
than the 1.12-Mev gamma. It is concluded, therefore, that 
the decay scheme is that shown in Fig. 1. This has been 
substantiated in part by Miller and Deutsch.* 

Further investigation reveals the presence of a weak 
conversion line at 0.86 Mev. This suggests strongly that the 
0.88-Mev gamma-ray is internally converted. Assuming 
this to be the case, one finds for the conversion coefficient 
a value of 0.0001. However, a careful check on the decay of 
the line, and further chemical separation are necessary to 
verify this conclusion. 

It is to be noted that the values of the energies given here 
differ somewhat from those of Walke*® and Meitner,* who 
have previously reported a 0.26-Mev beta-ray and a 
1.5-Mev gamma-ray. In addition, Walke suggested a 1.5- 
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Mev beta-ray and a 1.25-Mev gamma-ray. Meitner was 
unable to find a beta-group with a 1.5-Mev end pointand 
attributed Walke’s result to scattering. Careful study has 
convinced us that the 1.49-Mev group cannot be due to 
scattering inherent in the geometry of our instrument or to 
secondary electrons from the source or backing. Further 
confirmation of this group is found in the fact that it decays 
with the same half-life as the 0.36-Mev group. The decay 
has been followed for about one‘half-life and is found to be 
approximately 85 days as reported by Walke. 


* This research was supported by a grant from the Office of Naval 


Research. 
1L. H. Gray, Proc. Camb. Phil. Soc. 27, 103 (1931). 
2 A. Miller and M. Deutsch, paper presented at Montreal meeting of 
American Physical Society, abstract to be published in Physical Review. 
+H. Walke, Phys. Rev. 57, 163 (1940). 
4L. Meitner, Arkiv f. Mat. Astron., och Fysik, 32A, No. 6. 





The (4n+1) Radioactive Series: The Decay 
Products of U?* 


F. HAGEMANN, L. I. ae M. H. Strupier,* A. GHIoRSO,** AND - 


SEABORG** 
Argonne National Laboratory, Chicago,t Illinois 
June 19, 1947 


URING 1944-1946 we studied the chain of decay 

products of U™*, the new isotope of uranium which 
was first separated and examined (1941-1942) by Seaborg, 
Gofman, and Stoughton.'! These decay products, which 
constitute a substantial fraction of the entire missing, 
4n+1, radioactive series are listed, together with their 
radioactive properties, in Table I. 

The radioactivity of the Tl has not yet actually been 
observed, its existence in the chain is inferred from the 
partial alpha-decay of Bi**. The isotope Pb™ has been 
previously reported, as a result of its production by the 
(d, p),? (n, y),* and (n, p)* reactions. 

A number of the preceding members of this (4n+1) 
radioactive series have been previously reported as follows: 

B~ a a B~ 
94Pu**—>,,Am**!—+93 N p57 9) Pa*8®—+ U8, 
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Taare I. 
Type of Bneray of 

an A (Mev) 
wTh™ a 7108 yr, Te 
seRa® A 14.8 day: ~0.2 
Ace ee 10.0 day: 5.80 
siFr™t a 4.8 min. 6.30 
ssAtl? a 0.018 sec, 7.00 
s3 Bi2!3 &: ao} 47 min. ~13 ey 
ssPo%8 a very short 8.30 
(si T 12) (8~) ? ? 
s2Pb20 B- 3.3 hr. 0.7 
s3Bi2? stable 

== 


The 27.4-day Pa® was first reported by Meitner, Strass. 
mann, and Hahn,‘ and the doubts as to this isotopic assign- 
ment which later arose as a result of the discovery of 
fission were cleared up by the work of v. Grosse, Booth 
and Dunning* and Seaborg, Gofman, and Kennedy.* The 
2.25 10° year Np** was first identified by Wahl and 
Seaborg,’ while Pu™ and 500-year Am™! were first reported 
by Seaborg, James, and Morgan.® Also of interest are two 
previously reported beta-emitting radioactive isotopes, 
23-minute Th*** from the reaction Th* (n, y) and 7-day 
U4" from the reaction U** (m,2n), which may be 
referred to as “‘collateral’’ members of the series. 

As the name for the (4n+1) radioactive decay family 
we suggest ‘‘neptunium”’ series or family ; thus the longest- 
lived member would give its name to the family in a manner 
similar to the naming of the uranium and thorium decay 
series. 

Another independent study of the decay products of U™ 
was carried on simultaneously by A. C. English, T. E. 
Cranshaw, P. Demers, J. A. Harvey, E. P. Hincks, J. V. 
Jelley, and A. N. May of the Division of Atomic Energy 
of the National Research Council of Canada, which has 
resulted in essentially similar findings.” 


* Present address, Institute of Nuclear Studies, University of Chicago, 
Chicago, Illinois. 

** Present address, Radiation Laboratory, University of California, 
Berkeley, California. 

t Work performed under auspices of Manhattan District, Contract 
Number W-7401-eng-37. 

1G. T. Seaborg, J. W. Gofman, and R. W. Stoughton, Phys. Rev. 71, 
378 (1947). 

2R. L. Thornton and J. M. Cork, Phys. Rev. 51, 383 (1937). 

*W. Maurer and W. Ramm, Zeits. f. Physik 119, 602 (1942). 
( 4L. Meitner, F. Strassmann, and O. Hahn, Zeits. f. Physik 109, 538 

1938). 
( Das v. Grosse, E. T. Booth, and J. R. Dunning, Phys. Rev. 59, 322 
1941). 

$ i T. Seaborg, J. W. Gofman, and R. W. Stoughton, Phys. Rev. 59, 
321 (1941). 

7A. C. Wahl and G. T. Seaborg, reported in Chem. Eng. News 23, 
2190 (1945). 

8G. T. Seaborg, R. A. James, and L. O. Morgan, reported in Science 
104, 379 (1946) and Chem. Eng. News 25, 358 (1947). 

®* E. Fermi, E. Amaldi, O. D'Agostino, F. Rasetti, and E. Segré, Proc. 
Roy. Soc. A146, 483 (1934). 

10 Y, Nishina, T. Yasaki, H. Ezoe, K. Kimura, and M. Ikawa, Phys. 
Rev. 57, 1182 (1940). 

ut E, M. McMillan, Phys. Rev. 58, 178 (1940). 

2 English, Cranshaw, Demers, Harvey, Hincks, Jelley, and May, 
Phys. Rev. 72, 253 (1947). 
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LETTERS TO 


The (4n+1) Radioactive Series* 


sH, T. E. CRaNSHAW, P. Demers, J. A. Harv EY, 
A. C. = . Hincks, J. V. JELLEY, AND A. N. May 


fo Energy, National Research Council of Canada, 
piinen e Chalk River, Ontario, Canada 


June 20, 1947 


HE three radioactive series known to occur naturally 
T are characterized by members whose mass numbers 
are represented by the integers 4n (thorium series), (4n+2) 
(uranium series), and (4n+3) (actinium series). W hile the 
existence of a fourth series comprising the ‘‘missing’’ 
radioactive (4 +1)-isotopes has been postulated at various 
times,'~* early searches which were made for its members 
in natural minerals gave negative results.** More recently 
Wahl’ has reported the observation of a line at mass 237 
in mass spectrograms from certain minerals, and attributes 
it toa member of the new series. Since 1935 a- and 8-active 
(4n+1)-species have been produced by artificial transmu- 
tation, and among the heavy elements almost a dozen 
have been reported in the literature. (See for example 
Seaborg’s Table of Isotopes* and his recent review.’) In 
1942 Seaborg, Gofman, and Stoughton*® discovered and 
investigated the properties of U™, the long-lived a-emitting 
product of the 8-decay of Pa** and the third member of a 
(4n+1)-chain which starts with Th*. 

The investigation of the ensuing route taken by the chain 
which was undertaken by the authors started with a search 
for the decay products of milligram amounts of U™®. 
A chemical extraction of radium resulted in an a-activity 
which decayed in a complex manner with periods of the 
order of days. Direct evidence for a chain including at 
least four a-emitters was first provided by the observation 
of 4-track stars in a photographic emulsion which had 
been soaked in the radium solution. Subsequent measure- 
ment of decay periods and a-particle energies showed 
that these formed part of a new disintegration series, 
distinct from the thorium, the uranium, and the actinium 
series. Making use of both chemical and physical evidence 
it has been possible to deduce the main route of the series, 
which, after the a-decay of U**, passes through seven new 
nuclear species, and ends with the known’ 3.3-hr. 8-decay 
of Pb*®* to stable Bi?®. 

The path of the series which passes through isotopes of 
the type (4n+1) is shown in Fig. 1. Both U* and its 
daughter Th”® are fairly long-lived, and, as far as this 
work was concerned, it has been convenient to treat the 
former as the parent of the series. Since the radium isotope 
is B-active, the main chain does not include an emanation 
followed by the usual “A” and “B” products. Instead, the 
main a-sequence starts with Ac™5, and the series thus 
includes isotopes of the recently named elements, francium 
(Z =87) and astatine (Z = 85). The end product is Bi?** and 
not a lead isotope as in the case of the other three series. 
(It is interesting to note that the route is almost exactly 
that predicted by Turner‘ in 1940.) A 2 percent a-branching 
of Bi** was detected, but the disintegration of TI? could 
not be investigated because of the low activity. A search 
was made for a-branching of Ra®* but none could be found. 
it is probably less than one percent. 

The disintegration properties of the members of the 
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Fic. 1. The path of the (4n +1)-series following U™, 


series are summarized in Table I. Previously published 
data included for completeness are shown in light-faced 
type, while results obtained by the authors are in heavy 
type. The a-particle energies were measured with a grid- 
type ionization chamber, biased amplifier, and pulse 
analyzer. The 21-millisecond half-life of At®*!? was measured 
by photographing pulses on an oscilloscope with a triggered 
time base, and the 3.2-microsecond half-life of Po*®*® by a 
coincidence experiment using a variable resolving time. 

The a-emitters of the series are found to lie on a good 
straight line in a Geiger-Nuttall plot (logarithm of dis- 
integration constant against logarithm of a-particle energy ), 
with the exception of Po*® which falls well below the line 
in the manner of the very short-lived C’ bodies. This line 
lies between and is parallel to those representing the 4n 
and (4n+-3)-series. 

The authors feel that the name ‘‘Neptunium Series,"’ 
after the longest-lived known (4n+1)-member, Np*’,’ 
should be adopted. 

We wish to thank those members of the Montreal 
Laboratory who gave us assistance during the above 
investigation, and especially Professor F. A. Paneth and 
Dr. B. L. Goldschmidt whose interést and advice were 
invaluable. 


TABLE I. The (4n +1) radioactive series. 








Energy of particles 





in Mev 

Type of 
Species radiation Half-life a 8 
2J™ ale 1.63 X10 y® 4.825 +-0.003 
ooTh™ a {5 X108 y®) 

{ S10ty | ~S 

ssRa™™ 6 14 <0.05 
spAc™ a 104 5.801 +0.010 
a7Fr™1 a 6.31 +0.02 
ssAt?!7 a 2.1X107 5s 7.023 +0.010 
Bit (2%), B(98%) 46 m 5.86 40.03 ~1.3 
ssP0%3 a 3.2 X10-' s 8.336 +0.005 
a2Pb™ feud 3.3 h® 0.70° 
ss Bi™® stable 
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A detailed report of the work will be published later. 

Note added in proof: A recalculation of the results from 
the measurement of the half-life of polonium 213 gives a 
better value of 4.4 10~* seconds. 


* The work here was carried out in the Montreal Laboratory 
of the Division of Atomic Ayo Research Council of Canada, 
between November 1944 and 

1A. ‘ Russell, Phil. Mag. 46, r% 91923). 

2 W. P. Widdowson and A. S. Russell, Phil. Beg. , 293 (1924). 
crag cures H. von Halban, ‘and P. Preiswerk, J e phys. et rad. 6, 361 

4L. A. Turner, Phys. Rev. 57, 950 (1940). 

5 W. Wahl, —~y ng) 16 (1941). 

*G. T. Seaborg, Rev. Mod. Phys. 16, 1 (1944). 

7G. T. Seaborg, Science 104, 379 (1946). 

378 (947). J. W. Gofman, and R. W. Stoughton, Phys. Rev. 71, 
®*R. L. Thornton and J. M. Cork, ty, Rev. 51, 383 (1937); K 

i z= V Phys B Rev. 60, 619 (1941). 

ioactivity and Nuclear Physics (D. Van Nostrand 

Company, q ‘New York, 1947), information attributed to G. T. Sea- 





Thresholds for Creation of Particles 
H. FESHBACH 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
AND 
L. I. ScuirF* 
University of Pennsylvania, Philadelphia, Pennsylvania 
June 9, 1947 
N the design of very high energy particle accelerators, 
it is of interest to know the smallest energy that a 
bombarding particle need have in order that a given com- 
bination of particles can be created in a collision. The 
calculation of such a threshold on the basis of purely 
energetic considerations, under the assumption that a 
mechanism of creation exists, was pointed out by Dr. M. S. 
Livingston to be a straightforward application of rela- 
tivistic dynamics. This note presents the pertinent for- 
mulas, and discusses some limitations on their application. 
We consider a collision in which a particle of rest mass, 
m,, and kinetic energy, £, collides with an initially sta- 
tionary particle of rest mass, m2, and new particles of 
aggregate rest mass M (in addition to m, and mz) are 
created. E is smallest when all particles are at rest after 
the collision in the center of mass-coordinate system (in 
which the center of mass of m and mz is at rest before the 
collision). Application to this problem of the relativistic 
relation between energy and momentum leads to an ex- 
pression for the threshold kinetic energy: 
E= Me(m+m1+4M)/m2. (1) 
Equation (1) can be solved for Mc*, which is the energy 
that is available in the center of mass system for particle 
creation, reactions, etc.: 
Me? =(2mac?E + (mi +-m2)°ct}!— (mi +ma)c. 


In the non-relativistic limit (E<mc* and mac*), this re- _ 


duces to the usual expression m2E/(m:+ m2). 

If the colliding particles are protons (m;=m:=m) and 
a pair of neutrons or a positive and negative proton are 
created (M=2m), Eq. (1) shows that the threshold energy 
is 6 m?®=25.6 Bev. If an electron or a photon (m,=0) 
collides with a proton, the threshold for nucleon-pair 
production is 4 mc*=23.7 Bev. The corresponding thresholds 
for meson-pair production are higher than would be ex- 
pected from non-relativistic theory by about 5 percent in 
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the case of protons and about 10 percent in the case of 
electrons or photons. 

Equation (1) indicates that a considerable reduction ; in 
the thresholds for nucleon-pair production would occur if 
the target particle were a heavy nucleus rather than a 
proton. This is actually true only if the target nucleus 
remains intact after the collision. We now show that the 
transfer of a large amount of momentum to the nucleus 
is very likely to break it up. The recoil momentum q of 
mz is given by: 


(q/cP= 2M(mi+4M)(m2+4M)(m, tmt4M) 
(m+m2+M/) 


In the case of a heavy nucleus (m2Km, m:<M), q~mc 
for nucleon production either by protons or electrons. The 
region in which the creation process takes place has dimen- 
sions of the order of k/mc~10-" cm. This is so small that 
only one of the constituent nucleons of the target nucleus 
is likely to be involved. Thus the nucleus remains intact 
only if this nucleon transfers the momentum g to the 
remainder of the nucleus. 

The fraction of the events in which this occurs is of the 
order of the fraction of all collisions between a nucleon of 
momentum g, and a nucleus in which the nucleon is 
scattered at a large angle and the nucleus absorbs most of 
the momentum. Such a scattering process can be treated 
by the Born approximation. The cross section for back 
scattering is roughly equal to 

(m?/h*)| fe*** V(r)dr|?, where k~(q/h), (2) 
and the total scattering cross section is approximately the 
geometrical nuclear cross section R*, where R is the nuclear 
radius; here, V(r) is the interaction potential between 
nucleon and nucleus. Equation (2) is roughly equal to 
(mV /h®)*(h/q)*~(1/Ro')(h/g)*, where V is the average 
nuclear potential, and Rp is the range of nuclear forces. 
We are thus led to expect that the target nucleus remains 
intact in about (A/qRo)*(k/gR)*~10~7 of those collisions in 
which nucleon creation occurs, and is energetically possible 
with a proton target (the corresponding fraction for 
meson creation is considerably larger). This indicates that 
the effect is extremely improbable, and especially difficult 
to observe since newly created neutrons (although perhaps 
not negative protons) would tend to be masked by frag- 
ments of nuclei broken up by ordinary collisions. 

McMillan and Teller! have pointed out that the kinetic 
energy of the nucleons inside a nucleus lowers the threshold. 
If this energy is assumed to be not greater than 30 Mev, 
it can be shown that the thresholds for nucleon-pair 
production are reduced to 4.1 Bev for protons and 2.6 Bev 
for electrons or photons. Since very few of the nucleons 
have the maximum kinetic energy, the cross section will 
be small just above the threshold. A simple estimate based 
on the Thomas-Fermi model indicates that the effective 
threshold is half-way between the figures given here and 
those given in the third paragraph. However, the difficulty 
of observing the effect in the presence of a background of 
nuclear fragments might make it worth while to work with 
a hydrogen target at the higher thresholds. 


* Now at Stanford University 
1W. G. McMillan or E. Teller, Phys. Rev. 72, 155(A) (1947). 
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The Mean Life of Mesotrons at an 
Altitude of 11,500 Ft. 


Harowtp K. TicuHo 
Department of Physics, U niversity of Chicago, Chicago, Illinois 
June 23, 1947 


N 1943 Nereson and Rossi' described an experiment for 
determining the mean life of cosmic-ray mesotrons at 
sea-level in various solid absorbers. The purpose of this 
letter is to report on a measurement of the mean life which 
was carried out at Climax, Colorado (11,500 ft. above sea- 
level) to determine whether the mean life is altitude- 
dependent. The method consisted of a direct comparison of 
the lifetime of each individual mesotron, stopping in a 
10-cm Al absorber with the period of vibration of an 
accurately calibrated quartz crystal.* Figure 1 shows the 
counter-tube assembly, and the integral-decay curve ob- 
tained by this method. The upper set of points between 0.05 
and 0.5 ysec. represents the total number of delays observed 
with the Al absorber in its place, whereas the lower set was 
obtained by subtracting all delays occuring with the Al 
absorber removed. When delays of this kind (produced by 
natural delays in the counters accompanying a simultane- 
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DECAY CURVE at CLIMAX,COLO. 
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ous event) are subtracted, the curve remains a true 
exponential to time delays as small as 0.05 ysec. This fact 
indicates that if any contribution to the decay curve of 
delays arising from the disintegration of of short-life 
(ro<10~* sec.) mesotrons should exist, it must be quite 
negligible for time intervals larger than 0.05 sec. 

The data in Fig. 1 were analyzed according to the 
statistical method of Peierls? which permits an evaluation 
of the mean life and its standard deviation from a set of 
individual lifetime determinations. The value obtained by 
this method for the mean life of mesotrons at rest is 
ro=1.78+%0.10 ysec. For the calculation of the standard 
deviation the points between 0.05 and 0.5 usec. were 
omitted because of the natural delays which occur oc- 
casionally in the G. M. counters in that interval. Assuming 
that the mean range of the disintegration electrons, having 
an energy of 50 Mev, is 7 cm in Al, the ratio of disintegra- 
tion electrons to mesotrons stopped in the absorber is 
calculated to be 0.50+0.06. 

The value for the mean life determined at Climax is not 
in complete agreement with that obtained by Nereson and 
Rossi' at sea-level. These authors compiled their best value 
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Fic. 1. The counter-tube assembly and the integral-disintegration curve obtained at Climax, Colorado. The upper set of points in the time interval 


from 0.05 to 0.5 usec. is due to counter-pulse delays accompanying simultaneous events. 
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for the mean life ro=2.15+0.07 usec. by adding all their 
data collected in various materials. It seems difficult to 
explain the difference between the two values by a system- 
atic error in the experimental method since the disintegra- 
tion curve appears linear, and hence the error would have 
to increase linearly with the time interval which has actually 
been recorded. A shortening of the mean life, however, may 
be due to a contribution by negative mesotrons decaying in 
the absorber. The following letter explains why such a 
contribution might result in an apparent change of the 
mean life. 

The author wishes to express hf gratitude to Dr. Walter 
O. Roberts of the Fremont Pass Station of the Harvard 
College Observatory at Climax, and to the Climax Molyb- 
denum Company for making available the facilities required 
for carrying out this investigation. 

1N. Nereson and B. Rossi, Phys. Rev. 64, 199 (1943). 


2H. Ticho, Rev. Sci. Inst. 18, 271 (1947). 
3R. Peierls, Proc. Roy. Soc. 149, 467 (1935). 





Errata: The Double Focusing 
Beta-Ray Spectrometer 


[Phys. Rev. 71, 681 (1947)] 
FRANKLIN B. SHULL AND Davip M. DENNISON 
Randall Laboratory of Physics, University of Michigan, 
Ann Arbor, Michigan 


R. EDWIN M. McMILLAN has kindly pointed out 

an error in our article which he discovered by com- 
paring our equations with the results of some earlier 
unpublished calculations of his own. The error has its 
origin in the expression for H, given on page 682. The last 
term in this equation should read [—(8—a/2)z*Ho/a*]. 
The subsequent calculations are, we believe, all correct, 
but this correction introduces changes in certain of the 
coefficients in the equations. The final result, the expression 
for r* given in Eq. (25), must be modified as follows. The 
coefficient of (6s)? should be [(48—3)/3a] rather than 
[(48—2)/3a]. The coefficient of ¢,* should be [(168/3—2)a] 
rather than [(168/3—2/3)a]. The remaining terms, as well 
as Eq. (26), are correct as they stand. 

The new expression for r* is less favorable, since the ¢, 
defocusing may be eliminated for 8=} but not the ¢, 
defocusing. To eliminate the latter, 8 must equal j.The 
correct choice of 8 will depend upon the baffle system to 
be employed. It may often be more convenient to allow 
a wider variation in ¢, than in ¢, in which case 8 should 
be }. Although the focused image will not be as perfect as 
that shown in Fig. 1, the conclusion still stands that, with 
the double focusing spectrometer, the image may be made 
both more intense and also sharper than with the usual 
semicircular spectrometer. 

Another advantage pointed out to us by Dr. McMillan 
is that the dispersion (pér/rép where p is the electron mo- 
mentum) is twice as great as in the semicircular case. 

Since submitting our paper we have received a reprint 
from Dr. N. Svartholm! in which the image formed by a 
point source is discussed. His results are in agreement with 
our corrected equations. 


1N. Svartholm, Ark. f. Mat. Astron. och Fys. 33A [24] (1946). 
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On the Magnetic Exchange Moment 
for H*® and He’ 
FeLtx VILLARS 


Swiss Federal Institute of Technology, Zurich, Switserland 
June 23, 1947 


ECENT experiments of Bloch and others! on the 
magnetic moment of H# gave the value 1.0666 for the 
ratio of the magnetic moments of H’ and proton. Assuming 
a value of 2.789 n.m. for u(p), we find a(H*)=2.975 =y(p) 
+0.186 n.m. This result seemed to be in contradiction to 
the results of theoretical investigations of Sachs and 
Schwinger,? which require 4(H*)< u(p), unless very arti- 
ficial assumptions on the *P and ‘P admixtures to the 
S-component of the ground-state eigenfunction are made. 
However, Schwinger’s ansatz‘ for the nuclear. Hamil. 
tonian does not allow for taking into account the charge- 
exchange phenomena connected with the interaction of 
nucleons. On the other hand, it is well known that these 
phenomena give rise to exchange moments.’ Whereas the 
magnetic exchange dipole moment vanishes in the case of 
the deuteron (on account of the symmetry properties of the 
de-eigenfunction), this is not the case for, H*® and He’, 
provided that the quantum number T of the total isotopic 
spin is 4. (It vanishes for T= }.) 

A calculation has been carried out on the basis of the 
symmetrical pseudoscalar meson theory. According to this 
theory, the H* and He? ground states are doublet states 
both with respect to spin and isotopic spin (S= 7'=4) and 
symmetrical with respect to permutations of the space 
coordinates of the particles, if we neglect the influence of 
the tensor force. The latter is responsible for small admix- 
tures of higher states, the influence of which may be 
neglected here, since there is a non-vanishing expectation 
value of the exchange moment in the above described 
S state. 

The exchange moment operator is given by M=M 
+M®; 


ef? 
MOoO= -£, _ (rA X Bef 548(sAl oA xo8) / 
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ran (1+ . -)— (oxo) bem, 


MT AB 


2 
M®= Paci > (r4Xr®)3- (24x28). V(AB), 
2° 4<B 
where V(AB) is the interaction energy of the nucleon pair 
AB in the pseudoscalar theory: 


V(AB)= { }(o4e®) + (3(04248)(o8248)/r45°— (0408) 


x“ (+ 1 + 1 = 
3 wrap (uras)/) rap 
On account of its symmetry properties, the expectation 


value of M® vanishes, whereas M“ gives, in units of 
nuclear magnetons 


1 
M® = (8/3)y(fu)®-Ts dve(——— 2)e-oas, 


(y is the ratio of the masses of proton and meson: y=10, 
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T:+} is the charge of the nucleus, ¢ the orbital part of the 
ground state eigenfunction, 248 =s4—2, rap=|s4*|, (fu) 
the dimensionless coupling constant of the meson field: 
(fu):2%5, and »=Mc/h.) A rough evaluation of M“ has 
been made with the help of gauss functions g~ exp(—ar’), 
f=h(rit tri? t+rss*), with the following result (J is the 


volume integral in 4). 


ple 1.0 1.5 2.0 2.5 0.75 
J -0.14 —0.21 -023 -0.23  —0.058 
M 40.18 40.28 +4031 +031  +0.077 


Thus, with reasonable values of , (fu), and u?/a we obtain 
both the right sign and right order of magnitude of the 


correction to be added. 

It should be noted that for He’ the correction is equal in 
magnitude but opposite in sign. We would, therefore, 
expect for He* a total magnetic moment py(N)—M 
22.1 n.m. Experimental evidence would be very 


interesting. 


1F, Bloch, A. C. Graves, M. Packard, and R. W. Spence, Phys. Rev. 

71, 373 and 551 (1947); H. L. Anderson and A. Novick, Phys. Rev. 71, 
1947). 

mt G. Sachs and J. Schwinger, Phys. Rev. 70, 41 (1946). 

2R. G. Sachs, Phys. Rev. 71, 457 (1947). 

+E. Gerjuoy and J. Schwinger, Phys. Rev. 61, 138 (1942). 

8S. T. Maand F. C. Yu, Phys. Rev. 62, 118 (1942); C. Moller and L. 
Rosenfeld, Kungl. Danske Vidensk. Selsk. 20, No. 12 (1943); W. Pauli 
and S. Kusaka, Phys. Rev. 63, 400 (1943). 





Errata: Theory of Dipole Interaction in Crystals 


[Phys. Rev. 70, 954 (1946)] 
J. M. LuttTinGer AND Laszo TIsza 
Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 


EVERAL misprints have been noticed in the above 
paper. These are the following: 
P. 956, line 7 should read p,”, py’, p2”, v=1, ---, 8. 
P. 956, Eq. 7 should read 
8 
U=- ts 2 , 2 Fy px" py’. (7) 
val zy 


B, vl 2 
P. 957, Eq. 12 should read 
Z,=(—)mtbiktyils fFo4, .--, 8, (12) 


P. 960, last equation. The denominator should be raised 
to the 5/2 power. : 
P. 960, Table II, first line should read. 


fre= —4[5S,(0, 3, 4)—S,(4, 0, 0)). 


P. 960. The small table under Table II contains several 
inversions and a sign error. It is correctly given by: 


S(t} 0 0)=—15.040 S,( }4 3)=31.521 
S0 4 4)= 4.334 S,(h 4 4)= 2.599 
Sy(4 4 = 10.620 S(O 34 4)=12.329 


Lastly, in Table V, p. 963, lines 4 and 5 should be 
exchanged (which moves a minus sign down one line), and 
line 12 should read —2X3—Yg+Zs. 

The authors would like to thank Professor L. W. 
McKeehan for having pointed out several of the above 
misprints. 
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Burst Production by Penetrating 
Cosmic-Ray Particles* 


HerBerRT BripGe, Bruno Rossi, AND ROBERT WILLIAMS 


Laboratory for Nuclear Science and Engineering, Massachusetts Instituic 
of Technology, Cambridge, Massachusetts 


May 23, 1947 


TRAY of Geiger-Mueller tubes, G, and an ionization 
chanrber, C, were arranged, respectively, above and 
below a lead block 6 in. thick, as shown in Fig. 1. The 
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FG. 1. Schematic arrangement of equipment. 


Geiger-Mueller tubes were connected in parallel. Each was 
1 in. in diameter and 20 in. long. The chamber was cylin- 
drical in shape, 3 in. in diameter, 20 in. long, and was filled 
to 7.3 atmospheres with highly purified argon so that “‘fast"’ 
electron pulses would be recorded quantitatively. The 
pulses of the ion chamber were applied to the vertical 
deflecting plates of a cathode-ray oscilloscope through a 
linear amplifier and a delay line. The oscilloscope was 
provided with a fast horizontal sweep (5 microseconds per 
inch) which was triggered by the coincidences between the 
(undelayed ) pulses of the ionization chamber and the pulses 
of the Geiger-Mueller tubes. The oscilloscope screen was 
photographed on a moving film. The individual counting 
rates of the chamber (N,) and of the Geiger-Mueller tubes 
(N,) were also recorded. 

A polonium source of a-particles was placed on the wall 
of the chamber for the purpose of calibration. The resolving 
time (r:+72) for the selection of coincident pulses was 
determined both by direct observation of the pulses on the 
oscilloscope screen and by counting chance coincidences 
between pulses in the Geiger-Mueller tubes and a-particle 
pulses in the ionization chamber. Its value was found to be 
50 microseconds. 

For the main experiments, the circuits were adjusted so 
as to record only pulses greater than 1.1 times a Po 








258 LETTERS TO 


a-particle pulse. (The a-particle pulse corresponds to about 
35 fast particles incident on the chamber from the vertical 
direction.) Measurements were taken at sea level in a light 
wooden building, and at 30,000 feet in the rear cabin of a 
B-29. The results are as shown in Table I. 


TABLE I. Results of measurements. 











Time of Ne (Counts Ny, (Counts  Coinci- 

obser- in the intheG.M. dence 

Altitude vation chamber, counter tray, (per (per hour) 
(feet) (hours) per hour) per hour) hour) NeN¢(ri+172) 


0 292 1.94 0.72 X105 0.25 0.002 
30,000 1.6 500 14.7 105 96 10 


Chance 
coincidence 











At both elevations only a small fraction of the coincident 
records obtained is accounted for by chance coincidences. 
We also believe-that air showers have a negligible effect on 
our results because of the small number of showers with a 
particle density high enough to be recorded in these 
measurements, and thus conclude that most of the coinci- 
dences observed are caused by ionizing particles capable of 
producing a burst in the ionization chamber after traversing 
6 in. of lead. 

It does not seem possible to assume that these particles 
are electrons since an electron capable of producing a suffi- 
ciently large shower under 6 in. of lead must have an energy 
in excess of 10" ev. Hence they must be of the “‘penetrating”’ 
type. On the other hand, the total number of penetrating 
cosmic-ray particles only increases by a factor of about 6 
from sea level to 30,000 feet, while the observed effect 
increases by a factor of several hundreds. We are thus led to 
the conclusion that the penetrating component at 30,000 
feet contains particles which are much more effective in 
producing bursts than are ordinary mesons. 

As has been shown by one of us,' the shape of the pulse 
from an ionization chamber enables one to decide in most 
cases whether the pulse is caused by heavily ionizing par- 
ticles from a disintegration or by an electron shower. 
Classification of our records according to pulse shape gives 
the results shown in Table II. 


TABLE II. Classification of records according to pulse shape. 








Percent of pulses from 





Altitude Heavily ionizing Uncertain 
(feet) Showers particles origin 
0 78 _ 22 
30,000 57 30 13 








It is likely that most of the showers observed at sea level 
are initiated by collision or radiation processes of ordinary 
mesons. These processes can only account for a very small 
fraction of the showers observed at 30,000 feet. It is natural 
to assume that the majority of the showers at this altitude 
are produced by electrons or photons generated by primary 
cosmic-ray particles (protons?) either directly or through 
the intermediary of short-lived mesons. 

‘At 30,000 feet our data show evidence for nuclear 
disintegrations produced by penetrating ionizing rays. It is 
likely that most of the particles responsible for this effect 
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are also primary “protons” (even though some of them may 
be negative mesons undergoing nuclear capture after bein 
brought to rest). . 

The shower production by primary cosmic-ray particles 
for which evidence is found in the present experiments, may 
account for that part of the soft component which cannot 
be explained by the disintegration of ordinary mesons.? 


* The research described in this letter was i 
Contract NS ORL-78U; 5. Navy Deariment, ice of Naval Rt 
» e B- u or the mea t i i 
provided by the U. S. Army Air Force. “ish altitude was 
1H. B . Phys. Rev. 72, 172(A) (1947). 
2H. Bridge and B. Rossi, Phys. Rev. 71, 379 (1947). 





Erratum: The Band Theory of Graphite 


[Phys. Rev. 71, 622 (1947)] 


P. R. WALLACE 
National Research Council of Canada, Chalk River Laboratory, 
Chalk River, Ontario 


Figure 13, the plot of the trend of pL/pll (ratio of 
resistivity across graphite“planes to that¥in the planes), 





° 100 200 300 400 
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Fic. 1. 


appeared in the paper in incorrect form. The corrected 
graph is as follows: 





Nuclear Moments of the Bromine Isotopes* 


S. B. Bropy, W. A. NIERENBERG, AND N. F. RAMSEY 
Columbia University, New York, New York 
June 13, 1947 


HE molecular beam! resonance method has been used 

to study the radiofrequency spectrum? of Br’ and 

Br* in CsBr and LiBr. Typical results are shown in Fig. 1. 
The positions of the resonance minima are the same in 
CsBr and LiBr, indicating that the minima are due to the 
bromine isotopes. The positions of the resonance minima 
shift proportionately to the magnetic field. This is the usual 
criterion! for assuming that the observed resonance fre- 
quency is at the Larmor frequency of the nucleus in the 
external magnetic field, which gives for the gyromagnetic 
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REDUCTION OF BEAM INTENSITY 
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Fic. 1. Reduction in beam intensity as a function of 
oscillator frequency. 


ratio 
g1=hv/unH, (1) 
where uw is the nuclear magneton. 

As will be discussed ina later paper by one of us (N.F.R.), 
there are certain special circumstances in the case of very 
large nuclear electrical quadrupole-moment interactions 
when the frequency of an observed resonance minimum 
may be double that of Eq. (1). This possibility is excluded 
in the present experiments by comparison with optical fine- 
structure experiments of Tolansky’ interpreted by Schmidt* 
which indicate a value of 2.6 nuclear magnetons for 
bromine in fair agreement with our observed values using 
Eq. (1) and in marked disagreement if values of the fre- 
quency double that given in (1) are used. Furthermore, the 
integrated intensity of the resonance minima in Fig. 1 is 
greater than would be expected for a resonance of double 
the frequency of Eq. (1) and is approximately correct for a 
normal resonance. 

With (1) being used to determine the gyromagnetic 
ratios and with the spins of the bromine isotopes taken to 
be*® 3, the magnetic moments of the bromine isotopes are 
2.110+0.021 and 2.271+0.023 nuclear magnetons. Since 
each isotope is present in about 50 percent abundance, it is 
impossible from the present experiments to assign specific 
isotopes to the different observed moments. 

Although the resonance minima are quite broad they do 
not possess the characteristic shape for a nuclear electrical 
quadrupole interaction with spin** §. This is somewhat 
surprising in view of the large quadrupole interaction of Br 
in BrCN found by Townes.® The spectra correspond more 
nearly to those of an arbitrarily assumed I-J interaction 
which has been found useful in accounting for some other 
resonance-minima curves. Further studies of the shape and 
interpretation of the bromine resonance are being made. 

* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 
os twee S. Millman, P. Kusch, and J. R. Zacharias, Phys. Rev. 55, 
wiW,Ac Nirenberg N. F. Ramsey, and S. B. Brody, Phys. Rev. 70, 

*Tolansky, Proc. Roy. Soc. A136, 585 (1932). 

‘T. Schmidt, Zeits. f. Physik 108, 408 (1938). 

'C. H. Townes, A. N. Holden, J. Bardeen, and F. R. Merritt, Phys. 


Rev. 71, 644 (1947). 
*B. T. Feld and W. E. Lamb, Jr., Phys. Rev. 67, 15 (1945). 
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Millimeter-Wave Spectra: Hyperfine Structure 
of BrCN and ICN* 
WaLtTer Gorpy, WILLIAM V. Sattu, A. G. SMITH, AND HAROLD RING** 


Department of Physics, Duke University, Durham, North Carolina 
June 7, 1947 


SING the method prevously described,! we have 
examined the J=3 to J=4 transition of BrCN and 
the J=4 to J=5 transition of ICN, both of which occur 
in the region of 9 mm. As was discovered by Townes? in 
the J=2 to J=3 transition of BrCN, Br has a rather 
large quadrupole moment which splits the pure rotational 
lines into a number of components. In our laboratory the 
quadrupole moment of the iodine nucleus has been ob- 
served through the splitting of the rotational lines of 
methyl iodide’ as well as through the hyperfine structure of 
ICN reported here. 

The quadrupole-moment effects were first noticed in 
atomic spectra, and the interaction was shown by Schuler 
and Schmidt‘ to be proportional to the square of the cosine 
of the angle between the angular momentum vector of the 
orbital electron and the nuclear spin vector. Casimer,* 
and also Bethe* have derived a theory for the quadrupole 
effects in atomic spectra. This theory as well as the in- 
tensity rules derived for atomic fine structure’ has been 
adapted to the analysis of the molecular hyperfine struc- 
ture reported here. 

In Figs. 1 and 2 are shown the calculated and observed 
hyperfine spectra for BrCN and ICN. The calculated lines 
are given for the ground vibrational states only, although 
in both these molecules lines caused by excited vibrational 
states were observed. These are indicated in the experi- 
mental diagrams by lighter lines. The lines on the low 
frequency side are caused by the vibration which corre- 
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Fic. 1. Theoretical and observed hyperfine structure of the BrCN 


spectrum for the rotational transition J =3 to J =4. The heavy lines are 
for the ground vibrational state; the light ones in the observed spectrum 
are due to excited vibrational states. 
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Fic. 2. Theoretical and observed hyperfine structure of the ICN 
spectrum for the rotational transition J =4 to J =5. The heavy lines are 
for the ground vibrational state; the light ones in the observed spectrum 
are due to excited vibrational states. 


sponds to the stretching of the carbon-halogen bond, and 
in each case the stronger lines in the ground state are 
repeated with identical spacing in the excited vibrational 
state. From the separation of these groups it is possible to 
calculate the increase in the average bond length in the 
excited state over that for the ground state. They are: 
for BrC 0.0029A (0.17 percent increase); for IC 0.0037A 
(0.18 percent increase). This increase in length is due to 
the anharmonicity of the bond-stretching potential func- 
tion. In the bending vibrational state the molecule is on 
the average shorter and the moment of inertia smaller. 
Hence, these lines fall on the high frequency side. Because 
of the interaction of the two modes of bending (/-type 
doubling) the hyperfine pattern of the ground state is not 
repeated exactly, as it is for the stretching vibration. 

The moment of inertia and the structure of the BrCN 
molecule have been determined by Townes? from the 
rotational lines occurring at the lower frequencies. Our 
- values for this molecule agree with his. The moment of 
inertia of ICN in its ground vibrational state is deter- 
mined from our data as 260X10~* g cm*. Assuming the 
CN distance to be that observed for BrCN and CICN 
(1.15A), the CI distance is 2.00A. No measurements from 
other sources are available for comparison, but this value 
should be accurate to 1 percent at least. To determine the 
structure of this molecule completely, we plan to search 
for the carbon isotope which should be detectable by our 
method. 

The spin of the iodine nucleus is definitely established 
as 5/2. The quadrupole coupling factors, eQd?V/d2, are: 
for Br7*, 670415 Mc; for Br*', 567415 Mc; for I’, 
2070+20 Mc. The quadrupole moment of I is negative, 
indicating that the iodine nucleus is flattened in the 
direction of the spin axis, i.e., it is oblate. That of Br is 
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positive, indicating that the Br nucleus is elongated j 
the direction of the axis of spin, i.e., it is prolate. Whee 
the coupling coefficient @V?/d=* is evaluated, the quad- 
rupole moments of I and of Br can be determined from 
this data. We are now attempting to evaluate 8 V /azt 
from data of other types. . 

We wish to acknowledge a helpful discussion with Dr. Cc 
H. Townes, of the Bell Telephone Laboratories. ; 


* This research was supported by Contract No. W- 
with the Army Air Forces, Watson Laboratories, Air Maree 
wont and by a grant-in-aid from the Research Corporation of New 

ork. 

** Frederick Gardner Cottrell Fellow. 
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On the Inelastic Scattering of Neutrons 
in Crystals* 


Orro HALPERN 
Columbia University, New York, New York 
June 20, 1947 


HE inelastic scattering of a monochromatic neutron 
beam by crystal powders has been studied exten- 
sively by C. G. Shull and E. O. Wollan'; the experimental 
technique followed closely the well-known procedure in 
x-ray experiments. The very remarkable results obtained 
by these authors with a variety of scatterers seem to be in 
quantitative, and even qualitative disagreement with the 
presently accepted theories. It is, therefore, perhaps of 
interest to indicate an experimental approach to this 
problem which has no analogy in x-ray technique and 
promises to give valuable independent information. We 
suggest the study of the transmission of monochromatic 
neutron beams through single crystals of iron in a trans- 
verse magnetic field. Theory* then predicts very precisely 
that (apart from special angles of incidence here left out of 
consideration) no coherent scattering will occur. The 
weakening of the beam is due to incoherent nuclear scat- 
tering, inelastic lattice scattering, isotope disorder effects, 
and capture. Of these effects only inelastic lattice scattering 
will be influenced by a magnetic field. This special polariza- 
tion effect has its origin in the fact that the incoherence of 
the inelastic scattering is due to energy exchange with the 
lattice while no element of incoherence enters into the 
interaction with the individual atom. Here the amplitudes 
of the wave scattered by the nucleusand by the magnetically 
active electron shell, superpose in the same way as we know 
them to interfere in the well-studied polarization effects 
obtained with polycrystalline iron. 

In the present experiment the difference in transmission 
in the unmagnetized and saturated state of the scatterer is 
measured. This difference, AJ/J, will again be given in the 
customary approximation by the well-known formula 


Al/1=}p"@. (1) 


In (1) p’ denotes the excess cross section arising from the 
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interference term of nuclear and magnetic scattering in the 
inelastic process ; d is the thickness of the crystal. The 
corresponding quantity p for coherent scattering has been 
assumed to be, and most likely is, much larger than p’; for 
this reason an arrangement eliminating coherent effects is 
here described. : | 

Any observed “single transmission effect”’ is, therefore, 
determined by the inelastic scattering and can in turn be 
used to obtain quantitative information about its magni- 
tude. One advantage of the arrangement here described lies 
in the replacement of a scattering by a transmission experi- 
ment. It is furthermore not difficult to obtain good satura- 
tion of the single crystal provided the magnetizing field is 
parallel to an easy axis of magnetization. Thus depolariza- 
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tion effects* can here be eliminated more easily than in 
previous experiments. 

The presently accepted theoretical estimates would lead 
to moderate but easily measurable effects; if, on the other 
hand, as suggested by the experiments quoted,' inelastic 
scattering should be much larger than expected, the 
“inelastic single-transmission effect’’ would become quite 
large. The isolation of the inelastic scattering would permit 
one to test the theory‘ rather precisely. 

* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 
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